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Wo present a measurement of time-dependent CP- violating asymmetries in neutral B meson decays 
collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at the Stanford 
Linear Accelerator Center. The data sample consists of 29.7 fb~^ recorded at the T(4S') resonance 
and 3.9 fb~^ off-resonance. One of the neutral B mesons, which are produced in pairs at the T(4S'), 
is fully reconstructed in the CP decay modes J/ipK°s, ip{2S)K°s, XciK°s, J/ipK''^ {K'° Jf^Tr") 
and J/ipK^, or in flavor-oigenstatc modes involving /?'•*•' 7r/p/ai and J/tp K*'^ K^tv^). The 

flavor of the other neutral B meson is tagged at the time of its decay, mainly with the charge of 
identihed leptons and kaons. A neural network tagging algorithm is used to recover events without 
a clear lepton or kaon tag. The proper time elapsed between the decays is determined by measuring 
the distance between the decay vertices. Wrong-tag probabilities, the time-difference resolution 
function, and the B^-B^ oscillation frequency Anid are measured with a sample of about 6350 
fully-reconstructed decays in hadronic flavor-eigenstate modes. A maximum-likelihood fit to 
this flavor eigenstate sample finds Am^ = 0.516 ± 0.016 (stat) ± 0.010 (syst) ps^^. The value of 
the asymmetry amplitude sin2/3 is determined from a simultaneous maximum-likelihood fit to the 
time-difference distribution of the flavor-eigenstate sample and about 642 tagged B^ decays in the 
CP-eigenstate modes. We find sin2/3 = 0.59 ± 0.14 (stat) ± 0.05 (syst), demonstrating that CP 
violation exists in the neutral B meson system. We also determine the value of the CP violation 
parameter |A| — 0.93 ± 0.09 (stat) ± 0.03 (syst), which is consistent with the expectation of |A| = 1 
for no direct CP violation. 

PACS numbers; 13.25.Hw, 12.15.Hh, ll.30.Er 
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I. INTRODUCTION 

CP violation has been a central concern of parti- 
cle physics since its discovery in 1964 0]. Interest 
was heightened by Sakharov's observation in 1967 
that without CP violation, a universe that began as 
matter-anti-matter symmetric could not have evolved 
into the asymmetric one we now see. An elegant ex- 
planation of the CP- violating effects in decays is pro- 
vided by the CP-violating phase of the three-generation 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix However, existing studies of CP violation in 
neutral kaon decays and the resulting experimental con- 
straints on the parameters of the CKM matrix [Q do not 
provide a stringent test of whether the CKM phase de- 
scribes CP violation [||. Moreover, the Standard Model 
does not, through the CKM phase, incorporate enough 
CP violation to explain the current matter-anti-matter 
asymmetry Q . Understanding CP violation thus remains 
a pressing challenge. 

An excellent testing ground for CP violation is pro- 
vided by B mesons through particle-anti-particle mixing. 
A particle that is purely at time t = will oscillate 
between that state and with a frequency Am^, the 
difference between the masses of the two neutral B mass 
eigenstates. If decays to a CP eigenstate / are observed, 
any difference between the rates when starting with a 
or with a is a manifestation of CP violation. In 
some circumstances, including those in the experiment 
described here, the fundamental parameters of CP vio- 
lation in the CKM model can be measured from such 
time-dependent rate asymmetries, unobscured by strong 
interactions. For example, a state initially produced as a 
B'^ (P*^) can decay to J/ipK'^ directly or can oscillate into 
a B° {B°) and then decay to J/i{jK°. With little theoret- 
ical uncertainty in the Standard Model, the phase differ- 
ence between these amplitudes is equal to twice the angle 
/3 = arg [ -V^dV;jVtaVtl ] of the Unitarity Triangle. The 
CP- violating asymmetry can thus provide a crucial test 
of the Standard Model. 

The unitarity of the three-generation CKM matrix can 
be expressed in geometric form by six triangles of equal 
area in the complex plane. A nonzero area @] directly 
implies the existence of a CP-violating CKM phase. The 
most experimentally accessible of the unitarity relations, 
involving the two smallest elements of the CKM matrix, 
Vub and Vtd, has come to be known as the Unitarity Tri- 
angle. Because the lengths of the sides of the Unitarity 
Triangle are comparable, the angles can be large, lead- 
ing to potentially large CP-violating asymmetries from 
relative phases between CKM matrix elements. 

In e+e~ storage rings operating at the T(45') reso- 
nance, a B'^B^ pair produced in an T{4S) decay evolves 
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in a coherent P-wave state. If one of the B mesons, re- 
ferred to as Ptag, can be ascertained to decay to a state 
of known flavor, i.e. B^ or i?'^, at a certain time itag, 
the other B, referred to as iJrocj «^ that time must be of 
the opposite flavor as a consequence of Bose symmetry. 
Consequently, the oscillatory probabilities for observing 
B°B^, B°B° and B°B° pairs produced in T{4:S) decays 
are a function of At = tree — ^tag, allowing mixing fre- 
quencies and CP asymmetries to be determined if At is 
known. The charges of identified leptons and kaons are 
the primary indicators of the flavor of the tagging B, but 
other particles also carry flavor information that can be 
exploited with a neural network algorithm. The recon- 
structed neutral B is found either in a flavor eigenstate 
(Picc — -Bfiav) or a CP mode (Proc = Bcp) by full recon- 
struction of its observed long-lived daughters. 

At the PEP-II asymmetric e+e^ collider resonant 
production of the T(4S') provides a copious source of 
B^B'^ pairs moving along the beam axis {z direction) with 
an average Lorentz boost of {(3^) = 0.55. Therefore, the 
proper decay-time difference At is, to an excellent ap- 
proximation, proportional to the distance Az between 
the two iJ^-decay vertices along the axis of the boost. 
At K, Az/c{/3j). The average separation between the 
two B decay vertices is Az = {/3j) ctb — 260 /im, while 
the RMS Az resolution of the detector is about 180 /im. 



A. Measurement of flavor oscillations 

The phenomenon of particle-anti-particle mixing in 
the neutral B meson system was first observed almost 
fifteen years ago 10 . The oscillation frequency in P°- 
B*^ mixing has been extensively studied with both time- 
integrated and time-dependent techniques |]ll] . By inter- 
changing bd with bd, B'^-B'^ mixing changes the additive 
bottom quantum number by two units, i.e., \AB\ — 2. 
In the Standard Model, such a process is the result of 
ordinary |Ai3| = 1 weak interactions in second order in- 
volving the exchange of charge-2/3 quarks, with the top 
quark contributing the dominant amplitude. A measure- 
ment of Anid is therefore sensitive to the value of the 
CKM matrix element Vtd- At present the sensitivity to 
Vtd is not limited by experimental precision on Am^, but 
by theoretical uncertainties in the calculation, in partic- 
ular the quantity J^Bb, where fs is the B^ decay con- 
stant, and Bb is the so-called bag factor, representing 
the AP = 2 strong- interaction matrix element. There 
may also be contributions from interactions outside the 
Standard Model 

Beyond these questions of fundamental interest, since 
the measurement of Anid incorporates all elements of the 
analysis for time-dependent CP asymmetries, including 
B reconstruction, tagging, and At determination and res- 
olution, it is a essential test of our understanding of these 
aspects of the sin2/3 measurement. 

For the measurement of Anid, one neutral B (Pflav) 
is fully reconstructed in a flavor eigenstate as 
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FIG. 1: Expected distribution for mixed and unmixed 
events a) with perfect tagging and resolution, and b) with 
typical mistag rates and At resolution. 



D^*^- Tx^ I p-^ I aX or Jl^K*^ {K*° K+n-), while the 
second is tagged by its decay products. The probability 
for B'^-B'^ mixing is a function of Ato^ and the proper 
time difference At between the two B decays: 

Prob(S°B" ^ or B°B°,B°B°) = 

£e-r|^*l(lTcosAmrfAt), (1) 

where tbo = l/F is the B" lifetime. The observed B°B" 
system produced in an T(45) decay can be classified 
as mixed or unmixed depending on whether the recon- 



the opposite flavor as the tagging B = Stag- If the At 
resolution and flavor tagging were perfect, the asymme- 
try as a function of At 



,(At) 



jV„nmix(At) - A^mix(At) 
7Vu„^i,(At) + 7V„,i,(At) 



(2) 



would describe a cosine function with unit amplitude. 
The asymmetry goes through zero near 2.1 B'^ proper 
lifetimes and the sensitivity to Am^, which is propor- 
tional to At^e~^l^*l sin^ Am^At, reaches a maximum in 
this region. If the tagging algorithm incorrectly identi- 
fies the tag with a probability w, the amplitude of the 
oscillation is reduced by a dilution factor T) — {1 ~ 2w). 
When more than one type of flavor tag is used, each has 
its own mistag rate wi. 

Neglecting any background contributions, the proba- 
bility density functions (PDFs) for the mixed (— ) and 
unmixed events, 7i±, can be expressed as the convo- 
lution of the underlying oscillatory physics distribution 



/i±(At;r, ATO<i,w) = 4 e"^!'^*! [1 ± Pcos AwdAt] (3) 

with a time-difference resolution function TZ{5^_ = At — 
Attt.uc;a) to give 



Ti± (At; r, Amd, w, a) = 

ft.±(Attruo; r, Amd, w) (g) 7^((5t; a). 



(4) 



structed flavor-eigenstate B^c 



has the same or 



where At and Attruo are the measured and the true time 
differences, and a are parameters of the resolution func- 
tion. Figure ^ illustrates the impact of typical mistag and 
At resolution effects on the At distributions for mixed 
and unmixed events. 

A full likelihood function is then constructed by sum- 
ming 'H± over all mixed and unmixed events in a given 
uniquely-assigned tagging category i and over all tagging 
categories 



In Cn 



tagging 

E 



. unmixed 



(5) 



This can be maximized to extract the mistag fractions 
Wi and resolution parameters ai and, simultaneously, the 
mixing rate Am^. The correlation between Wi and Amd 
is small, because the rate of mixed events near At = 0, 
where the B^-B'^ mixing probability is small, is princi- 
pally governed by the mistag rate. Conversely, the sensi- 
tivity to increases at larger values of At; when At 
is approximately twice the B lifetime, half of the neutral 
B mesons will have oscillated. 



B. Measurement of CP asymmetries 



For the measurement of CP asymmetries, one B {Bcp) 
is fully reconstructed in a CP eigenstate with eigenvalue 
VCP = -1 ^{2S)Kl or xaK^s) or +1 (J/^if°), 

while the second is tagged with its decay products just as 
for the mixing measurement. The Bqp sample is further 
enlarged by including the mode J/ipK*'^ {K*^ K^;!^^)- 
However, due to the presence of even [L = 0, 2) and odd 
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(L = 1) orbital angular momenta in the J/tpK*^ system, 
there are rjcp = +1 and —1 contributions to its decay 
rate, respectively. When the angular information in the 
decay is ignored, the measured CP asymmetry in J/ip K*^ 
is reduced by a dilution factor D± — 1~2R±, where R± is 
the fraction of the L — 1 component. We have measured 
R± = 0.160 ± 0.032 ± 0.014 ||l3| which, after acceptance 
corrections, leads to an effective r]cp = +0.65 ± 0.07 for 
the J/i;K*° mode. 

The expected time evolution for the tagged Be p sam- 
ple depends both on B'^-B^ mixing and on the decay am- 
plitudes of B'^ and B^ to the final state / through a single 
complex parameter A. Mixing generates a lifetime differ- 
ence as well as a mass difference between the two neutral 
B meson mass eigenstates, but the lifetime difference is 
expected to be small since it is a consequence of common 



final states in B^ and B^ decays. Such common states, 
which include the CP eigenstates studies here, make up 
a very small fraction of the decay width because they are 
quite suppressed by CKM matrix elements. Dropping 
these, and thus ignoring any lifetime difference, results 
in a simple expression for A in terms of the | AS ~ 1| and 
IAS = 21 interactions, 



mH^B=2m{f\n^B=i\B°) 

(S°|Hab=2|S") {f\-H^B=i\B^)' 



(6) 



Redefining the states for B^ and B^ by multiplying them 
by two different phases has no effect on A, which is thus 
phase-convention independent, as every physical observ- 
able must be. The decay distributions are 



J 



^r|At| 



2/mA 



■ sin Am^ At - 



1-|A|^ 
1 + |A|^ 



cos Amd At 



(7) 



r 



where the + or — sign indicates whether the Stag is 
tagged as a i?° or a S", respectively. The dilution factor 
T) — 1~ 2w accounts for the probability w that the flavor 
of the tagging B is identified incorrectly. 

The distributions are much simpler when |A| — 1, 
which is the expectation of the Standard Model for de- 
cays like B^ J/tpK'^. If all the mechanisms that 
contribute to the decay have the same weak phase then 
the ratio of the weak decay amplitudes in Eq. ^ is just 
rycpe^''^'^™ , where 0dcc is the weak phase for B" /; 0dcc 
is convention dependent and unobservable. The remain- 
ing factor introduces a phase due to B'^-B'^ mixing. The 
combination of these phases is convention independent 
and observable. 

For decays such as B^ J/ipK'^, or more generally 
{cc)K^ and {cc)Kl, an explicit representation for A can 
be found from the ratio of the amplitude for {cc)K^ 
to the interfering process B^ ^ B^ ^ {cc)K^ {cc)K'''. 
The decay B^ — > {cc)K'^ involves a 6 ^ ccs transition 
with an amplitude proportional to [V^JJ^F^g], while 
{cc)K^ provides analogously a factor VcpiVcb^cs]- Be- 
cause B'^ —> B^ mixing is dominated by the loop diagram 
with a t quark, it introduces a factor [V^^Vj^/F^^Vj^], 
while ~-> mixing, being dominated by the c-quark 
loop, adds a factor [VcdKs/KdKsl- Altogether, for tran- 
sitions of the type h — > ccs. 



VCP 



VCP 



Vcpe 



-2i/3 



V jV* 

* ch cs 



v*,v 
V .,v* 



The time-dependent rate for decay of the Bcp final state 
is then given by 

/±(At;r,ATOd,w,sin2/3) = 

_g-r|At| _p ^^^p sin 2/3 sin Am^At] . (9) 

In the limit of perfect determination of the flavor of the 
fully-reconstructed B in the Sflav sample, which we as- 
sume throughout, the dilution here and in the mixed and 
unmixed amplitudes of Eq. |^ arise solely from the -Stag 
side, allowing the values of the mistag fractions Wi to 
be determined by studying the time-dependent rate of 
B^-B^ oscillations. 

To account for the finite resolution of the detector, the 
time-dependent distributions f± for i?'^ and B^ tagged 
events (Eq. |^) must be convolved with a time resolution 
function TZ{5t — At — Attme', a) as described above for 
mixing, 

T± {At; r, Am^, w, sin 2/3, a) = 

/±(Attruo; r, Amrf, w, sin 2/3) (g> n{St;a), (10) 

where d represents the set of parameters that describe 
the resolution function. In practice, events are separated 
into the same tagging categories as in mixing, each of 
which has a different mistag fraction Wi, determined in- 
dividually for each category. Figure ^ illustrates the im- 
pact of typical mistag and At resolution effects on the 
At distributions for B'^- and B"-tagged CP events. 
It is possible to construct a CP- violating observable 



(8) 



^cp(At) = 



J^+(At) - JP_(At) 
T+{At)+T_{At)'- 



(11) 
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-5 5 

At (ps) 

FIG. 2: Expected At distribution for B°- and B°-tagged CP 
events a) with perfect tagging and At resolution, and b) with 
typical mistag rates and At resolution. 



which, neglecting resolution effects, is proportional to 
sin2/?: 



ylcp(At) cx -?7cp25sin2/3sin Am^Ai. (12) 



Since no time-integrated CP asymmetry effect is ex- 
pected, an analysis of the time-dependent asymmetry 
is necessary. The interference between the two ampli- 
tudes, and hence the CP asymmetry, is maximal af- 
ter approximately 2.1 B'^ proper lifetimes, when the 
mixing asymmetry goes through zero. However, the 
maximum sensitivity to sin2/3, which is proportional to 
e-r|At| gjj^2 ^jyi^^f^ occurs in the region of 1.4 lifetimes. 

The value of the free parameter sin2/3 can be extracted 
with the tagged Bcp sample by maximizing the likeli- 
hood function 



CP 



tagging 

E 



In J^+(At; r, Am^, a, Wi, sin2/3) 

B° tag 



E 

3" tag 



In JF_(At; F, Am^, a, Wi, sin 2/3) 



(13) 



where the outer summation is over tagging categories i 
and the inner summations are over the B'^ and B^ tags 
within a given uniquely-assigned tagging category. In 
practice, the fit for sin2/? is performed on the combined 
flavor-eigenstate and CP samples with a likelihood con- 
structed from the sum of Eq. I and |l|, in order to de- 
termine sin2/3, the mistag fraction Wi for each tagging 
category, and the vertex resolution parameters a^. Addi- 
tional terms are included in the likelihood to account for 
backgrounds and their time dependence. 

The mistag rates can also be extracted with a time- 
integrated analysis as a cross check. Neglecting possi- 
ble background contributions and assuming the flavor of 
Bfiav is correctly identified, the observed time-integrated 
fraction of mixed events Xobs can be expressed as a func- 
tion of the B'^-B^ mixing probability Xd'- 



Xobs = Xd + (1 - 2xd) w, 



(14) 



where Xd = ^xy{l + xj) = 0.174 ± 0.009 |Tl| and 
Xd = Anid/T. Taking advantage of the available decay 
time information, the statistical precision on w can be 
improved by selecting only events that fall into an opti- 
mized time interval \At\ < to, where to is chosen so that 



the integrated number of mixed and unmixed events are 
equal outside this range. With the use of such an op- 
timized At interval the time-integrated method achieves 
nearly the same statistical precision for the mistag rates 
as a full time-dependent likelihood fit. 



C. Overview of the analysis 

This article provides a detailed description of our pub- 
lished measurement of flavor oscillations |15|1 and CP- 



violating asymmetry 
These measurements 



I in the neutral B meson system, 
lave six main components: 



• Selection of the Bcp sample of signal events for 
neutral B decays to CP modes J/i^K^, i}{2S)K^, 
XciK°, J/i^K*^ {K*^ ^sTi-O), and J/^jKI; se- 
lection of the Bflav sample of signal events for 
neutral flavor-eigenstate decays to D^*'>~'k'^ / p+/a^ 
and J/i}K*° {K*° ^ if+Tr-); _selection of the B+ 
control sample in the modes D^*^°tt'^ , J/ip K^*'^'^ , 
■(/'(25)i4r+, XciK^; and selection of a semileptonic 
neutral B sample in the mode D*^£~V, as de- 



scribed in Section [II 
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Determination of the flavor of the -Stag, as de- 



scribed in Section IV 



as discussed in Section VF, while the latter results in 



• Measurement of the distance Az between the two 
B° decay vertices along the T{AS) boost axis, and 
its conversion to At, as described in Section 

• Construction of a log-likelihood function to de- 
scribe the time evolution of signal and background 
events in the presence of mixing and CP asymme- 
tries, as described in Section VI; 



• Measurement of the mixing rate Am^, mistag frac- 
tions Wi, and vertex resolution parameters di for 
the different tagging categories i, with an unbinned 
maximum-likelihood fit to the i?flav sample, as de- 
scribed in Section VII; 



• Extraction of a value of sin2/3, or more generally 
/mA/|A| and |A|, from the amplitude of the CP 
asymmetry, the mistag fractions Wi , and the vertex 
resolution parameters di for the different tagging 
categories i, with an unbinned maximum- likelihood 
fit to the combined -Bflav and Bcp samples, as de- 
scribed in Section VIII. 



Whenever possible, we determine time and mass reso- 
lutions, efficiencies and mistag fractions from the data. 
The measurement of Am^ is performed with a slightly re- 
duced subset of the full i?flav sample, which is optimized 
for such a precision measurement. The Bcp sample is 
not included, since this would add additional assump- 
tions about the resolution function without significantly 
improving the precision of Am^. The measurement of 
sin2/3 is performed with the full iJflav and Bcp samples, 
with a fixed value for Arud and the 5° lifetime. This 
strategy allows us to account correctly for the small cor- 
relations among the mistag rates. At resolutions param- 
eters, and sin2/3. The same Bflav sample and vertex sep- 
aration techniques have been used to determine precision 
values for the charged and neutral B lifetimes llj. 



substantially improved yields in Run 2 for reconstructed 
B mesons. 



A. The BABAR detector 

The BABAR detector is a charged and neutral spec- 
trometer with large solid-angle coverage. For this anal- 
ysis, the most important detector capabilities include 
charged-particle tracking, vertex reconstruction, and par- 
ticle identification. Charged particles are detected and 
their momenta measured by a combination of a 40-layer, 
small-celled drift chamber (DCH) filled with a 80:20 he- 
liumiisobutane gas mixture, and a five-layer silicon ver- 
tex tracker (SVT), consisting of 340 AC-coupled double- 
sided silicon microstrip sensors. The cells of the DCH 
are organized into 10 superlayers within which the sense 
wires all have the same orientation, thereby allowing 
segment-based tracking. Both the DCH and the SVT 
lie inside a 1.5-T solenoidal magnetic field. Beyond the 
outer radius of the DCH is a detector of internally re- 
fiected Cherenkov radiation (DIRC), which is used pri- 
marily for charged-hadron identification. The device con- 
sists of 144 fused silica quartz bars in which relativistic 
charged particles above the Cherenkov threshold radiate 
photons while traversing the material. The light is trans- 
ported by total internal reflection down the length of the 
bars to an array of 10752 photomultiplier tubes mounted 
on the rear of the detector, where the opening angle of the 
Cherenkov ring is measured. A flnely segmented electro- 
magnetic calorimeter (EMC), consisting of 6580 CsI(Tl) 
crystals, is used to detect photons and neutral hadrons, 
and also to identify electrons. The EMC is surrounded by 
a thin cylindrical superconducting coil and a segmented 
iron flux return, organized into a hexagonal barrel and 
two endcaps. The instrumented flux return (IFR) con- 
sists of multiple layers of resistive plate chambers (RPCs) 
interleaved with the flux-return iron and is used in the 
identification of muons and neutral hadrons. 



II. THE BABAR DETECTOR AND DATA SETS 

The data used in this analysis were recorded with the 
BABAR detector Q at the PEP-II colhdcr g in the pe- 
riod October 1999-June 2001. The total integrated lumi- 
nosity of the data set is equivalent to 29.7 fb^^ collected 
near the T(4S') resonance and 3.9 Ad"^ collected 40MeV 
below the T(4S') resonance (off-resonance data). The cor- 
responding number of produced BB pairs is estimated to 
be about 32 million. The T{AS) sample is sometimes di- 
vided into two subsamplcs for comparison purposes: data 
recorded in 1999-2000, about 20.7 fb^^ and referred to as 
"Run 1", and data recorded in 2001, about 9.0 fb^^ and 
referred to as "Run 2" . These subsamples differ primar- 
ily in the quality of the tracking system alignment and on 
the track-finding efficiency. The former requires a sepa- 
rate treatment of the At resolution for the two periods, 



B. Charged particle reconstruction 

Charged track finding starts with pattern recognition 
in the DCH, based on three different algorithms. The 
first uses the same fast algorithm employed by the Level-3 
trigger for finding and linking superlayer-based track seg- 
ments from moderate-to-high px tracks originating from 
the interaction point. Two subsequent track finders then 
work on superlayer segments not already attached to a 
reconstructed track. They are designed to find tracks 
with lower px, passing through fewer than the full ten 
superlayers of the chamber, or originating away from the 
interaction point. At the end of this process, all tracks 
are refit with a Kalman-filter fitter that takes into 
account the detailed distribution of material in the de- 
tector and the non-uniformities in the detector magnetic 
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field. These tracks are then projected into the SVT, and 
siUcon-strip hits are added if they are consistent within 
the extrapolation errors through the intervening mate- 
rial and field. A search is performed for tracks that are 
reconstructed with the remaining unused silicon clusters, 
again with two different algorithms. At the end of the 
SVT-only track finding, an attempt is made to match 
SVT- and DCH-only track segments, which may result 
when a hard scatter occurs in the support tube material 
between the two devices. 

Charged-particle transverse momenta px are deter- 
mined with a resolution parameterized by ct(pt)Ipt = 
0.0013(pT/GeV/c) + 0.0045. The SVT, with typical 
single-hit resolution of 10 /im, provides vertex informa- 
tion in both the transverse plane and in z, as well as 
the decay angles at the interaction point. Decay vertices 
for B meson candidates are typically reconstructed with 
a resolution of 50 /xm in z for fully reconstructed modes 
and about 100 to 150 /im for the vertex of the (unrecon- 
structed) tagging B meson in the event. The efficiency 
for finding tracks in hadronic events that traverse the full 
DCH radius {pr > 200 MeV/c) is about 90% for Run 1 
and 95% for Run 2. 



C. Neutral reconstruction 

EMC clusters are formed around initial seed crystals 
containing at least lOMeV of deposited energy. Neigh- 
boring crystals are added to the cluster if their energy 
exceeds 1 MeV. If the newly added crystal has energy 
greater than 3 MeV, its contiguous neighbors (including 
corners) are also considered for inclusion in the cluster. 
In order to identify cases where several showers are in 
close proximity, such as unresolved photons from high- 
energy 7r° decays, local maxima within a cluster are iden- 
tified. These local maxima are defined as candidate crys- 
tals that have an energy exceeding each of its neighbors 
by a fraction that depends on the number of crystals in 
the local neighborhood. Clusters are then divided into 
as many "bumps" as there are local maxima. The divi- 
sion is based on a two-dimensional weighting scheme that 
assumes electromagnetic shower shapes to divide up the 
cluster energy. The position of each bump is calculated 
with a logarithmic weighting of crystal energies. 

We determine whether a bump is associated with a 
charged or neutral particle by projecting all tracks in the 
event to the inner face of the calorimeter. A bump is de- 
termined to be neutral, and therefore a photon candidate, 
if no track intersects any of its crystals. A track inter- 
section is determined by computing the two-dimensional 
distance on the face of the calorimeter from the projected 
track impact point to the bump centroid. A require- 
ment is made on the difference between the measured 
intersection distance and the Monte Carlo expectation 
for different particle species based on the measured track 
parameters. 

The energy resolution in the EMC is measured directly 



with a radioactive source at low energy under ideal low- 
background conditions and with electrons from Bhabha 
scattering at high energy, from which we determine 
aiE)/E = (5.0 ±0.8)% at 6.13 MeV and (1.9 ±0.07)% at 
7.5 GeV. The energy resolution can also be extracted 
from the observed mass resolutions for vr" and i] de- 
cays to two photons, which are measured to be around 
7 MeV and 16 MeV, respectively. A fit to the observed 
resolutions obtained from the 7r°, r], and Bhabha sam- 
ples gives a photon energy resolution parameterized by 
aiE)/E = 0.023{E/ GeVy^^^ © 0.019. 

D. Particle identification 

Identification of electrons, muons and kaons is an es- 
sential ingredient in both B reconstruction and flavor 
tagging. Particle species can be distinguished by mea- 
surements of the specific energy loss {dE / dx) in the SVT 
layers and in the DCH gas along the particle trajectory, 
the number of Cherenkov photons and the Cherenkov 
angle in the DIRC, the electromagnetic shower energy 
in the EMC, and the particle penetration length in the 
IFR. Selection criteria are based on these quantities, on 
likelihood ratios derived from them, or on neural network 
algorithms combining different detector likelihoods. Typ- 
ically, looser selection criteria are applied for B recon- 
struction than for _B-fiavor tagging. Efficiencies and par- 
ticle misidentification probabilities are determined from 
data control samples with similar characteristics. 

1. Electron identification 

Electron candidates are identified primarily by the ra- 
tio of the bump energy in the electromagnetic calorime- 
ter to the track momentum, E /p. They also must have 
a measured mean AE/Ax in the DCH that is consistent 
with the electron hypothesis. In addition, for some ap- 
plications, the lateral and azimuthal shape of the EMC 
shower |2^ and the consistency of the observed and 
expected Cherenkov angle in the DIRC are used for iden- 
tification. Four different categories of electron candidates 
(VeryLoose, Loose, Tight, and VeryTight) are defined 
with the criteria listed in Table ^ Candidates that are 
not matched to an EMC bump are retained as noCal 
electron candidates if their measured dE / da; satisfies the 
same requirements as the VeryTight selection. Elec- 
tron identification efficiencies in the momentum range 
0.5 < p < 3.0GeV/c vary between 88% and 98% for the 
criteria in Table |, while the pion misidentification rates 
are below 0.3% for the VeryTight selection. 

2. Muon identification 

Muon candidates are primarily identified by the mea- 
sured number of hadronic interaction lengths n\ tra- 
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TABLE I: Criteria used for selecting the available categories of electron candidates. The difference between the measured 
mean dE/dx and the expectation for an electron is required to lie within the interval specified in terms of the expected dE/dx 
resolution a. 



Category 


dE/dx 


E/p 


Cumulative additional requirements 


VeryLoose 
Loose 
Tight 
VeryTight 


[-3a, 7a] 
[-3a, 7a] 
[-3a, 7a] 
[-2.2(7, 4(t] 


> 0.50 

> 0.65 
[0.75, 1.3] 
[0.89, 1.2] 


Lateral shower shape 

Azimuthal shower shape; consistency 

of DIRC Cherenkov angle {3a) 



TABLE II: Criteria used for selecting available categories of muon candidates. 



Category 


n\ 




nuits 




'2 1 
Xtrk/ '^layers 


•2 1 
Xfit/IT-layers 


Eemc [GeV] 


VeryLoose 


> 2.0 


< 2.5 


< 10 


< 6 






< 0.5 


Loose 


> 2.0 


< 2.0 


< 10 


< 6 


< 7 


< 4 


< 0.5 


Tight 


> 2.2 


< 1.0 


< 8 


< 4 


< 5 


< 3 


[0.05, 0.4] 


VeryTight 


> 2.2 


< 0.8 


< 8 


< 4 


< 5 


< 3 


[0.05, 0.4] 



versed from the outside radius of the DCH through the 
IFR iron, and the difference AnA between nx and the 
predicted penetration depth for a muon of the same mo- 
mentum and angle. Contamination from hadronic show- 
ers is rejected by a combination of the average number 
rihits and the variance (7n^._^^ of hits per RPC layer, the 
for the geometric match between the track extrapolation 
into the IFR and the RPC hits, Xtrk' ^^'^ of a poly- 

nomial fit to the RPC hits, Xfu- addition, for those 
muons within the acceptance of the EMC, we require the 
calorimeter bump energy E to be consistent with a mini- 
mum ionizing particle. Four different categories of muon 
candidates (VeryLoose, Loose, Tight, and VeryTight) 
are selected with the criteria listed in Table ||. In the 
forward region, which suffers from some machine back- 
ground, additional requirements are made on the fraction 
of RPC layers with hits. Muon identification efficiencies 
in the momentum range 1.1 < p < 3.0GeV/c vary be- 
tween 60% to 92% for the criteria in Table ^ while pion 
misidentification rates are about 3% for the Tight selec- 
tion. 



3. Kaon identification 

Kaons are distinguished from pions and protons on the 
basis of specific energy- loss measurements dE/dx in SVT 
and DCH and the number of Cherenkov photons and 
the Cherenkov angle in the DIRC. The difference be- 
tween the measured truncated- mean dE/dx in the DCH 
and the expected mean for the pion, kaon and proton 
hypothesis, with typical resolution of 7.5%, is used to 
compute likelihoods £,r, J~-k and Cp assuming Gaussian 
distributions. Similarly, the difference between the mea- 
sured 60% truncated-mean dE/dx in the SVT and the 
expected dE/dx is described by an asymmetric Gaussian 
distribution. For minimum-ionizing particles the resolu- 



tion on the SVT truncated mean is about 14%. In the 
DIRC, a likelihood is obtained for each particle hypoth- 
esis from the product of two components: the expected 
number of Cherenkov photons, with a Poisson distribu- 
tion, and the difference between the measured average 
Cherenkov angle to the expected angle for a given mass 
hypothesis, assuming a Gaussian distribution. 

For i?-flavor tagging the likelihood variables from SVT, 
DCH and DIRC are combined as inputs to a neural net- 
work whose output is a single discriminating variable for 
kaon selection. The network is trained with Monte Carlo 
simulation of generic B decays. The average efficiency 
of the selection is about 85% for a pion-misidentification 
probabilit y of a bout 2.5%. Further details are described 
in Section IV B. 



The exclusive reconstruction of many B meson final 
states does not generally require explicit kaon identi- 
fication. For some channels a VeryLoose kaon selec- 
tion based on likelihood ratios is imposed to reduce 
backgrounds to acceptable levels. The combined likeli- 
hood uses the individual likelihoods from SVT and DCH 
for momenta below 0.5GeV/c, from DCH only for mo- 
menta between 0.5 and 0.6GeV/c, and from DIRC only 
for momenta above 0.6GeV/c. Kaon candidates are re- 
jected if the likelihood ratios satisfy Ck/^tt < r and 
Cx/Cp < r, where r = 0.1 for p < 0.5GeV/c and r — 1 
for p > O.b GeV/c. Tracks with no particle information 
are assumed to be kaons. This VeryLoose kaon require- 
ment has a nearly constant kaon-identification efficiency 
of about 96% and a pion-misidentification probability 
of at most 15% for tracks in the transverse momentum 
range 1 to 2.5 GeV/c. Tighter kaon selections require 
^k/^-k > 1", with r typically greater than one. For a 
loose pion selection, candidates are rejected if they sat- 
isfy tighter kaon or lepton criteria. 
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III. RECONSTRUCTION OF B MESONS 



Neutral B mesons in flavor eigenstates are recon- 
structed in the hadronic final states i?° Z?^*^~7r+, 
D^*^-p+, £)(*)-a+, and J/^P K*°{K+tt-), and the 
semileptonic decay mode D*^i^v. The CP sample 

is reconstructed in the channels — > J/^K°, 'ip{2S)K°, 
XciK'^, J/i'K*° {K*° K°n°) and J/ipK°. In some 
cases, control samples of charged B decays are stud- 
ied, where the hadronic final states B~^ — + _D(*)''7r+, 
J/^K^*'>+, ^P{2S)K+ and XciK+ are used. All final- 
state particles, with the exception of the neutrino in the 
semileptonic decay, are reconstructed. A number of Z?" 
and D~ decay modes are used to achieve reasonable re- 
construction efficiency despite the typically small branch- 
ing fractions for any given B oi D decay channel. A sum- 
mary of the various reconstructed B samples and purities 
is provided in Table III. 



TABLE III: Event yields for the different samples used in this 
analysis, before any tagging or tagging vertex requirements. 
The yields, purity, and signal size for B decays to hadronic 
final states are obtai ned f rom a fit to the tties distribution 
described in Section [II C, after selection on /S.E. Purities 
are quoted for jties > 5.27MeV/c^. The results for J/i> K° 
are obt ained from a fit to the AE distribution described in 
Section [II ^ The purity for J/tp Kl is quoted for events with 
Ai5 < lOMeV. The results for D'^l'^v are obtained from a 
distribution described in Section 
T.l < coses-D't < 1.1. 



fit to the cos 6 B~D* 
Purity is quoted for 



[HE 



Sample 


Final state 


Signal 


Purity 








(%) 


Bcp 




461 ± 22 


99 






113 ± 12 


93 






86± 17 


96 






44 ±8 


98 






64± 10 


74 






257 ± 24 


60 




Total 


1025 ±41 


83 






2380 ± 57 


92 




D*-p+ 


1438 ± 52 


84 






1146 ±45 


80 




D--K+ 


2685 ± 65 


83 




D-p+ 


1421 ± 57 


74 




D'at 


845 ± 44 


67 






1013 ±36 


95 




Total 


10941 ± 133 


83 


B+ 




6850 ± 102 


83 




D*\+ 


1708 ± 51 


91 




J/^K+ 


1921 ±46 


97 




iP{2S)K+ 


292 ± 18 


98 




XciK+ 


195 ± 29 


95 






384 ± 25 


87 




Total 


11343 ± 129 


86 


Semi- 




29042 ± 1500 


78 


leptonic 5° 









A. Event selection 

Multihadron events are selected by demanding a min- 
imum of three reconstructed charged tracks in the polar 
angle range 0.41 < 6iab < 2.54 rad. Charged tracks must 
be reconstructed in the DCH and are required to origi- 
nate within 1.5 cm in xy and 10 cm in z of the nominal 
beamspot. A primary vertex is formed on an event-by- 
event basis from a vertex fit to all charged tracks in the 
fiducial volume. Tracks with a large contribution to 
the vertex fit are removed until an overall probability 
greater than 1% is obtained or only two tracks remain. 
The resolution achieved by this method is about 70 fim 
in X and y for hadronic events. Events are required to 
have a primary vertex within 0.5 cm of the average po- 
sition of the interaction point in the plane transverse 
to the beamline, and 6 cm longitudinally. Electromag- 
netic bumps in the calorimeter in the polar angle range 
0.410 < Blab < 2.409 rad that are not associated with 
charged tracks, have an energy greater than 30MeV, and 
a shower shape consistent with a photon interaction are 
taken as neutrals. A total energy greater than 4.5 GeV in 
the fiducial regions for charged tracks and neutrals is re- 
quired. To reduce continuum background, we require the 
normalized second Fox- Wolfram moment i?2 of the 
event, calculated with both charged tracks and neutrals, 
to be less than 0.5 (0.45) in hadronic (semileptonic) decay 
modes. The Fox- Wolfram moment is the momentum- 
weighted sum of Legendre polynomial of the order 
computed from the cosine of the angle between all pairs 
of tracks. The ratio i?2 provides good separation between 
jet-like continuum events and more spherical BB events. 



B. Reconstruction of decay daughters 

The reconstruction of B mesons typically involves the 
summation of a set of related decay modes, with multi- 
ple decay chains for the charm daughters or other short- 
lived decay products. To simplify analysis of such com- 
plex decay chains, virtual composite particles and their 
error matrices are constructed from the original daugh- 
ter particles. The composite particle then replaces the 
daughters in subsequent fits and analysis. The three- 
momentum of the virtual particle is fit directly, rather 
than computed from the updated daughters, improving 
speed and numerical accuracy. 

Vertex and kinematic fitting is used to improve four- 
momenta and position measurements, as well as to mea- 
sure the time difference between decaying B hadrons in 
the T{AS) BB decay. For example, in the case of 
J/ipK'^, the position measurement of the B^ can 
be improved with the constraint that the line-of-fiight 
of the if" intersects the J/tp vertex. The energy resolu- 
tion of the B'^ can also be improved by applying a mass 
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constraint to the J/V' and iiT^ daughters. Generahzed 
procedures have been developed and tested with con- 
straints implemented by the Lagrange-multiplier tech- 
nique. Possible constraints include a common decay ver- 
tex, mass, energy, momentum, beam energy (with and 
without smearing), beam-spot position and line-of- flight. 

Non-linearities in the fits require the use of an iterative 
procedure, where convergence is defined by demanding 
that the change in between two successive iterations 
is less than 0.01, within a maximum of six iterations. 
Simple fits involving only vertex constraints (except long- 
lived particles) are, however, accurate enough with a sin- 
gle iteration. 



1. 7r selection 

Neutral pion candidates are formed from pairs of EMC 
bumps with energy greater than 30 MeV, assumed to be 
photons originating from the interaction point. The in- 
variant mass of the photon pair is required to be within 
±20MeV/c2 (2.5a) of the nominal 7r° mass, with a min- 
imum summed energy of 200 MeV. Selected candidates 
are subjected to a kinematic fit with a tt*^ mass constraint. 
Within the acceptance of the EMC, efficiencies for this 
selection vary from about 55 to 65% for 7r° energies from 
0.3 to 2.5 GeV, typical of B decays. 



2. K*^ selection 

Candidates in the tt+tt^ mode are selected by 

requiring an invariant tt+tt" mass, computed at the ver- 
tex of the two tracks, between 462 and 534MeV/c^. The 
of the vertex fit must have a probability greater than 
0.1%. The angle between the flight direction and the 
momentum vector for the ifg candidate is required to be 
smaller than 200 mrad. Finally, the transverse flight dis- 
tance from the primary vertex in the event, r^y^ must be 
greater than 2 mm. 

Optimization for the reconstruction of the CP sample 
has produced slightly different selection criteria. The 
TT+TT^ invariant mass, determined at the vertex of the 
two tracks, is required to lie between 489 and 507MeV/c^ 
and the three-dimensional flight length with respect to 
the vertex of the charmonium candidate is required to be 
greater than 1 mm. 

Pairs of 7r° candidates, each in the mass range 
100-155 MeV/c^ (—5(7, +ia) and formed from non- 
overlapping EMC bumps, are combined to construct 
— > vr^vr" candidates. For each K^^ candidate with 
an energy greater than 800 MeV and a mass between 300 
and 700MeV/c^ at the interaction point, we determine 
the most probable decay point along the path de- 
fined by the initial momentum vector and the J/'ip 
vertex by maximizing the product of probabilities for the 
daughter tt*^ mass-constrained fits. Allowing for vertex 
resolution, we require the distance from the decay point 



to the J/V' vertex to be between —10 and -1-40 cm and the 
K^; mass, using the measured decay point, to be between 
470 and 536MeV/c2. 



3. K1 selection 

Candidates for mesons are identified in the EMC 
and IFR detectors as reconstructed clusters that can- 
not be associated with any charged track in the event. 
EMC candidates must have a cluster energy between 
200 MeV and 2 GeV and a polar angle 9 that satisfies 
cos9 < 0.935. To suppress backgrounds from tt" decay, 
candidates consistent with a photon are paired with 
other neutrals with > 30 MeV. Any candidate with 
100 < m(77) < 150 MeV/ is rejected. Likewise, clusters 
with more than 1 GeV energy that contain two bumps are 
rejected if the bump energies and shower shapes are con- 
sistent with two photons from a 7r° decay. Monte Carlo 
simulation shows that clusters due to true mesons are 
easily distinguished from candidates by these criteria. 
The remaining background consists primarily of photons 
and overlapping showers. Isolated clusters produced by 
charged hadrons are removed by the basic clustering al- 
gorithm, which requires a minimum separation of about 
20 cm between clusters. 

IFR candidates are defined as clusters with hits in two 
or more RPC layers that are not matched to any recon- 
structed charged track. To reduce beam-related back- 
grounds and to avoid regions where the charged tracking 
efficiency is low, we require that the polar angle 9 of 
the IFR cluster satisfy —0.75 < cos 6 < 0.93, and elimi- 
nate clusters that begin in the outer 25% of the forward 
IFR endcap. Due to the irregular structure of hadronic 
showers, some hits from charged tracks are missed by the 
tracking association. We suppress these clusters by re- 
jecting if ° candidates that lie within ±350 mrad in polar 
angle and in the range — 750(— 300) to -|-300(-|-750) mrad 
in azimuth of the EMC intersection of any positively 
(negatively) charged track in the event. The remaining 
background is predominantly from charged particles and 
detector noise. 

Some candidates satisfy both the EMC and IFR 
selection requirements. In the reconstruction of B'^ — > 

are 



J/ipKl, additional criteria described in Section HID 



applied to resolve the classification of the corresponding 
B candidates. Extensive studies of detection effi- 
ciencies have been conducted with a control sample of 
radiatively produced 4> mesons, decaying to K^K^. 



4- Selection of light resonances 



For p~ candidates, the tt^tt^ mass is required to lie 
within ±150MeV/c^ of the nominal p~ mass. The 7r° 
from the p~ decay is required to have an energy greater 
than 300 MeV. We reconstruct if*" candidates in the 
K~^TT~ and K'^tt'^ modes, while _fi'*+ candidates are re- 
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constructed in the K'^n'^ and i^^Tr"*" modes. The invari- 
ant mass of the two daughters is required to be within 
±100MeV/c^ of the nominal K* mass. Candidates in 



the mode 



TT+TT TT"*" are reconstructed by combin- 



ing three charged pions, with invariant mass in the range 
of 1.0 to 1.6GcV/c^. In addition, the probabihty of 
a vertex fit of the a'^ candidate is required to be greater 
than 0.1%. 



5. Charmed meson and charmonium selection 

The decay channels if+Tr", i4r+7r~7r° 
and ifg7r+7r~ are used to reconstruct D° candidates, 
while D~ candidates are selected in the K^tt^tt"' and 
ifgTT" modes. Charged and neutral kaons are required 
to have a momentum greater than 200 MeV/ c. The same 
criterion is applied to the pion in B° — > D(*)~7r+, 
D^*^~p'^ decay. For the decay modes B° D'-*^^ai, 
the pions are required to have momentum larger than 
150MeV/c. We require and D~ candidates to lie 
within ±3(7 of the nominal masses, where the error a 
is calculated event-by-event. The distributions of the 
difference between measured and nominal and 
meson masses, normalized by the measured error on the 
candidate masses, are found to have an RMS in the 
range 1.1-1.2 when fit with a Gaussian distribution. For 
D'^ _ftr+7r~7r°, we only reconstruct the dominant reso- 
nant mode D° K~^p', followed by p~ 7r~7r°. The 
angle ^'Ijo^r between the tt" and in the p~ rest frame 
must satisiy Icos^Jjo^l > 0.4. Finally, all and 
candidates are required to have a momentum greater 
than 1.3GeV/c in the T(4S') frame and a probabil- 
ity for the topological vertex fit greater than 0.1%. A 
mass-constrained fit is applied to candidates satisfying 
these requirements. 

We form D*^ candidates in the decay D*^ D^tt^ by 
combining a D'^ with a pion that has momentum greater 
than 70 MeV/ c. The soft pion is constrained to originate 
from the beamspot when the D*^ vertex is computed. 
To account for the small energy release in the decay 
T(4S') — > BB (resulting in a small transverse fiight of the 
B candidates), the effective vertical size of the beam spot 
is increased to 40 /Ltm. Monte Carlo simulation was used 
to verify that this does not introduce any significant bias 
in the selection or in the At measurement. After applying 
a mass-constrained fit to the D'^ daughter, D*~ candi- 
dates are required to have m(£'°7r~) within ±1.1 MeV/c^ 
of the nominal D*~ mass for the — > ii'+7r~7r° mode 
and ±0.8MeV/c^ for all other modes. This corresponds 
to about ±2.5 times the RMS width of the signal distri- 
bution, which is estimated by taking a weighted average 
of the core and broad Gaussian components of the ob- 
served m{D^'iT^) distributions. 

We form D*^ candidates by combining a with a 
7r° with momentum less than 450MeV/c in the T(4S') 
frame. D*^ candidates are required to have m{D^TT^) 
within ±4 MeV/c^ of the nominal value, after applying a 



mass-constrained fit to the daughter. 

Candidates for J/ip and ^^(25) mesons are recon- 
structed in their e^e~ and p.^ pT decay modes, while 
ijj{2S) mesons are also reconstructed in the J/ip7T^7T~ 
channel. Table |^ shows the particle identification and 
invariant mass requirements for the e^e~ and p^p^ 
daughters. These vary with reconstructed B decay chan- 
nel due to the differing levels of background encountered. 
For J/ip e+e~ and ip{2S) — > e+e~ decays, where the 
electron may have radiated Bremsstrahlung photons, the 
missing energy is recovered by identifying clusters with 
more than 30 MeV lying within 35 mrad in polar angle 
and 50 mrad in azimuth of the electron direction pro- 
jected onto the EMC. 

For the ip{2S) J/i/jtt+tt" mode, J/ip candidates 
are constrained to the nominal mass and then combined 
with pairs of oppositely-charged tracks considered as pi- 
ons, with invariant mass between 400 and 600 MeV/ c^. 
Candidates with 0.574 < m( J/^'tt+tt") — m{J/ip) < 
0.604 GeV/c^ are retained. 

Photon candidates used for the reconstruction of 
Xci ~^ J/ipl 3're required to lie within the calorimeter 
fiducial volume (0.41 < 9j < 2.41 rad) and have an en- 
ergy greater than 150 MeV. In addition, the candidate 
should not form, in combination with any other photon 
in the event having at least 70 MeV of energy, a tt^ candi- 
date with mass between 120 and 150 MeV/ c^. The invari- 
ant mass of the Xci candidates is required to be greater 
than 3.476 and smaller than 3.546 GeV/c^. 



C. B meson selection in fully-reconstructed modes 

We reconstruct B candidates in all modes except B° 
J/tpKl and B'^ D*~t'^v using a pair of nearly uncor- 
related kinematic variables, the difference A_B between 
the energy of the B candidate and the beam energy in 
the T{AS) center-of-mass frame, and the beam-energy 
substituted mass, mES, defined as 



P 



(15) 



where s is the square of the center-of-mass energy, Ei and 
Pi are the total energy and the three momentum of the 
initial state in the laboratory frame, and p is the three 
momentum of the B candidate in the same frame. For the 
purpose of determining event counts and purities, a signal 
region is defined in the (mEs, Ai5) plane as 5.27 < mss < 
5.29GeV/c2 and |A£:| < 3cr(A£'), where cr(A£') is the 
resolution on Ai?. Likewise, a sideband region is defined 
as 5.20 < mES < 5.26GeV/c2 and |A£:| < 3cr(A£;). The 
value of (T(Ai?) is mode-dependent and varies between 
7 to 40 MeV as measured in the data. When multiple 
B candidates (with toes > 5.20MeV/c^) are found in the 
same event, the candidate with the smallest value of | Ai5| 
is selected. 

Two types of background in the sample of selected 
B^ candidates are distinguished. The first background. 
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TABLE IV: Particle identification and invariant mass requirements for J/ip and ipi^S) f candidates. The minimal 

particle identification criteria are applied to both daughters, wh ile o nly one daughter must pass the restrictive requirement. 
Electron and muon selection requirements are defined in Section [ID Mass ranges are quoted in GeV/c'^ and MIP refers to a 
minimum-ionizing particle. 



B channel 


e+e candidates 


candidates 


Minimal 


Restrictive 


m(e^e ) 


Minimal 


Restrictive 




XciKl (J/^7) 
J/tPK* 


None 
VeryLoose 
VeryLoose 
Loose 
Tight 
Loose 


Tight or noCal 
Tight 
Tight 
Tight 
Tight 
VeryTight 


2.95-3.14 
3.436-3.736 
2.95-3.14 
2.95-3.14 
2.95-3.14 
3.00-3.13 


MIP 
VeryLoose 
VeryLoose 
VeryLoose 
Loose 
Loose 


Loose 
Loose 
Loose 
Loose 
Loose 
Tight 


3.06-3.14 
3.06-3.14 
3.06-3.14 
3.06-3.14 
3.06-3.14 
3.06-3.13 



called combinatorial, arises from random combinations of 
charged tracks and neutral showers from both B mesons 
in BB events or from continuum events. This back- 
ground is smoothly distributed in rngs and does not peak 
near the B mass. The second, so-called "peaking" back- 
ground, consists of events in which, for example, a slow 
pion from the reconstructed B meson is replaced by a 
slow pion from the tagging i3, causing an enhancement 
near the nominal B mass. The peaking background from 
charged B decays is considered as a specific background 
source in the construction of the full likelihood function 
for B^-B'^ mixing, since these events have a particular 
time structure and set of effective dilutions. In this case, 
the peaking background from other neutral B decays has 
time-dependent properties and dilutions that are essen- 
tially identical to the signal and is treated as such. For 
the likelihood describing the CP sample, the peaking 
background is simply assumed to have zero effective CP. 

Suppression of continuum background, in addition to 
a general requirement on i?2, is typically provided by 
restricting the thrust angle 0th, defined as the angle be- 
tween the thrust axis of the particles that form the re- 
constructed Brec candidate and the thrust axis of the 
remaining tracks and unmatched clusters in the event, 
computed in the T{AS) frame. The two thrust axes are al- 
most uncorrelated in BB events, because the mesons 
are nearly at rest in the T(AS) rest frame. In continuum 
events, which are more jet-like, the two thrust axes tend 
to have small opening angles. Thus, a requirement on 
the maximum value of | cos 0th | is effective in continuum 
rejection. 

Signal yields and sample purities are extracted from 
fits to the Toes distributions of B candidates with a 
Gaussian distribution for the signal and an ARGUS back- 
ground shape |23j for the combinatorial background with 
a functional form given by 



^(toes ; mo , = ^-B Toes V 1 



'-ES 



(16) 



for xes = toes/toq < 1, where toq represents the kine- 
matic upper limit and is held fixed at the center-of-mass 
beam energy = 5.291 GeV, and ^ and Ab are free 
parameters. 



We assign background and signal probabilities to each 
event included in the likelihood fit based on the measured 
value for toes- However, it is the toes sideband region, 
where the background probabilities are essentially 100%, 
that dominates the determination of the combinatorial 
background fraction and At structure for background 
events under the signal peak. Monte Carlo simu- 
lation shows a modest TOes dependence on the composi- 
tion of the combinatorial background over the sideband 
range toes > 5.2 GeV/c^ through the B signal region, 
due to variation of the fraction of continuum versus BB 
contributions. Since these two sources have different At 
behaviors, the changing composition leads to a small cor- 
rection and systematic error on the precision mixing mea- 
surement, but is negligible for the sin2/3 extraction. The 
fraction of peaking backgrounds from charged B decays 
are estimated with Monte Carlo simulation as described 
in the following sections. 



1. B decays to flavor-eigenstates 

Candidates in the -Bfiav sample of neutral flavor- 
eigenstate B mesons are formed by combining a, D*~ or 

( 



D with a 7r+, p+ (p^ 



-tt"), at {a+ 



or by combining a J/ip candida te w ith a K {K 
A'+TT^). As described in Section IIP , kaon identification 
is used to reject background. For most B'^ modes, it is 
possible to achieve signal purities of at least 90% with the 
VeryLoose selection, or no particle identification at all. 
However, for the mode B^ D~at, the tighter kaon 
identification is required to reduce large combinatorial 
backgrounds. 

For final states with a D* and 2 (3) pions we require 
Icos^thl < 0.9 (0.8) forjhe D" K+ir^ and K+tt'tt° 
modes and 0.8 (0.7) for D" -> K+n+n-n- and is:°7i-+7r-, 
while no requirement is made for the — s- D*~Tr~^ mode. 
In modes which contain a D~ and a 7r+, or in 
the final state, we require |cos0th| < 0.9, 0.8, or 0.7, 
respectively. 

The -B° signal yield and sample purity extracted from 
fits to the Toes distribution are summarized in Table IH . 
The net signal sample, before applying any decay ver- 
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FIG. 3: Distribution of mES for all selected flavor-eigenstate 
B° candidates in hadronic decays to (a) open charm and (b) 
charmonium final states. Overlaid in both cases is the result 
of a fit with a Gaussian distribution for the signal and an 
ARGUS function for the background. 



tex requirements, consists of 9922 ±129 signal candidates 
in open charm decays with a purity of about 82%, and 
1013 ± 36 in the decay 5° J/iPK*^ {K*" K+n- ), 
with a purity of about 95%. Figure |^ shows the combined 
mES distribution for all the hadronic B'^ modes. Super- 
imposed is the result of a fit with a Gaussian distribution 
for the signal and an ARGUS background form . 

The signal obtained by this method includes a small 
fraction of peaking background from other charged and 
neutral B decay modes. However, only the charged B 
component needs to be determined, since it alone has a 
time structure that differs from the signal events. There- 
fore, the fraction of peaking background is estimated with 
a sample of r{AS) B+B' Monte Carlo events. The 
_B+ mesons are forced to decay in the decay modes D*^ 
or with a tt"*", p"*", or a^, since the main source is 



decay channels that have one more or one fewer pion in 
the final state than the signal modes of interest. We then 
attempt to reconstruct neutral B mesons in the channels 
used for the -Bflav sample in data. A small peak at the B^ 
mass, obtained with the charged B Monte Carlo sample, 
leads to an estimate of (1.3 ± Q.2>'t^'^)% as the peaking 
component in the -Bflav signal. This result is obtained 
from a fit with a Gaussian distribution, whose mean and 



width are fixed by the signal parameters. The At 
structure of the peaking background in Monte Carlo is 
found to be consistent with the lifetime of the as 
expected. 



2. B'^ control samples 

The _B+ control sample of charged B candidates is 
formed by combining a D*'^, D", J/ip, or ■ip{2S) candi- 
date with a 7r+ or K^. For the £'°7r+ final state, we 
require |cos0th| < 0.9 for the K^tt^ mode and 

0.8 for all other channels. In modes that contain 
a D*'^ D°Tr°, the requirement is |cos0th| < 0.9 for 
D° K+TT-, 0.8 for the /f+7r^7r° and if+7r"7r"7r+, and 
0.7 for K^tt+t:-. 

The B^ signal yield and sample purity extracted from 



fits to the m^s distribution are summarized in Table II] 



The net signal sample in open charm modes, be 
fore applying any decay-vertex requirements, consists of 
2797 ± 62 signal candidates in charmonium modes, with 
a purity of about 94%, and 8547 ± 115 signal candidates 
in open charm modes, with a purity of about 84%. Fig- 
ure I shows the combined toes distribution for all the 
hadronic B^ modes. Superimposed on the data is the 
result of a fit with a Gaussian distribution for the signal 
and an ARGUS background form . 



3. -B" decays to CP modes involving K*^ 

We form the Bcp sample of neutral B candidates 
in charmonium modes with a K^^ by combining mass- 
constrained J/ip, il){2S) or Xci candidates with mass- 
constrained candidates, following the techniques of 
our recent branching- fraction study p^ . The helicity 
angle Oh of the J/ij] daughters with respect to the J/ip 
flight direction in the B candidate rest frame should 
have a sin^ 9h distribution. Therefore, we require that 
I cos9h\ < 0.8 for the e+e~ mode and 0.9 for the /i^/i^ 
mode, as an efficient way of rejecting backgrounds. For 
the 'ijj{2S)K^ candidates, |cos0/i| of the ^j{2S) must be 
smaller than 0.9 for both leptonic modes. 

Distributions of toes are shown in Fig. || for the CP 
samples. Signal event yields and purities, determined 
from a fit to the toes distributions after selection on AE, 



are summarized in Table III 



The fraction of peaking background has been esti- 
mated with a sample of i? ^ J/ipX Monte Carlo events. 
The main source is decay channels that have one more or 
one less pion in the final state than the signal mode. The 
fractions are obtained by fitting the misreconstructed 
B —^ J/ipX sample with a Gaussian distribution, whose 
mean and width are fixed by the B^ signal parame- 
ters. The estimated contributions are (0.41 ± 0.09)%, 
(1.2±0.2)%, (2.9±1.7)%, and (1.1±1.1)% for the Jf-ipK^ 
{K"^ TT+TT-), J/ipK° {K° -> Tr^TrO), V'(2S')ii'" and 
Xci^s channels respectively. 
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FIG. 4: Distribution of tties for all selected flavor-eigenstate 
candidates in hadronic decays to (a) open charm and (b) 
charmonium final states. Overlaid in both cases is the result 
of a fit with a Gaussian distribution for the signal and an 
ARGUS function for the background. 



In the case of the XciK^ mode we have also explored 
the possibility of contamination from Xc2K'^ events. 
These would have a very similar final-state signature, but 
opposite CP. However, this decay mode has never been 
observed and the rate is expected to be highly suppressed 
due to angular momentum considerations. Figure ^ 
shows the invariant mass difference, m{£ij) — m{£i), for 
the Xci daughters of the XaK^ and XciK^ 

candidates. The distribution is background subtracted 
with the TTiEs sideband and a fit with two Crystal Ball 
distributions 1 25 is superimposed, where the means have 
been fixed to the known Xci and Xc2 masses and the 
widths are forced to be equal. The fraction of Xc2K 
events in the selected sample is found to be consistent 
with zero and, from the fit, an upper limit of 3.5% at 
95% C.L. is set on the fraction of i? ^ Xc2K candidates 
in the selected sample. 
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4. B° decays to the CP mode J/ipK*° 

The Bcp sample is further enlarged by the addition 
of B° candidates in the mode J/il}K*° {K*^ K°n"). 
For this purpose, mass-constrained J/tp candidates are 



FIG. 5: Candidates for B° J/ipK^ where K°s decays to 
a) TT+TT" or b) TT^Tr"; Candidates for c) B° ^ ip{2S)K°s and 
d) _B° Xci^s — > tt+tt"). Overlaid in each case is the 
result of a fit with a Gaussian distribution for the signal and 
an ARGUS function for the background. 
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FIG. 6: Distribution of m{££y) for the Xci daughters of 
XciK^ and B° Xci^s candidates. The expected 
location of a Xc2 signal is indicated by the arrow. 



combined with K*'^ — > K^ir^ candidates to form a 
candidate. To reduce the combinatorial background, the 
angle between the flight direction of the and the vec- 
tor connecting the reconstructed vertices of the J/ip and 
the candidates is required to be less than 200mrad. 
Cross- feed background from other B J/^p X modes in- 
volving a tt" (which includes cross-feed from the CP mode 
itself) is suppressed by requiring the cosine of the helic- 
ity angle of the K*'^ in the 5° meson rest frame to be 
smaller than 0.95. Further details of the selection and 
analysis of this sample can be found in Rcf . [0 . 
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FIG. 7: Distribution of tties for selected J/ipK*^ combina- 
tions, where K*° K%it° . The arrows indicate the region 
between 5.273 and 5.288 GeV/c^ that is used to define the sam- 
ple of B'^ candidates. Monte Carlo estimates of the various 
background contributions are also indicated. 



The mss distribution for J/il)K*° (A'*" K^tt^) com- 
binations in data is shown in Fig. 0. Given the rela- 
tively tight criteria applied in the lepton identification of 
the daughters of the J/-0 candidates (see Tabic IV), the 
background is dominated by true J/ip mesons from B 
decays. Its composition can therefore only be estimated 



with Monte Carlo simulation and the rriEs distribution 
is not expected to follow the phase-space form of Eq. [ij. 
Monte Carlo simulation of events with true J/ijj candi- 
dates has been adjusted to match recent results for char- 
monium branching fractions in B decays and takes into 
account the indication of S'-wave B^ — > J/ijjK^n^ decays 
and contributions due to higher K* resonances reported 
in Ref. |l|. 

As a result, backgrounds are not estimated with a fit 
to the observed tties distribution (Fig. |7|), but rather 
by extrapolation of Monte Carlo background distribu- 
tions, normalized to the number of produced B mesons 
in the data. All J/ipK'^TT^ combinations in the range 
5.273 < Toes < 5.288 GeV/c^ are considered as candi- 
dates for this purpose. Estimates of the signal and back- 
ground contributions in the candidate sample, and the 
corresponding effective CP, after acceptance correction 
for the signal selection, is provided in Table ^ while the 
signal yields and purities in data are listed in Table III. 
The dominant source of cross-feed background, with zero 
effective CP, is B+ J/ijK*+ {K*+ K'^Tr+), where 
the daughter 7r+ is exchanged for a background tt". 



TABLE V: Signal and background estimates for the selected 
B° J/ipK*" {K*° K^Tv°) sample. All the events within 
the range 5.273 < m-ES < 5.288 GeV/c^ are considered as B^ 
candidates and the background contributions are estimated 
with Monte Carlo simulation. The quoted errors are derived 
from conservative bound on the branching fractions and rep- 
resent the size of the variation used to estimate the systematic 
error on sin2/3 due to backgrounds. 



Event type 


Fraction (%) 


Effective CP 


Signal 


73.6 ± 7.4 


±0.65 ±0.07 


B+ ^ J/ijK'+(Kl'K+) 


17.4 ± 1.7 







2.4 ±0.7 


-1 


Higher K* resonances 


2.6 ± 1.3 


0± 1 


Non-resonant B° J/ifjKlir^ 


1.8 ±0.9 


0± 1 


Other B° J/^jjX 


2.4 ± 1.2 


0± 1 


Non B ^ J/xjjX 









D. B° decays to the CP mode J/iiKl 

Candidates for the Bcp sample in the mode B^ 
J/il^ Kl are obtained by combining mass-constrained J/^p 
and candidates, following the methods in Ref. [ p4| . 
The J/4' candidates are required to have a momentum 
in the T{4S) frame between 1.4 and 2.0GeV/c. As the 
K1 energy is not well measured by the EMC or IFR de- 
tectors, the laboratory momentum of the is deter- 
mined by its flight direction as measured from the EMC 
or IFR cluster and the constraint that the invariant mass 
of the J/tpK^ system has the known B° mass. The pro- 
duction angle 9b oi a. B meson with respect to the z 
axis in the T(4S') frame follows a sin^ 9b distribution. 
We require that |cos6'b| < 0.9. The J/V' helicity an- 
gle is required to satisfy |cos0;i| < 0.9 and the sum of 
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TABLE VI: Monte Carlo prediction for the composition of 
background channels containing a true J/tj) that pass the 
B" — > J/'tpK'i selection criteria. Events are required to have 
|AiJ| < lOMeV. The quoted errors are derived from con- 
servative bound on the branching fractions and represent the 
size of the variation used to estimate the systematic error on 
sin2/3 due to backgrounds. 



Event type 


EMC (%) 


IFR (%) 


Effective CP 




23 ±3 


26 ±3 


-0.68 ±0.07 


B+ -* J/i>K*^{K°Tv^) 


28 ±4 


45 ±6 







13 ±2 


2± 1 


-1 




3± 1 


5± 1 


+ 1 


B J/ipKl-K 


It? 







Other B° J/i^X 


32± 16 


21 ± 10 


0±0.25 



I cos^bI and | cos6'h| must be less than 1.3. Events with 
a reconstructed charged or neutral B decay to J/ipK^ 
(ifO ^ TT+TT- or ttOttO), J/?/>if*o [K*° K+n- or 
K^-K^), J/^pK+, or J/^PK*+ {K*+ K'^Tr+ or K+n°) 
are explicitly removed. The total missing transverse mo- 
mentum projected along the direction, where the to- 
tal momentum is calculated with all charged tracks and 
neutral clusters (without the A'°), must be no more than 
0.25 (0.40) GeV/c lower than the calculated transverse 
momentum for EMC (IFR) ifj^ candidates. 

Events where multiple J/ipK'^ combinations with 
AE < 80MeV satisfy these requirements are treated as 
a special case. A hierarchy is imposed where the highest 
energy EMC cluster for multiple EMC combinations, or 
the IFR cluster with the largest number of layers for mul- 
tiple IFR combinations is selected. In cases where there 
are both an EMC and IFR combination, the EMC com- 
bination is selected because it is expected to have better 
angular resolution. Although the EMC information is 
used, such events are counted as IFR events, since they 
have the same relatively high signal purity. 



TABLE VII: Results of the binned Ukelihood fit to the full AE 
distribution of the B" — > J/ip K1 combinations. All signal 
yields and background estimates are reported for the Bcp 
candidate range |Ai5| < lOMeV. 





Kl reconstruction type 




EMC & IFR 


EMC 


IFR 


Data events 


427 


228 


199 


Signal 


257 ± 24 


128 ± 17 


129 ± 17 


J/^X bkgd 


154 ± 15 


89 ± 11 


65 ± 10 


non- J/ip bkgd 


19 ± 2 


14 ± 2 


5 ± 1 


Signal fraction 


0.60 ± 0.04 


0.56 ± 0.05 


0.65 ± 0.05 



The difference AE between the energy of the J/tp 
system and the beam energy in the T{AS) frame is used 
to discriminate between signal and backgrounds. The 
AE distribution of selected B'^ JhpKl combinations 
for the T{AS) data is shown in Fig. ||. Signal events are 
peaked within ztlOMeV of AE = while background 



events extend towards positive values of AE. The small 
signal width and the asymmetric distribution of the back- 
ground in comparison with the modes are both con- 
sequences of the mass constraint used to determine the 
momentum. 

The purity of the B^ — > J/ipK'^ sample is the low- 
est of the CP modes (60%). Irreducible backgrounds are 
dominantly from B J/ipKlX modes, which cannot 
be distinguished from signal due to imposition of the ttib 
mass constraint in determining the momentum of the K'l 
candidate. The largest single background contribution is 
from B — > J/ipK*, where the K* decays to K^^tt. This 
mode and backgrounds from other B — > J/ip X decays are 
studied with Monte Carlo simulation. The composition 
of the events that are included in the J/-0 sample is 
given in Tabic 0. The effective CP, after acceptance cor- 
rection for the signal selection, is also provided. The ad- 
ditional background from events with a misreconstructed 
J/'ip U candidate is studied with the m{i£) sidebands. 

A binned likelihood fit to the AE distribution is per- 
formed separately for the EMC and IFR categories to 
determine the composition of the — > J/ipK^ sample. 
There are three fit components: the fraction of the data 
that is signal, the number of S ^ J/ipX background 
events, and the number of non-J/ip background events. 
The AE shapes for the signal and the J/ipX background 
are determined from Monte Carlo simulation. The AE 
shape of the non- J/t/i component is taken from the m,{U) 
sideband in the data. A Poisson term, with mean given 
by the expected number of non- J/i/; events in the m{£i) 
signal region, is included in the likelihood to constrain 
the normalization of the non-J/ip component. The result 
of the fit is shown in Fig. and the correspondin g sig nal 
and background fractions are reported in Table VII for 
the i?° J/4>K1 combinations that are selected as Bcp 
candidates in the interval |Ai5| < lOMeV. As expected 
from Monte Carlo studies, the purity of the IFR sam- 
ple is significantly better than the EMC sample, mainly 
because the B° J/tpK° {K° tt^tt^) background is 
significantly larger in the EMC sample. Since the puri- 
ties of the two subsamples are quite different, we obtain 
better statistical sensitivity in the sin2/3 fit by treating 
the EMC and IFR categories separately. 



E. Semileptonic decays 



The semileptonic decay B^ ^ D* 



wi 



th a 



sured branching fraction of (4.60 ± 0.27)% [|llj, is poten- 
tially a copious source of reconstructed B^ mesons. How- 
ever, since the neutrino cannot be detected, the back- 
ground levels in selected samples are generally larger and 
more difficult to characterize. Likewise, the Az deter- 
mination cannot take advantage of the beam spot and 
reconstructed i3° direction. As a consequence, we use a 
large sample of D*~i^V candidates only as a cross check 
on our determination of the mistag rates from the time 
structure of the i?flav and Bcp events. The selection 
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D* through its decay to -D°7r , where the three decay 
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FIG. 8: Distribution of AE for selected _B° J/^KI com- 
binations where the is identified a) in the EMC, b) in 
the IFR, or c) either subsample combined. The points with 
error bars are the data. The open histogram is the result of a 
three-component binned likelihood fit, where the three com- 
ponents are signal (open), inclusive J/i/' background (cross 
hatched), and non-J/tp combinatorial background (dark shad- 
ing). The shapes of the signal and inclusive J/ip background 
are taken from Monte Carlo. The shape of the non-J/^p com- 
binatorial background is taken from the m{££) sideband in 
data. Candidates for the Bcp sample are selected in the re- 
gion \AE\ < 10 MeV. 



criteria for this control sample and the characterization 
of backgrounds is described here, while the anal ysis of 

The 



the niistag rates is reported in Section VIII C 4 



modes K^tt , K'^tt^tt 
construct the D'^. 



TT and K^TT 7r° are used to re- 



1. Event selection 

All reconstructed D'-' candidates are required to have 
an invariant mass within ±2.5it of the nominal D'^ mass, 
based on the observed RMS width of the signal. A vertex 
fit to the Z?" candidate is required to have a probabil- 
ity greater than 0.001. There are no additional require- 
ments for D° K~^'K~. For — > ii'+7r+7r~7r~ and 
we require a VeryLoose kaon and a Loose 



pion particle identification as described in Section [ID, 



and a minimum tt*^ momentum of 200 MeV/ c in the labo- 
ratory frame. In addition, the K and tt candidates are re- 
quired to have momenta greater than 200 and 150 MeV/c, 
respectively, for the mode _fir+7r+7r~7r~. The decay 

D'^ — !• if+TT^Tr'^ occurs mostly through quasi-two-body 
channels. The p and K* resonances dominate and we 
use weights calculated from the Dalitz-plot position for 
each candidate |2^ to construct a probability per and 
select events using this quantity as a way of suppressing 
combinatorial background. 

D'^ candidates satisfying these requirements are com- 
bined with all charged tracks with a minimum transverse 
momentum of 50MeV/c and charge opposite to that of 
the kaon from the 13° to form D*~ candidates. The mass 
difference m{D^i:^) — m{D^) is required to lie within 
±2.5ct of the nominal value, based on the observed RMS 
width of the signal. 

Finally, D*~ candidates are combined with elec- 
tron or muon candidates satisfying the Tight lepton- 
identification requirements described in Section [ID. 



The lepton is required to have charge opposite to that 
of the D* and momentum greater than 1.2 GeV/c in the 
T{AS) frame. A vertex fit to the D*~£'^ candidate is 
required to converge and have a probability greater 
than 0.01. The D*~ and lepton from a true decay 
tend to be back-to-back in the rest frame, so we re- 
quire cosOo'-e < where 9d*~i is the angle between the 
D*^ and the lepton in the T{AS) frame. The cosine of the 
angle between the thrust axes of the (£>*"£+) pair and 
the rest of the event is required to satisfy | cos 0th| < 0.85, 
in order to reduce background from e^e~ —f cc events. 

The neutrino cannot be reconstructed in the detec- 
tor, but we can determine whether the missing four- 
momentum of the candidate is consistent with a particle 
of zero mass: 



[PB - PD* - Pef =pI = 0. 



(17) 



Applying this relation in the T{AS) frame, we obtain a 
constraint on the angle between the and the D*^l^ 
system: 



semileptonic sample is obtained by reconstructing the 



cos Ob-D'I = 



2|Pi3||PL>* 



(18) 
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The energy Eb and magnitude of the momentum jp^j of 
the initial-state i?° are known in the T(45) frame from 
the boosted beam energies. The energy E]j*(, the mag- 
nitude of the momentum |pd*«|, and the invariant mass 
TO£i«£ of the Z3*^£+ system are obtained from the four- 
momenta of the D*^ and the Icpton. The cosine of the 
angle cos 6b-D'1 should lie in the physical region [—1, -1-1] 
for true D*~£^V events. Allowing for detector resolution 
effects in the reconstructed momenta and angles, we re- 
quire \cos9B-D*e\ < l-l- 

After applying these selection criteria, we obtain a 
sample of about 37500 B° D*-£+V candidates. The 
cos9B~D'e distribution of these candidates is shown in 
Fig. I 




B-D*l 



FIG. 9: Extraction of BB and D** backgrounds from the 
data. The points show the cos9b-d'i distribution for the 
reconstructed signal. The histogram shows the result of the fit 
described in the text. The different background contributions 
are indicated by the various hatchings. 



2. Sample composition 

The final sample contains a fraction fsig of 
D*~l^v signal events, as well as fractions /comb of D*~ 
combinatorial background, /og of continuum background 
(a true D*~ and an identified lepton), fB+-tD*xeu of 
B 0*71 £V events, where the D*it can come from a ra- 
dially or orbitally-excited L = 1 state or non-resonant de- 
cay, and finally /^^g- of background from BB events with 
a true D*~ and an identified lepton {BB D*Y£). Ex- 
amples in this last category are cases where the D*~ and 
the lepton come from two different B mesons or where 
the lepton and the D*~ are from the same B but the 
lepton is produced in a charm decay. 

The combinatorial and continuum backgrounds can be 
extracted directly from the data. The former is deter- 
mined with the ra{D^TT^) — ra{D^) distribution itself. 
The latter is estimated with the off-resonance data sam- 
ple, weighted by the ratio of the relative integrated lu- 
minosities for on- and off-resonance data. The remain- 
ing three contributions can be distinguished by their dif- 



ferent distributions in cosOb-D'I- The B^ D*~£'^v 
signal events are expected to lie in the region —1 < 
cosOb-D'I < 1, while contributions from B — > D*tt£17 
semilcptonic decay, due to missing particles, must extend 
below the kinematic threshold cos6b-D'1 < Finally, 
BB background events populate the hdl cos6b-d*i dis- 
tribution. Thus, the region cos6b-D'1 > 1 contains 
mainly BB background, while the region cos9b-d*i < 
— 1 is mostly populated by i? D*it£I'. The shape of 
the cosOb-d'I distributions for these three components 
is obtained from Monte Carlo simulation and a fit to the 
full cosOb-d*^ range is used to determine the two back- 
ground fractions in the signal region. 

The orbitally-excited resonances that can decay to 
D*n are the two narrow states Di, D2 (observed with 
masses around 2420 and 2460MeV/c2) and the broad 
state D'l (not yet seen, but with mass expected to be 
about 2420MeV/c2 and F > 250 MeV). Contributions 
from B D* X£T' decays with more than one pion are ex- 
pected to be small and are more easily separated from the 
signal with cos9b-d*i- Isospin symmetry requires that 
the charged B contribution be 2/3 of the total D*~7r 
pairs from B D*tt£U decays, either from orbitally- 
excited states or non-resonant decays. The cos9b-d*i 
distribution obtained from Monte Carlo simulation for 
the different channels is modeled with a general func- 
tion that is sufhciently flexible to describe both individ- 
ual channels as well as a superposition of excited charm 
modes. 

After subtraction of continuum and combinatorial 
backgrounds, a fit is performed to the resulting 
cos9b-d*i distribution over the full observed range 
[—8, -|-5]. However, it is only the relative fraction of the 
various backgrounds in the signal window cos9B-D'e. € 
[—1.1,1.1] that we require. Furthermore, in the case of 
£)**, only the charged B contribution is a background for 
the measurement of the mistag fraction and is assumed to 
be 2/3 of the total D** contribution. The fitted fractions 
are defined by: 



9* 



9BB 



N{B+ D*X£V)/[N{B° D*-£+V) + 

N{B^ D*X£+V) + N{B+ D*X£j7)] 

N{BB D*Y£)/[N{B° ^ D*-£+V) + 

N{B" D*X£+V) + iV(BS D*Y£)] (19) 



where the N is the number of events from a given process 
that survives the selection requirements. 

The result of the fit to the full untagged sample is 
shown in Fig. ^, along with the Monte Carlo model 
for the D** component. The of the fit in the full 
cos9B-D'e range is 82 for 69 degrees of freedom. The 
fitted contributions are g** = (4.5 ± 0.3 ± 2.2)% and 
QBB = (4.8 ± 0.4 ± 2.2)%, where the first error is sta- 
tistical and the second systematic. To estimate the sys- 
tematic error on these fractions, three extreme assump- 
tions have been made concerning the B D*X£V back- 
ground: all narrow D** states, all broad D** states, or all 
non-resonant decays. The largest deviation comes from 
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the non-resonant model. Another source of systematic 
uncertainty is the assumed form for the cos9b-d*i dis- 
tributions in the Monte Carlo simulation. The contribu- 
tion from this effect has been estimated by incorporating 
a 30% fraction with a uniform distribution. 

The absolute backgrou nd fr actions in the untagged 
sample are given in Table VIII, where the uncertainties 
include both statistical and systematic contributions. 

The sum of the fractions of signal and background con- 
tributions is constrained to unity. On this basis the signal 
component is found to be /sig = (78 ±4)% leading to an 
estimated yield of 29042 ± 1500 5° D*-l+V signal 
events. 



TABLE VIII: Sample composition in data as determined from 
fits to the cos9B-D*e distributions. The fractions have been 
computed without a requiring tagging information. The dom- 
inant errors are systematic except for fee, which is limited by 
the statistics. 



comb 



fo 



ff 



0.139 ± 0.028 0.008 ± 0.002 0.039 ± 0.018 0.037 ± 0.018 



IV. FLAVOR TAGGING 

After the daughter tracks of the reconstructed B are 
removed, the remaining tracks are analyzed to determine 
the flavor of the Btag , and this ensemble is assigned a tag 
flavor, either B'^ or B'-\ 

We use four different types of flavor tag, or tagging 
categories, in this analysis. The first two tagging cate- 
gories rely upon the presence of a prompt lepton, or one 
or more charged kaons in the event. The next two cat- 
egories exploit a variety of inputs with a neural-network 
algorithm. These tagging categories are hierarchical and 
mutually exclusive. 

To quantify the discriminating power of each tagging 
category, we use as a figure of merit the effective tagging 
efficiency Qi = (1 — 2wif' , where is the fraction of 
events associated to the tagging category i and Wi is the 
mistag fraction, the probability of incorrectly assigning 
the tag to an event in this category. The statistical er- 
rors in the measurements of sin2/3 and Amd are inversely 
proportional to \J^2~Qi- 

The mistag fractions are measured with the -Bflav data 
sample. The results are shown in Section VII. The per- 



formance of the tagging algorithm, \/^~Qi, was opti- 
mized and the neural networks were trained with Monte 
Carlo simulations only. Differences between the tagging 
inputs in data and in simulation may make the actual 
tagging algorithm somewhat non-optimal, but would not 
lead to a bias because the wrong-tag fractions Wi are 
measured directly with data, both for the mixing and 
CP-violation measurements. 



A. Lepton and kaon tags 

The Lepton and Kaon tagging categories use the cor- 
relation between the flavor of the decaying h quark and 
the charge of a primary lepton from a semileptonic decay 
or the charge of a kaon from the chain 6 ^ c — > s. 

For the Lepton category we use both electrons and 
muons, which are required to pass the VeryTight and 
Tight selection, respectively (see Tables | and ||). A 
minimum requirement of 1.0 (l.l)GeV/c on the elec- 
tron (muon) center-of-mass momentum is applied to re- 
duce the contamination from softer, opposite-sign lep- 
tons coming from cascade semileptonic decays of charm 
mesons. The center-of-mass momentum spectra for elec- 
trons and muons are compared to simulation in Fig. |l^ 
for the Sflav sample, after background subtraction based 
on the mEs sideband events. In each event, the electron 
or muon with the greatest center-of-mass momentum is 
used for flavor tagging; for the tiny fraction of events with 
both an electron and muon, the electron is used due to 
its smaller misidentification rate. 

The kaon content of the event is evaluated by taking 
the sum of the charges of all kaons identified with a neu- 
ral n etwor k algorithm (K subnet described below in Sec- 
tion IV B ). The kaon identification algorithm has been 



set to maximize the effective tagging efficiency Q. There 
are 0.8 charged kaons per B decay, and roughly 15% of 
these have the wrong sign (e.g. K~ from B^ , rather than 
the expected K'^). Wrong-sign kaons occur primarily in 
B decays to a charmed-anti-charmed pair of mesons. The 
momentum distributions are quite similar for right- and 
wrong-sign kaons; we find no kinematic quantity that 
distinguishes between them. The center-of-mass momen- 
tum spectrum for charged kaons and the distribution of 
charged kaon multiplicity are shown in Fig. 11 for the 
Bflav sample. 

An event with an identified high-momentum lepton is 
assigned to the Lepton category unless the sum of the 
charges of any kaons present has the opposite sign. Next, 
events are assigned to the Kaon category if the sum of 
the kaon charges is non-zero. The charge of the lepton 
or sum of kaon charges is used to assign the fiavor of 
the -Stag- AH remaining events, approximately 55% of 
the total including those with inconsistent lepton and 
kaon charge (0.5% of all events in simulation) and those 
with two oppositely-charged kaons (4.6% of all events 
in simulation), are passed to the neural-network-based 
categories. 



B. Neural- network tags 

Besides identified high-momentum leptons and charged 
kaons, there are other features that can be used to de- 
termine the flavor of the Btag, although not as easily or 
cleanly distinguishing. These include soft pions from D* 
decays, high-momentum primary leptons that are not se- 
lected by the electron or muon identification algorithms, 
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FIG. 10: Center-of-mass momentum distribution for a) elec- 
trons and b) muons. Data from the Bflav sample, after back- 
ground subtraction based on the niES sideband, are shown 
as points. The open histogram shows primary leptons, the 
cross-hatched histogram cascade leptons, and the diagonally- 
hatched histogram fake leptons, all from simulation. The sim- 
ulation is normalized, with a residual overall systematic error 
of 5%, to the total number of 5° decays in data after back- 
ground subtraction, not to the number of observed leptons. 
The vertical lines at 1.0 GeV/c for electrons and 1.1 GeV/c for 
muons indicate the requirement on center-of-mass momentum 
for the Lepton category. 



lower-momentum primary leptons, and charged kaons 
that are not selected by the kaon identification algorithm. 
These sources are combined with a multivariate method; 
we use a sequence of neural networks to flavor-tag those 
events not assigned to the Lepton or Kaon categories. 

Three different track-based neural networks, called 
"subnets", are trained, each with a specific goal. The 
L, K, and Sof tPi subnets are sensitive to the presence of 
primary leptons, charged kaons and soft pions from D* 
decays, respectively. Each subnet is applied to all tracks 
from the i?tag- 

The L subnet uses the binary output of the electron 
and muon identification algorithms on the input track, 
the center-of-mass momentum of the input track, and a 
pair of kinematic variables, and cos 9i,j, that separate 
primary leptons from cascade leptons and other tracks. 

The isolation variable, E^^, is given by the sum of the 
energies of all tracks within 90° of the W direction. The 



FIG. 11: a) Center-of-mass momentum distribution for kaons 
and b) kaon multiplicity per event. Data from the -Bfiav sam- 
ple, after background subtraction based on the rriES sideband, 
are shown as points. The histograms are from simulation. In 
a), the diagonally-hatched histogram is from fake kaons, the 
cross-hatched histogram is from kaons that have the wrong- 
sign charge, and the open histogram is from kaons with the 
right-sign charge, all from simulation. The simulation is nor- 
malized to the total number of B^' flavor candidates after 
background subtraction, not to the number of observed kaons. 



W momentum is inferred as the sum of the input track 
momentum and the neutrino momentum, which we take 
to be the missing momentum in the center-of-mass frame 
using all charged tracks in the Stag- This variable is ef- 
fective because in a semileptonic decay the hadrons re- 
coiling against the virtual W would generally go off in 
the opposite direction. 

The other kinematic variable used, cos 9i^, is the co- 
sine of the angle between the input track and the neutrino 
direction. The distributions in the Bflav sample and sim- 
ulation of Ey^ and cos Bi^^ are shown in Fig. |lja and b, for 
all events not in the Lepton or Kaon category. The cor- 
responding distribution of the L subnet output is shown 
in Fig. [l3|a. 

The K subnet uses the input track momentum in the 
laboratory frame, together with the three relative like- 
lihoods Ck/{C„ + Ck) for the SVT, the DCH and the 
DIRC. The SVT and DCH likehhoods are derived from 
dE/dx measurements, and the DIRC likelihood is calcu- 
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FIG. 12: Inputs to the subnets: a) E^, b) cos9ii,, c) cos^th 
for low center-of-mass momentum tracks {p* < 0.18GeV/c), 
and d) the center-of-mass momentum for all tracks. The 
points are data from the Bflav sample after background sub- 
traction based on the thes sideband, and the histogram is sim- 
ulation. For cos^th the diagonally-hatched histogram shows 
the contribution from soft-vr coming from D* decays, and for 
the other distributions shows the component from primary 
leptons. The simulation is normalized to the total number of 
_B° flavor candidates after the background subtraction. 



lated from a global fit to the number of photons detected, 
their positions and arrival times relative to the corre- 
sponding track. The distribution of the K subnet output, 
again for events not in the Lepton or Kaon category, is 
shown in Fig. [l^. 

The SoftPi subnet uses the center-of-mass momen- 
tum of the input track, the cosine of the angle of the 
input track with the thrust axis cos 0th, and the center- 
of-mass momentum of the minimum momentum track. 
The thrust axis is determined from all charged tracks and 
neutral clusters in the -Stag- The direction of any D* in 
the decay of the Stag is approximated by the direction of 
the thrust axis. Thus soft pions from D* decays, which 
are aligned with the D* direction in the center-of-mass 
frame, tend to be correlated with the thrust axis. The 
distribution of cos 9th is shown for low center-of-mass mo- 
mentum tracks in Fig. ^2|c, comparing the i?fiav sample 
with simulation for all events not in the Lepton or Kaon 
category. The corresponding distribution of the SoftPi 
subnet output is shown in Fig. p^ . 

The outputs of the three subnets are among the inputs 
to a final neural network, which is trained to distinguish 
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FIG. 13: Output of the subnets for events not assigned to 
the Lepton or Kaon categories: a) L subnet, b) K subnet, and 
c) SoftPi subnet. The points are data from the Bflav sam- 
ple after background subtraction based on the tues sideband, 
and the histogram is simulation. For each distribution, the 
filled portion of the histogram shows the component with a 
(singly-hatched) or (cross-hatched) tag from the full neu- 
ral network algorithm that arises from true primary leptons, 
kaons, or soft-pions respectively. Note that the latter has the 
opposite charge correlation with the B'^ tag. The simulation 
is normalized to the total number of -B" flavor candidates after 
background subtraction. 



between i?" and B'^. The variables used as inputs to the 
final network include the maximal values of the L and 
SoftPi subnet outputs, each multiplied by the charge of 
the corresponding input track, and the highest and the 
second-highest values of the K subnet output again mul- 
tiplied by the charge of the corresponding input tracks. 
The two other inputs to the final neural network are the 
center-of-mass momentum of the maximum momentum 
track multiplied by its charge, and the number of tracks 
with significant impact parameters. The latter is an in- 
dicator of the presence of K'^ mesons. The distribution 
of the center-of-mass momentum for all tracks is shown 
in Fig. |l|d. 

The output from the final neural network, xnt, is 
mapped onto the interval [—1,1]. The assigned flavor tag 
is if xnt is negative, and otherwise. Events with 
\xnt\ > 0.5 are assigned to the NTl tagging category and 
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events with 0.2 < I^attI < 0.5 to the NT2 tagging cate- 
gory. Events with \xnt\ < 0.2, approximately 30% of 
the total, have very little tagging power and are rejected. 
The distribution of x nt for all events not assigned to the 
Lepton or Kaon category is shown in Fig. [l4|a . 

Most of the separation between and in the NTl 
and NT2 tagging categories derives from primary Icptons 
that are not identified as electrons or muons and from soft 
pions from D* decays. Simulation studies indicate that 
roughly 37% of the effective tagging efficiency Q is due to 
events with unidentified primary leptons, 28% is due to 
events with a soft pion, a further 11% from events with 
a lower momentum primary lepton, and the remainder 
from a mixture of the various inputs. This classification 
is shown for the NTl and NT2 categories in Fig. [l|b for a 
simulation of 5° decays. The modest disagreements be- 
tween Monte Carlo simulation and data that are evident 
in the distributions of the subnet output variables shown 
in Fig. |l^ lead to a difference in the predicted value of 
Q = (3.0 ±0.1)% and (1.4 ±0.1)% for NTl and NT2 cate- 
gories in simulation versus (2.5 ±0.4)% and (1.2 ±0.3)% 
as measured in data. 



V. TIME DIFFERENCE MEASUREMENT 

The difference between B decay times, At = tree ~^tagj 
is determined from the measured separation Az between 
the vertex of the reconstructed B meson (-Brcc) and the 
vertex of the flavor-tagging B meson (Stag) along the z 
axis. The Az resolution is dominated by the z position 
resolution for the -Btag vertex. 
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A. Az reconstruction 

In the reconstruction of the B^^c vertex, we use all 
charged daughter tracks. Daughter tracks from K'^ and 
D candidates are first fit to a separate vertex and the 
resulting parent momentum and position are used in the 
fit to the iJj-cc vertex. Mass constraints, which include 
neutral daughters, are used for D candidates but not for 
D*-, J/ip and i/'(25') candidates. The RMS resolution 
in z for the Brec vertex in Monte Carlo simulation is 
about 65 /Ltm for more than 99% of the B candidates, and 
40 fim f or ab out 80% of the candidates. As described in 
Section VE. the resolution is about 5% worse in data 



than in Monte Carlo simulation. 

The vertex for the -Btag decay is constructed from all 
tracks in the event except the daughters of i?roc- For 
fully reconstructed modes, an additional constraint is 
provided by the calculated i?tag production point and 
three-momentum, with its associated error matrix. This 
is determined from the knowledge of the three momen- 
tum of the fully reconstructed B^ec candidate, its decay 
vertex and error matrix, and from the knowledge of the 
average position of the interaction point and the T(45') 
average boost (see Fig. |l^). These Stag parameters are 



FIG. 14: a) Output of the final neural network for -Bflav events 
that are not assigned to the Lepton or Kaon category, where 
the points are from the -Bflav data after a background subtrac- 
tion and the histogram is simulation; b) Neural network out- 
put from simulation of single B'^ decays with no time evolu- 
tion, again for events not in the Lepton or Kaon category. The 
breakdown from bottom to top is events with two kaons or a 
kaon and lepton (KK-fKL), events with a soft pion (soft tt), 
events with a high momentum unidentified lepton (HPUL), 
events with a lower momentum lepton (LPL) , and all remain- 
ing events. The outermost bins correspond to the category 
NTl and the next to NT2. Entries for xmt > 0.0 represent 
correct tags, while those for xmt < are mistags in each of 
the categories. The center bin contains events for which no 
tag is assigned. 



used as input to a geometrical fit to a single vertex, in- 
cluding all other tracks in the event except those used to 
reconstruct -Broc- In order to reduce bias and tails due 
to long-lived particles, and candidates are used as 
input to the fit in place of their daughters. In addition, 
tracks consistent with photon conversions (7 e+e^) 
are excluded. To reduce contributions from charm decay 
products, which bias the determination of the vertex po- 
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sition, the track with the largest vertex contribution 
greater than 6 is removed and the fit is redone until no 
track fails the requirement. In Monte Carlo simula- 
tion, the RMS of the core and tail Gaussian components 
of the residual Az distribution (measured Az minus true 
Az) is 190 /im. We fit this residual distribution to the 
sum of three Gaussian distributions and find that the 
RMS of the narrowest Gaussian, which contains 70% of 
the area, is about 100 /zm. Only 1% of the area is in the 
widest Gaussian. 

The absolute scale of the measurement of Az depends 
on the assumed positions of the silicon wafers in the 
SVT. These positions are determined from a combina- 
tion of survey measurements made before the SVT was 
installed and measurements of the positions of individ- 
ual SVT modules, each containing several silicon wafers, 
made with high-momentum charged particles that tra- 
verse the SVT. We check the absolute scale comparing 
the known positions of distinct mechanical features at 
each end of the beampipe (about 18 cm apart) with the 
apparent position measured with charged tracks in the 
SVT. The locations of these mechanical features are mea- 
sured from the positions of track vertices at least 2 cm 
from the interaction point that contain a well-identified 
proton, which are mainly due to e^-nucleon interactions 
in material. The measured distance in z between these 
mechanical features is in agreement with the known dis- 
tance to a precision of 0.2%. We conservatively enlarge 
this to 0.4% to account for any additional uncertainty in 
extrapolating to the interaction point. 



B. At measurement 

By far the dominant limitation on the accuracy with 
which At is determined from the measured decay length 
difference Az is the experimental resolution on the Az 
measurement. The next most significant limitation is the 
B meson momentum of about 340MeV/c in the T(4S') 
rest frame. We partially correct for this effect, as de- 
scribed below. The impact on the At measurement of 
the spread in the two beam energies, which results in a 
distribution of T(4S') momenta with a Gaussian width of 
about 6MeV/c, is negligible. Finally, we correct for the 
20 mrad angle between the T{AS) boost direction (the z 
axis in the following discussion) and the axis of symmetry 
of the detector, along which we measure the separation 
between vertices. 

Neglecting the B momentum in the T(4S') frame, we 
can write 



Az 



PjcAt, 



(20) 



where is the T(45') boost factor. The average value 
for the boost factor is (3j = 0.55. The boost factor is 
calculated directly from the beam energies, which are 
monitored every 5 seconds, and has an uncertainty of 
0.1%. 



In the case of a fully reconstructed -Broc, we measure 
with good precision the momentum direction of the re- 
constructed candidate, which can be used to correct for 
the B momentum in the T(45) frame. However, the cor- 
rection depends on the sum of the decay times, trcc + ^tag, 
which can only be determined with very poor resolution. 



We use the estimate t,-. 



tag 



tb -I- lAil to correct for 



the measured i?rcc momentum direction and extract At 
from the following expression: 

Az = /377,;ccAt -t- 7/3;ec7;cc cose;,,c{rB + \At\), (21) 

where 0*^^, P*^^, and j*^^ are the polar angle with respect 
to the beam direction, the velocity, and the boost factor 
of the Bj-cc in the T(45') frame. The difference between 
At calculated with Eq. |2^ and Eq. ^ is very small be- 
cause 7*(,(. = 1.002 and (3*^^ = 0.064. The event-by-event 
difference in At calculated with the two methods has an 
RMS of 0.20 ps. Equation ^ improves the At resolu- 
tion by about 5%. In addition, it removes a correlation 
between the resolution on At and the true value of At. 
This correlation is due to the fact that the RMS of the 
second term in Eq. ^ depends on the expectation value 
of (tree + itag)^, which In turn depends on |At|. Equa- 
tion ^ is used for all B decays to hadronic final states, 
while Eq. |o| is used for semileptonic modes since the B 
direction cannot be measured for these decays. 



C. Vertex quality requirements 

A number of requirements are made in order to ensure 
a well-determined vertex separation. The fit for both 
the By-ec and i?tag vertex must converge. Also, the error 
on At determined from the vertex fit must be less than 
2.4 ps and |At| must be less than 20 ps. The efficiency 
for passing these requirements in data and Monte Carlo 
simulation is about 97% for all i3roc modes. From the 
Monte Carlo simulation, we find that the reconstruction 
efficiency does not depend on the true value of At. 

The Bj-cc sample is used both to extract the B^-B^ 
mixing frequency and to measure the mistag probabilities 
for the analysis of time-dependent CP- violating asymme- 
tries. While the CP measurement is statistically limited, 
the mixing measurement has a statistical precision of a 
few percent. Therefore, in order to reduce systematic 
uncertainties in the mixing measurement, more restric- 
tive vertex criteria are imposed for the Siec sample used 
for the mixing measurement than for the CP and Bicc 
samples used for t he CP m easurement. However, as de- 
scribed in Section VIII C 5, the more restrictive criteria 
are also used as a cross-check in the CP measurement. In 
order to reduce further the contributions from charm de- 
cay products in the mixing analysis, a track is not used in 
the reconstruction of the -Stag vertex if it is identified as a 
kaon according to the kaon identification algorithm used 
for tagging (see Section IVA). The maximum allowed 



error on At determined from the vertex fit is decreased 
from 2.4 ps for the samples used for the CP measurement 
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FIG. 15: Schematic view of the geometry in the yz plane for a T(4S') — > BB decay. For fully reconstructed decay modes, 
the line of flight of the -Btag can be estimated from the (reverse) momentum vector and the vertex position of Brcc, and from 
the beam spot position in the xy plane and the T(45') average boost. Note that the scale in the y direction is substantially 
magnified compared to that in the z direction. 



to 1.4 ps for the sample used for the mixing measure- 
ment. The efficiency to pass these two additional criteria 
is about 87% in data. All figures in this section are ob- 
tained with the vertex selection criteria applied in the 
CP analysis. 



tions (called the core, tail and outlier components) with 
different means and widths: 



D. At resolution function 

The At resolution function is represented in terms of 
dt — At — Attruc by a sum of three Gaussian distribu- 
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{St - bkCTAt) 



exp 



(22) 



r 



For the core and tail Gaussians, we use the measurement 
error a At derived from the vertex fit for each event but 
allow two separate scale factors 5*1 and 5*2 to accommo- 
date an overall underestimate {Sk > 1) or overestimate 
{Sk < 1) of the errors for all events. Figure |l6|a illus- 
trates the correlation between the RMS of 6t and (JAt in 
Monte Carlo simulation. The core and tail Gaussian dis- 
tributions are allowed to have a nonzero mean to account 
for residual charm decay products included in the i^tag 
vertex. In the resolution function, these mean offsets are 



scaled by the event-by-event measurement error aAt to 
account for an observed correlation shown in Fig. ^6|b 
between the mean of the St distribution and the mea- 
surement error aAt in Monte Carlo simulation. This cor- 
relation is due to the fact that, in B decays, the vertex 
error ellipse for the D decay products is oriented with 
its major axis along the D flight direction, leading to a 
correlation between the D flight direction and the calcu- 
lated uncertainty on the vertex position in z for the -Btag 
candidate. In addition, the flight length of the D in the z 
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direction is correlated with its flight direction. Therefore, 
the bias in the measured Stag position due to inclusion 
of D decay products is correlated with the D flight di- 
rection. Taking into account these two correlations, we 
conclude that D mesons that have a flight direction per- 
pendicular to the z axis in the laboratory frame will have 
the best z resolution and will introduce the least bias in 
a measurement of the z position of the -Btag vertex, while 
D mesons that travel forward in the laboratory will have 
poorer z resolution and will introduce a larger bias in the 
measurement of the i?tae vertex. 
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FIG. 16: Correlation between the event-by-event error on At 
(cAt) and a) the observed RMS and b) offset of the mean for 
5t = /S.t — Attruc from Monte Carlo simulation. 

Monte Carlo simulations confirm the expectation that 
the resolution function is less biased for events with a pri- 
mary lepton tag than those with a kaon tag. Therefore, 
the mean of the core Gaussian is allowed to be different 
for each tagging category. One common mean is used for 
the tail component. The third Gaussian has a fixed width 
of 8ps and no offset; it accounts for the fewer than 1% of 
events with incorrectly reconstructed vertices. The reso- 
lution parameters extracted from the full likelihood fits 
to the At distributions are shown in Table XIII for the 



mixing analysis and in Table XVI for the CP analysis. 

Figure [itI shows the distribution of the uncertainties on 
At calculated from the fit to A2; for the flavor-eigenstate 
sample, and the combined -qcp — —1 {J^K^, tp{2S)K^, 
XciK^s) ^'^'i J/i'K*° samples, and compares data with 
Monte Carlo predictions. Since the Stag vertex preci- 
sion dominates the At resolution, no significant differ- 



ences between the At resolution function for the flavor- 
eigenstate sample and the CP-eigenstate sample are ex- 
pected. Hence, identical resolution functions are used for 
all modes. 
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FIG. 17: Distribution of event-by-event uncertainties on At 
((TAt) for a) the sample of neutral B decays to flavor eigen- 
states other than J/i}} K*^{K^n~) and b) the combined rjcp = 
-1 {J/iliKl, '4){2S)Kl, XciKl) and J/xpK*^ samples. The 
histogram corresponds to Monte Carlo simulation and the 
points with error bars to data. All distributions have been 
background-subtracted with events from the rriES sideband. 
The Monte Carlo distribution has been normalized to the 
same area as the data distribution. 



E. Checks and control samples 

Two of the fundamental assumptions in this analy- 
sis are that the At resolution function for the sample 
of flavor-eigenstate modes is the same as that for CP 
events, and that the event-by-event vertex errors provide 
a good measure of the relative uncertainty on the Az 
measurement from event to event. In this section, we 
describe several studies that have been done to validate 
these assumptions. We compare various distributions for 
the CP and flavor-eigenstate samples, in both data and 
Monte Carlo simulation. We take advantage of the small 
vertical size of the beam to measure the resolution for 
the B-ccc and the i?tag vertices in the vertical direction. 
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We also study the vertex resolution for a sample of D*"*" 
candidates from cc events and for continuum events in 
data and Monte Carlo simulation. 



1. Comparison of flavor-eigenstate and CP samples 

In Fig. ^ we compare various properties of the flavor- 
eigenstate sample with the combined r]cp = —1 and 
J/ipK*^ samples. These include the probability for 
the vertex fits, the number of tracks used in the Stag 
vertex, and the momentum in the T(4S') rest frame and 
polar angle in the laboratory frame of tracks used in the 
Stag vertex. Good agreement in all variables is observed 
between the two data samples. 

A similar comparison of the momentum and polar- 
angle distribution of tracks in data and Monte Carlo sim- 
ulation also shows good agreement. However, there are 
modest discrepancies for the probability for the ver- 
tex fits and the number of tracks used in the Stag ver- 
tex. The agreement improves when we include residual 
misalignments between the SVT silicon modules in the 
Monte Carlo simulation. Systematic uncertainties due 
to resid ual SVT mis alignme nts in data are discussed in 
Sections IVII B l| and IVIIIB. 



As expected, there are no significant differences ob- 
served in comparisons between the CP modes used in 
the sin2/3 analysis. However, comparisons between the 
CP and fiavor-eigenstate samples, in data as well as in 
the Monte Carlo simulation, show that the CP events 
have a slightly better Az resolution. For example, in 
Monte Carlo simulation the most probable value for a/^t 
is about 0.017 ps (3%) worse for the Sfiav sample, as can 
be seen by comparing the distributions in Figures p7| a and 
b. This is due to the fact that the Bcp vertex is better 
determined because tracks in the lower-multiplicity CP fi- 
nal states generally have higher momentum. We account 
for this effect in the likelihood fit by using the calcu- 
lated event-by-event errors, as described in Section VD. 
Indeed, for Monte Carlo simulation, the pull distribu- 
tions for CTAt (defined as the difference between the fitted 
and generated value divided by the calculated error) are 
nearly Gaussian with unit width for both the Bcp and 
Sfiav samples. Any residual effect due to differences in 
the observed scale factors in da ta is in cluded as a sys- 
tematic uncertainty (see Section VIII B ) and found to be 
negligible. 
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2. Vertex resolution in vertical direction 

Since the size of the PEP-II beam is only about 10 /Ltm 
in the vertical (y) direction, the measured distance Ay 
between the Sicc or Stag vertex and the nominal beam 
spot position in the y direction can be used to compare 
the resolution for the CP and flavor-eigenstate samples, 
and to evaluate the accuracy of the event-by-event errors 
cr(Ay). The average beam-spot y position is determined 



FIG. 18: Distributions of a) probability of the -Btag vertex 
fit, b) number of charged tracks and candidates used in the 
Stag vertex, c) momentum in the center-of-mass frame, and 
d) polar angle in the laboratory frame for tracks in the Stag 
vertex, for the flavor-eigenstate (histograms) and the com- 
bined ricp = — 1 and J/^ K*'^ (points with errors bars) data 
samples. All distributions have been background-subtracted 
with events from the mES sideband. The flavor-eigenstate 
distributions have been normalized to the same area as the 
distributions from the combined r]cp = — 1 and J/tp K*° sam- 
ples. 
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with a precision of better than a few microns with two- 
track events for each data run (approximately one hour 
of recorded data) . There is a non-neghgible contribution 
to Ay of ~ 25 /im (RMS) due to the B Ufctimc and the 
transverse momentum of the B. 

The distance in y between the i?tag vertex and the 
beam spot is used to measure the i?tag vertex resolution 
and bias in y. In Fig. |l^, we show the distribution of 
Ay/a{Ay) for the i?tag vertex for the flavor-eigenstate 
and CP samples, in data and Monte Carlo simulation. 
The RMS of the Ay/a{Ay) distribution is 1.3 and 1.4 
for Monte Carlo simulation and data, respectively. No 
statistically significant biases are observed. 

Similar results are obtained for the B^^c vertex reso- 
lution. In addition, good agreement in the resolution on 
the y position is observed between the flavor-eigenstate 
sample and the CP sample. The resolution is about 5% 
worse in data than in Monte Carlo simulation. 



3. Vertex resolution in continuum events 

Two samples have been used to cross-check the relia- 
bility of the resolution function extracted from the like- 
lihood fit as well as the discrepancies between data and 
Monte Carlo simulation: a sample of 109, 000 candi- 
dates from cc events and a sample of off-resonance data. 

For the first sample, we reconstruct high-momentum 
D*+ candidates in the mode D*+ D"tt+, followed by 
Z?" K~7T~^, if~7r+7r°, or i^~7r+7r~7r+, and then use 
the remainder of the charged tracks in the event (frag- 
mentation particles and recoil charm decay products) to 
determine a vertex position with the standard Stag ver- 
tex algorithm. Since position information for the D*'^ 
vertex is poor, due to scattering of the slow pion, and 
the D*~^ decay point coincides with the e'^e~ interaction 
point, a beam-spot constraint is used for the D*^. In 
Fig. we show the distribution of the distance along 
the z axis between the D*^ vertex and the vertex formed 
from the rest of the tracks in the event Az, as well as Az 
divided by the event-by-event error on Az, for both data 
and Monte Carlo simulation. 

In Monte Carlo simulation, the resolution on z for the 
candidate is « 90 /xm, very similar to that for B^cc 
vertices. However, the momentum spectrum of fragmen- 
tation tracks in cc events is softer than that for tracks 
from B decays, while D mesons are more energetic in 
the D*~^ control sample than in B decays. Therefore, 
a slightly more asymmetric resolution function is ex- 
pected for the D*^ control sample compared to that for 
B events, as shown in Fig. |2C|a. Comparison of distribu- 
tions of several sensitive variables (such as the number of 
tracks used in the vertex, and the momentum and polar 
angle of the tracks) shows small differences between Z)*+ 
and B events. 

The RMS of the distance between the Z?*+ vertex and 
the vertex formed from the rest of the tracks in the events 
is about 220 /im in the Monte Carlo simulation. Fitting 




FIG. 19: Distributions of the measured distance in the verti- 
cal direction Ay between the -Btag vertex and the beam spot 
position, divided by the event-by-event error on the mea- 
sured distance o-{Ay) for each event: a) Monte Carlo sim- 
ulation (histogram) and data (points with error bars) for the 
flavor-eigenstate sample; b) flavor-eigenstate (histogram), and 
r]cp = — 1 and J/tpK*^ (points with error bars) samples in 
data. All distributions have been background-subtracted with 
events from the mES sideband. In a), the data distribution 
has been normalized to the same area as the Monte Carlo 
simulation distribution; in b) the combined t^cp = — 1 and 
J/ip K*" data distribution has been normalized to the same 
area as the flavor-eigenstate distribution. 

this distribution to the sum of three Gaussians, we find 
a resolution of about 140 fim for 97% of the events, com- 
pared to 150 /im for 99% of the B^B^ events. Only small 
differences are observed in the distribution of Az/(t(Az), 
as illustrated in Fig. |20|b. Therefore, the sample can be 
used to confirm the resolution and scale factors extracted 
from the likelihood fit, as well as to compare how well the 
Monte Carlo simulation reproduces the data. 

The distributions are fit to the sum of three Gaussian 
distributions with different widths and means. The width 
of the third Gaussian is fixed to 2.0 mm. From the fit 
results, we come to the following conclusions: 

• The event-by-event errors on Az are underesti- 
mated by about 10% in data (Fig. po|b). 

• The bias in the resolution function due to charm 
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decay products that is observed in data is well re- 
produced by the Monte Carlo simulation as shown 
in Fig. pO|a and b. 



• The resolution measured in the data is about 5% 
worse than that predicted by the Monte Carlo sim- 
ulation (Fig. 

These results will be compared with those obtained from 
the likelihood fit to the B events, described in Sec- 
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FIG. 20: a) Az and b) Az/cr{Az) distributions for the D*+ 
control sample in data (points with error bars) and Monte 
Carlo simulation (solid histogram). For comparison, the dif- 
ference between the measured Az and true Az for rjcp = — 1 
and J/'^K*^ events in Monte Carlo simulation is also shown 
(dashed histogram). All distributions are normalized to unit 
area. 



The second control sample is obtained from off- 
resonance data alone. Charged tracks from these con- 
tinuum events are randomly split into two sets, and the 
vertex of each set is found with the same algorithm used 
to determine the i?tag vertex. In this case the Stag ver- 
tex reconstruction strategy is applied to both vertices 
in the event, so that this sample provides an unbiased 
estimation of the resolution, suitable for comparisons be- 
tween data and Monte Carlo simulation. Results from 
this study are compatible with those reported above. 



F. Comparison of 1999-2000 and 2001 performance 

The internal alignment of the SVT has improved sig- 
nificantly for the reconstruction of the 2001 data set (Run 
2) compared to 1999-2000 (Run 1). Therefore, we expect 
better resolution and event-by-event errors on Ai for Run 
2, which requires the use of separate resolution functions 
for the two data sets. 

The differences in resolution and event-by-event errors 
for Run 1 and Run 2 are illustrated in Fig. where 
a comparison of the distributions for Az and Az/ct(Az) 



in the control sample described in Section VE 3 
shown. From the separate analysis of the two data sets, 
we conclude the following: 

• The event-by-event errors on Az are underesti- 
mated by 15% (5%) for the Run 1 (Run 2) data 
set (Fig. |l|b). 

• There is no statistically significant difference in the 
bias between the two data sets. 

• The resolution for the Run 1 data set is about 15% 
worse than for that for Run 2 (Fig. ^T|a). 
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FIG. 21: Comparison of the distributions of a) Az and 
b) Az/a(A z) for the D*'^ control sample described in Sec- 
tion V E 3| , for Run 1 (points) and Run 2 (histogram) data. 
All distributions are normalized to unit area. 
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The improved quality of the event-by-event errors in 
Run 2 is also illustrated in Fig. where we compare 
the distributions of probability for the i?tag vertex fit 
with the flavor-eigenstate data sample selected from the 
two different data periods. 
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FIG. 22: Comparison of the x probability distributions of 
the Btag vertex fit for the flavor-eigenstate data samples in 
Run 1 and Run 2. The distributions have been background- 
subtracted with events from the tties sideband. The area of 
each distribution equals the total number of events in the 
corresponding sample. 



VI. LIKELIHOOD FIT METHOD 

The value of sin2/3 is extracted from the tagged Bcp 
sample with an unbinned maximum-likelihood technique 
based on In Ccp and the probability density functions J-± 
of Eq. However, the dilutions Vi and Az resolution 
parameters are also needed for the measurement. As- 
suming that mistag rates and vertex resolutions do not 
depend on the particular channel used to reconstruct the 
B meson, these parameters are best determined with the 
much larger mixing sample, since they also appear in 
Cmix- In order to properly incorporate the correlations 
between these parameters and sin2/?, the ht is performed 
by simultaneously maximizing the sum 



ln£ 



CP 



In Cn 



(23) 



for the combined tagged Bflav and Bcp samples. The 
values of B^ lifetime and Am^ are kept hxed in extracting 
sin2/3. 

The value of Arud is obtained with an unbinned 
maximum-likelihood ht to the tagged Bfiav sample alone, 
where the log-likelihood In £mix is maximized while keep- 
ing the B^ lifetime hxed. 



example, the response of the detector to positive pions 
and kaons differs from its response to negative pions and 
kaons due to differences in total and charge- exchange 
cross sections. To account for any possible mistag dif- 
ferences, we introduce separate mistag probabilities w 
for B^ and w for with the conventions 



V 

{V) 



1/ -^ 

-(w + w); 
1 - 2w] 



Aw = 


[w 


— w) 


V = 


1 


- 2w 


AV ^ 


{V 


-V) 



The time distributions for the mixing and CP samples 
will thus depend on whether the tag was identified as a 
or a S*^ , resulting in modihcations to the expressions 
for mixing time development (Eq. 0) 



h±,taa=BO oc [{I + -AV) ± {V) cos ArudAt] 
h 



±M9=B« « [(l--AP)±(I?)cosAmdAt], (24) 



where the ± in the index refers to mixed (— ) and un- 
mixed (-I-) events as before, and for CP violation time 
development (Eq. 0) 

/± cx [(1 ± ]^AV) T {V)ricp sin2/3sin Am^Ai], (25) 

where the ± in the index refers to events where Stag is a 
5° {+) and B° (-) and we have taken |A| = 1. 

B. Background modeling 

In the presence of backgrounds, the probability distri- 
bution functions 7i± of Eq. ^ and J-± of Eq. |l^ must 
be extended to include a term for each significant back- 
ground source. The backgrounds for the flavor eigen- 
states and rjcp ~ —1 modes are quite small and are 
mostly combinatoric in nature. However, for the B'^ —> 
J/'ipKl and B^ — > J/ipK*'^ channels the backgrounds are 
substantial and originate mainly from other B — > J/il^X 
modes that have, to a very good approximation, the same 
flavor tagging and At resolution properties as the sig- 
nal. The background properties of the flavor eigenstates, 
rjcp — ~\ modes, and the non-J/ip background in the 
B^ —> J/ipKl channel are determined empirically from 
sideband events in the data. 



1. Background formulation for flavor eigenstates and 
r]cp = — 1 modes 



A. Mistag asymmetries 

The probabilities of mistagging a 5° or meson are 
expected to be very nearly, but not exactly, equal. For 



The background parameterizations are allowed to dif- 
fer for each tagging category. Each event belongs to a 
particular tagging category i. In addition, the event is 
classified as either mixed (— ) or unmixed (-I-) for a fiavor- 
eigenstate or by whether i?tag was a B'^ {+) or a S'^ (— ) 
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for a CP-eigenstate. Thus background distributions j the fuU hkehhood function becomes 
must be specified for each possibihty (+/— so that 



W±,, = /5;^±(Ai;r,Amrf,ii;,,aO + /!!p^:ak^ilpoak(Ai;a,)+ ^ /«f 6,) (26) 

j=bkgd 



for flavor-eigenstates, and 



•^±,. - /Sg-^±(A<; r, Ama, sin2/3, a,) + f^P,,^Bg%.A^t- a,) + ^ /^^Sg^^.(Ai; 6,) (27) 

j=bkgd 



for CP-cigenstates. 



The fraction of background events for each source and 
tagging category is a function of rngg and is given by 
fi^j. The peaking and combinatorial background PDFs, 
B±,i,i,eak and B±^i^j, provide an empirical description of 
the At distribution of the background events in the sam- 
ple, including a resolution function parameterized by di 
and bi, respectively. These distributions are normalized 
such that, for each i and j, 



B-,^j)dAt = 1. 



(28) 



The probability that a candidate is a signal or a 
background event is determined from a separate fit to the 
observed mEs distributions of Bflav or Bcp candidates 
with ricp — —1- We describe the mgs shape with a 
single Gaussian distribution 5(mEs) for the signal and 
an ARGUS parameterization ^(tties) for the background 
(Eq. Based on this fit, the event-by-event signal 

and background probabilities that appear as the relative 
wei ght s for the various signal and background terms in 
Eq.]2^ and ^ are given by 



/i,sig(mEs) 

/j,pcak("T-Es) 

X! /ij("*Es) 
j=bkgd 



(1 ~ Spcak)S{mEs) 

S{mEs) + -4(mEs) 

5(toes) + ^(toes) 

-4(toes) 
5(toes) + A{mEs) 



(29) 



The fraction Jpoak of the signal Gaussian distribution 
that is due to peaking backgrounds is determined from 
Monte Carlo simulation. 

Backgrounds arise from many different sources. Rather 
than attempting to determine the various physics con- 
tributions we use an empirical description in the likeli- 
hood fit, allowing for background components with var- 
ious time dependencies. For the Bflav sample, the back- 
ground time distributions considered, each with its own 
effective dilution factor Vi and either a common resolu- 
tion function TZ{At]bi) or the signal resolution function 



TZ{St = At ~ Attruo; ai), are 



Bi^li = i(l±I?^r) ^(A^true)®7^(<5t;6, 



K?flav 
lb, 2, peak 



irfj(i±i?«r)x 



-rfj|Att„o| 



«)7^((5t;6i), 



rf J(l ± cos Am,,3 Attruo) X 



-rfj|Att,uo 



\ ri^pcak(l ± T^f^pcak COS Awi^peak A^truc) 



tpe.J^*--l®7^(5t;aO, 



(30) 



corresponding to prompt, non-prompt, and mixing back- 
ground components, as well as a peaking contribution. 
For the 7]cp = — 1 sample, the possible background con- 
tributions are 



V.CP 



kCP 



j^CP 

^±,2, peak 



-<5(A^)®7^(5t;&^), 

i ri,2(l ± Vf^ sin AmrfAt) x 



-^??|At|^7^(,5t;M, 

r,,peak(l±2?fpeakSinAmrfAt) 



Atl 



(31) 



corresponding to prompt and CP background compo- 
nents, as well as a peaking contribution. The background 
resolution function parameters bi are common with the 
background resolution function of the i?flav sample. The 
likelihood fit includes as free parameters the fraction of 
each time component, as well as apparent lifetimes, res- 
olutions, mixing frequencies and dilutions that best de- 
scribe the events with high weights for being background. 
These parameters are described in Section VI D| below, 
along with additional assumptions. 
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TABLE IX: Parameters of the probability distribution func- 
tion for the non-J/tj) background contribution in the 
J/'tpK'i channel. 



As discussed in Section HID, the data are used to de- 



Parameter 


Fit result 




0.16 ±0.49 




1.25 ±0.45 


Si 


1.12 ±0.26 


bx 


-0.11 ±0.20 


S2 


3.9 ±0.8 


b2 


-1.2 ± 1.0 


/2 


0.23 ±0.14 


/3 


0.005 (fixed) 



2. Background formulation for B° — > J/i/j Kl 



termine the relative fraction of signal, background from 
B — > J/^X events, and events with a misreconstructed 
J/'0 — > U candidate. Along with a Monte Carlo simula- 
tion of the channels that contribute to the B — + J/tpX 
background, this information is used to formulate the 
PDF model. In addition, some of the J/tpX background 
modes, such as B" J/tpK*^ and -> JZ-^if", have a 
non-zero CP asymmetry (ijcp), as given in Table VI. The 
value of the asymmetry in i?° J/tjj K**^ {K'^tt'^) is taken 
from the measurement of Rj_ = 0.160 ± 0.032 ± 0.014 
in Ref. jl3|. The probability density functions T± of 
Eq. ^ are modified to include contributions for each of 
the B — > J/i/j X channels a specified in Table VI and the 
non-J/ip background component. The complete PDF is 
given by 



The higher background level in the i?'^ J/tp chan- 
nel requires a more extensive treatment of its properties. 



^±,z = /*,fc,sig(A£;).F±(Ai;r,Amd,w;„sin2^,a,) ± ^ f.^kA^E) T±{At;T, Anid J Vcp,ay sin2/3, Ui) 

a=J/ipX 

+ /^,fc,no„-j/^(A^) S±''(At;S). (32) 
I 



Each event is classified according to its flavor tagging 
category (i), flavor tag value (±), and the recon- 
struction category (fc), which is either EMC or IFR. The 
signal fraction /i,fc,sig and background fractions fi^k,a and 



-J/->P 



are determined as a function of AE and are 



the same for all tagging categories. The shape of the sig- 
nal and background AE functions are determined either 
from data (non-J/ip contribution) or from Monte Carlo 
samples (signal and J/tpX background). The normaliza- 



tion s f^ ii^^y^f / d{AE) are determined from Tables VI 
and [VI I| so that 

"lOMoV 



lOMcV 



/,,fe..ig(Ai?)±/,,fc 

,non— J/V' 

(A£;)± 

f,,kA^E)]diAE) ^1. (33) 



E 



The non-J/ip background PDF is given by 



Bf^ = F.=oS^-i±(l-F.=o)S-- 



CP 



(34) 



where the dilutions I'pl' — and the prompt fraction 
Fr^o, effective decay width ri^2, and At resolutions pa- 
rameters b are fixed to values obtained from an external 
fit to the m{U) sideband events as given in Table IX. 
The resolution function TZ{At; b) is defined in Eq. |2^ with 
/a = 0.005 and with core bias 5i equal for all tagging cat- 
egories. 



The J/tpK^ sample has significant background, pri- 
marily from other J/ip modes. The Monte Carlo sim- 
ulation is used to check the flavor tagging efficiency of 
the inclusive J/^ background relative to the signal for 
the Kl mode. The inclusive J/ip background fraction in 
the simulation is consistent across the flavor tagging cat- 
egories to within a few percent. The flavor tagging effl- 
ciency for the fake- J/V' background, determined from the 
J/ip sideband, is also roughly consistent with signal. The 
composition of the J/ipK^^ sample is determined from a 
flt of the AE spectrum before flavor tagging. We assume 
the inclusive J/ijj and fake- J/?/; background fractions are 
independent of flavor tag in the nominal fit and adjust 
the fractions as a function of tagging category, based on 
the Monte Carlo simulation and J/ip sideband, in order 
to determine systematic errors. 

Some of the decay modes in the inclusive J/tj} back- 
ground, such as J/^K*^ and J/tpK^, have an expected 
CP asymmetry. The mistag fractions for all CP modes 
in the inclusive J/ip background are determined with the 
Monte Carlo simulation and found to be consistent with 
the values for the signal. We assume that the signal 
mistag fractions apply to the CP modes in the inclusive 
J/^p background. 

The At resolution for the B J/ipX background 
should be very similar to the signal resolution. However, 
extra tracks associated with _B+ J/ipX^ decay, such 



36 



as the charged tt from the K*'^ decay in J/tpK*~^, 
could bias the measurement of At since they are not as- 
sociated with the Bcp vertex and therefore can be used 
in the -Btag vertex. In the Monte Carlo simulation, we 
find that extra tracks in the B —> J/tpX decay modes 
have a negligible effect on the At resolution. Therefore, 
we assume that all B — > J/4> X background has the same 
resolution as the signal. 

The At resolution of the non- J/ijj background is mea- 
sured with the J/tp sideband sample. The noii-J/ip At 
resolution parameters are varied by their statistical un- 
certainties to estimate the systematic uncertainty. 



channel. As shown in Table the background for this 
channel is due to true B — > J/ipX decays. Thus, we as- 
sume that the background has the same resolution func- 
tion and tagging performance as the signal. The prob- 
ability density functions J-± of Eq. ^ are modified to 
include contributions for each of the B J/ipX chan- 
nels a specified in Table 0. The complete PDF is given 

by 



3. Background formulation for B° J/tp K'°{K°sTt°) 

Monte Carlo simulation is used to construct the prob- 
ability density function for the S° J/V'X*"(X°7r°) 



^±A = /sig.f±(A<; r, Am<j, 77cp,signai, m, sin2/3, a^) 



fa ^±{At;r,Amd,rif^a,Wi,s'm2/3,ai), (35) 

a— bkgd 



r 



where each event is classified according to its fiavor tag- 
ging category (z) and flavor tag value (±). The signal 
and background fractions as well as rjcp are taken from 
Table 0. 



C. Extensions for direct CP search 

While the main likelihood fits are performed with the 
Standard Model expectation that |A| — 1, a search for the 
effects of direct CP violation is also made. Such a mea- 
surement is also particularly sensitive to possible differ- 



ences in the fraction of S*^ or B'^ meson that are tagged. 
Defining etag and etag as the tagging efficiencies for B^ 
and B^ , and and e^- as the reconstruction efficiencies for 
and in the i?flav sample, it is useful to construct 



eta 



eta 



(eta 



eta 



eta 



(36) 
(37) 



For the Bcp sample, the time-dependent decay rate 
(Eq. P) becomes 



J 



/±,(At) = ie-i^l^*l 



{l + \Xm±X,) 



2ImX 



i + |A|' + e, 
I 



sin Amrf At 



1-|A|= 



i + iAr+e 



cos AnidAt 



, (38) 



where = ^,(|A|2 - 1)/(1 + xj), X, = fi,{V), + AVJ2, where 
and X^ — {'D)i+ fiiA'Di/2. Likewise, for the -Bflav sample 
the time-dependent decay rate (Eq. |^) becomes 

si = 1(-1) if the reconstructed B is a B"{B°) 

S2 = 1(-1) for a B°{B°) tag. 



h{At) 



1 + xl 



S2Xi 



The parameters Vi, (etag)i, and can be extracted 
si(/^i + S2X[) cos ATO^At] (39) from time-integrated numbers of events in the -Bflav sam- 
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pie. Defining integrated samples of events by 

TVfs ^ Ar(BO/B° tag in category, Bflav = B") 
= N{B"/B" tag in i*'' category, Ba^^ = B°) 
= iV(no tag in i^^ category, Busn = B°) 
—no tag ^ ^^^^ category, Bflav = (40) 



it can be shown that 



iV^8 - TV ' 



no tag 



TV 



(eta 



jyno tag 

2iVfs7vrViVr^Af, 



^no tag 



(41) 



—tag ,^„o tag _|_ jytagjy-no tag 



Mi 



2(iVf*5 + iVf° t-g)(7V, 
(1 + a;2)(iVr^A^r - 



2^tag^tag ^tag tag _^ ^tag^no tag 



under the assumption that nearly all B mesons decay to 
final states that can be reached from either B'^ or 
but not both. The results for (etag)i and are shown in 
Table The value of Vi, averaged over all four tagging 
categories, is 0.004 ± 0.012. While there is no statisti- 
cally significant difference in the tagging efficiencies or 
the reconstruction efficiencies given by fii and Vi, we use 
the central values obtained from the i?flav sample in per- 
forming the fit for |A|. 



TABLE X: Values of {etag)i and fii for the four tagging cat- 
egories, as determined by counting numbers of tagged and 
untagged events in the Bflav sample. 



Tagging category 






Lepton 
Kaon 
NTl 
NT2 


0.095 ± 0.002 
0.358 ± 0.003 
0.080 ±0.002 
0.139 ±0.002 


0.069 ± 0.032 
-0.005 ±0.014 
0.061 ± 0.035 
0.017 ±0.026 



TABLE XI: Average mistag fractions {wi) and mistag dif- 
ferences Awi for each tagging category i from a maximum- 
likelihood fit to the distribution for the control sample. 



Tagging category 


(w). [%] 


Aw^ [%] 


Lepton 
Kaon 
NTl 
NT2 


4.6 ±0.6 
11.8 ±0.5 
21.3 ± 1.6 
37.2 ± 1.3 


1.1 ± 1.2 
-0.3 ± 1.0 
-5.9 ±3.2 
-0.7 ±2.7 



D. Free parameters for the sin2/3 and Arud fits 

The unbinned likelihood fit for sin2/3 has a total of 45 
free parameters: 

• Value of sin2/3 



• Signal resolution function: Sixteen parameters 
flj to describe the resolution function for the sig- 
nal. Due to improvements in the reconstruction 
algorithms, the Run 1 and Run 2 resolution func- 
tions are found to be different, as described in Sec- 
tion VT. Thus, wc allow for separate resolution 
function parameters for these data samples, each 
with eight free parameters, being a scale factor 
Si for the event-by-event Az resolution errors of 
the core Gaussian components, individual core bias 
scale factors bi i for the four tagging categories and 
a common tail bias 62, and the tail /2 and outlier /a 
fractions; the scale factor of the tail component is 
fixed to 3.0 and the width of the outlier component 
is fixed to 8 ps with zero bias. 

• Signal dilutions: Eight parameters to describe 
the measured average dilutions {'D)i and dilution 
differences AVi in each tagging category. 

• Background resolution function: Six parame- 
ters are used to describe a common resolution func- 
tion for all non-peaking backgrounds. As with the 
signal resolution function, we include separate res- 
olution function parameters for the Run 1 and Run 
2 data samples. The resolution function is taken as 
a single Gaussian distribution with a scale factor 
5*1 for the event-by-event Az errors and a common 
bias scale factor 61, and an outlier fraction f^; the 
width of the outlier component is taken to be a 
fixed 8ps with zero bias. 

• -Bflav background composition parameters: 

A total of 13 parameters describe the Bflav back- 
ground composition. We make several assumptions 
to simplify the parameterization shown in Eq. |3^, 
such as removing the mixing background contribu- 
tion by setting /^^|^ — 0, and assign a correspond- 
ing systematic uncertainty. The size of the peaking 
background is determined from Monte Carlo simu- 
lation to be S^^l^ = (1.5±0.5)% of the signal contri- 
bution in each tagging category. This contribution 
is predominately from B'^ events, so Ami peak = 0, 
rf^^cak = and D^l^^ are taken from the B+ 
data sample (Table |Xl| ) . The effective dilutions for 
the prompt {Dfi^, 4 parameters) and non-prompt 
i-^i.T^ 4 parameters) contributions are allowed to 
vary. The relative amount of these two contribu- 
tions is allowed to vary independently in each tag- 
ging category (4 parameters). For the non-prompt 
contribution, Tf^ is assumed to be the same for all 
tagging categories, giving one free parameter. 

• CP background composition parameters: 

One parameter, the fraction of prompt relative to 
non-prompt background, assumed to be the same 
for each tagging category, is allowed to float to de- 
scribe the CP background properties. The effective 
dilutions of the non-prompt and peaking contribu- 
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tion are set to zero [0^2 = -Cfp^ak = 0), corre- 
sponding to no CP-asymmetry in the background. 
The size and parameters of the peaking background 
are determined from Monte Carlo simulation. The 
fraction of peaking background is 5'^^^-^ — (1 ± 1)% 
of the signal contribution, independent of tagging 
category. This contribution is assumed to have life- 
time parameters in common with the signal. Fi- 
nally, the lifetime of the non-prompt background 
is assumed to be r^o in all tagging categories. 

The unbinncd likelihood fit for Am^ has 44 free param- 
eters, removing sin2/3 and the parameter for fraction of 
prompt background in the CP sample and leaving Am^ 
to float. 



E. Blind analysis 

A blind analysis technique was adopted for the extrac- 
tion of sin2/3 and Am^; in order to eliminate possible ex- 
perimenter's bias. We used a method that hides not only 
the central value for these parameters from the unbinned 
maximum- likelihood fit, but also the visual CP asym- 
metry in the At distribution. The error on both the 
asymmetry and Am^ is not hidden. 

The amplitude of the asymmetry Acp{^t) from the 
fit was hidden by a one-time choice of sign flip and arbi- 
trary offset based on a user-specified key word. The sign 
fiip hides whether a change in the analysis increases or 
decreases the resulting asymmetry. However, the magni- 
tude of the change is not hidden. The visual CP asym- 
metry in the Ai distribution is hidden by multiplying At 
by the sign of the tag and adding an arbitrary offset. 

With these techniques, systematic studies can be per- 
formed while keeping the numerical value of sin2/3 or 
Amd hidden. In particular, we can check that the hid- 
den At distributions are consistent for and tagged 
events. The same is true for all the other checks con- 
cerning tagging, vertex resolution and the correlations 
between them. For instance, fit results in the different 
tagging categories can be compared to each other, since 
each fit is hidden in the same way. The analysis proce- 
dure for extracting sin2/3 and Am^ were frozen prior to 
unblinding. 



VII. B° FLAVOR OSCILLATIONS AND 
MISTAG RATES 

A. Likelihood fit results for A.md 

We extract Arrid, the dilution factors T>i, the At res- 
olution parameters a^, and the background At param- 
eterization by fitting the At distributions of the flavor- 
eigenstate B^ sample with tties > 5.2GeV/c^ with the 
likelihood function described in Section VI . The selection 
of the -Bfiav data sample is described in Section [II C 1 . 



We also demand a valid tag and At determination for the 
event, based on the algorithms described in Sections IV 
and ^ The more restrictive requirements |At| < 20 ps 
and (TAt < 1.4 ps are applied to the proper time-difference 
measurement. In addition, identified kaons in the -Stag 
decay are rejected in the reconstruction of the tagging 
vertex. These requirements are intended to reduce sys- 
tematic errors on the precision Arud measurement. The 
final sample consists of 12310 fully-reconstructed and 
tagged candidates with toes > 5.2GeV/c^, of which 
7399 are in the signal region rngs > 5.27GeV/c^. 

The breakdown of this mixing i?fiav sample into indi- 
vidual tagging categories is shown in Fig. ^ as a function 
of TOES- Superimposed on the observed mass spectra are 
the results of the fits with a Gaussian distribution for the 
signal and the ARGUS background function for the back- 
ground. The tagging efficiency and signal purity for the 
individual tagging categories in data are extracted from 
fits to the Toes distributions and are listed in Table XII . 
The efficiency for each tagging category is defined as the 
ratio of the number of signal events for each tag over the 
total number of signal events after imposition of vertex 
cuts. 



TABLE XII: Tagging efficiencies for hadronic decays and 
signal purities in data, shown separately for the four tag- 
ging categories. Signal purities are estimated for tties > 
5.27GeV/c^ 



Tagging 
Category 


Efficiency 

[%] 


Signal 


Purity 

[%] 


Lepton 
Kaon 
NTl 
NT2 


11.8 ±0.3 

33.9 ±0.5 
8.6 ±0.3 
13.9 ±0.4 


1097 ± 34 
3156 ±63 
798 ± 31 
1293 ±43 


96.0 ±0.7 
84.6 ±0.7 
88.9 ± 1.2 
79.4 ± 1.3 


Full sample 


68.1 ± 0.4 


6347 ± 89 


85.8 ±0.5 



The results from the likelih ood fi t to the mixing sam- 
ple are summarized in Table XIII. The probability to 
obtain a likelihood smaller than the observed value, eval- 
uated with fast parameterized Monte Carlo simulation 
of a large number of similar experiments, is (44 ± 1)%. 
The At distributions of the signal (toes > 5.27GeV/c^) 
and background (toes < 5.27GeV/c^) candidates, over- 
laid with the projection of the likelihood fit, are shown 
in Fig. ^ In Fig. ^ the mixing asymmetry of Eq. ^ is 
plotted; the time-dependence of the mixing probability 
is clearly visible. 

The tagging separation Q = etagi^ — 2w)'^ is calcu- 
lated from the efficiencies and the mistag rates quoted in 
Tables XII and XIII respectively. Summing over all tag- 



ging categories, we measure a combined effective tagging 
efficiency Q « 27%. 

Two small corrections, which are described in more 



detail in Sections VII B 2 and VII B 4 together with their 
assigned systematic errors, are applied to the output of 
the fit. The value of Ato^ obtained after applying these 
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FIG. 23: Distribution of niES for mixing Bflav candidates in 
separate tagging categories (Lepton, Kaon, NTl and NT2), over- 
laid with the result of a fit with a Gaussian distribution for 
the signal and an ARGUS function for the background. 



corrections is 

Amrf = 0.516 ± 0.016 ± 0.010 ps^\ 

where the first error is statistical and the second system- 
atic. 

We have also examined the fitted value for Ato^ with 



various subsamples of the full data set, including indi- 
vidual B decay channels, separate tagging categories, the 
state of the reconstructed Brec or tagging Stag , and dif- 



ferent time periods. As can be seen from Table XIV , the 
values obtained from the subsample fits are all consistent 
with the global result for Am^. 



B. Systematic error estimation 

Systematic errors can be grouped into four categories: 
signal properties and description, background properties 
and description, fixed external parameters and statistical 
limitations of Monte Carlo validation tests of the fitting 
procedure. A summary of the se so urces for the hadronic 
sample is shown in Table XV. In the following, the 
individual contributions arc referenced by the lettered 
lines in this table. 



1. Signal properties and description 

For the signal events, the use of a double Gaussian plus 
outlier model for re-scaling the event-by-event At errors 
as part of the likelihood fit means that uncertainties in 
the vertex resolution are incorporated into the statistical 
error on Am^, including proper treatment of all corre- 
lations. Assuming that this model is sufficiently flexi- 
ble to accommodate the observed distribution in data, 
no additional systematic error need be assigned. The 
contribution to the total statistical error due to the ver- 
tex resolution can be extracted by fitting the data twice: 
once holding all parameters except Am^ fixed, and once 
allowing the resolution function parameters to vary in ad- 
dition to Am^. Subtracting in quadrature the respective 
errors on Anid from the two fits shows that ±0.005 ps~^ 
of the statistical error can be attributed to the resolution 
parameters. 

To determine the systematic error due to the assumed 
parameterization of the resolution model, we apply a 
number of possible misalignment scenarios to a sample 
of simulated events. By comparing the value of Arn^ 
derived from these misaligned samples to the case of 
perfect alignment, we derive a systematic uncertainty of 
±0.004 ps-i (a). 

An additional systematic error is attributed to uncer- 
tainties in the treatment of the small fraction of At out- 
liers that are the result of misreconstructed vertices. The 
stability of the Arud result is examined under variation of 
the width of the third Gaussian component in the resolu- 
tion function between 6 and 18 ps, and through replace- 
ment of the third Gaussian with a uniform distribution 
and varying the width between 8 and 40 ps. On this ba- 
sis, we attribute a systematic uncertainty of ±0.002 ps~^ 
to the outlier treatment (b). 

As described in detail in Section VA, the beamspot 



position is an integral part of the determination of At. 
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FIG. 24: Distributions of At for unmixed (upper panel) and mixed (lower panel) events in the hadronic B sample, divided 
into a signal region tues > 5.27GeV/c^ with a) a linear and c) logarithmic scale, and a sideband region tties < 5.27GeV/c'^ 
with b) a linear and d) logarithmic scale. In all cases, the data points are overlaid with the result from the global unbinned 
likelihood fit, projected on the basis of the individual signal and background probabilities, and event-by-event At resolutions, 
for candidates in the respective samples. In a) and c), the At distributions obtained from the likelihood fit to the full sample 
are overlaid, along with the simultaneously determined background distribution shown as the curve in b) and d). 



Increasing its vertical size by up to 80 /im, and system- 
atically biasing its vertical position by up to 80 ^m, re- 
sults in a corresponding variation of Am^ by less than 
0.001 ps-i (c). 



The requirement on the maximum allowed value of a At 
is varied between 1 and 2.4 ps, and the observed variation 
of 0.003 ps~^ in Am^ is assigned as a systematic uncer- 
tainty (d). The observed dependence is mainly due to 
correlations between tagging and vertexing, as described 
in Sec. IVIIB4 



2. Background properties 

A systematic uncertainty in Am^ arises from our abil- 
ity to separate signal from background as a function of 
niEs- We estimate this uncertainty by varying the width 
and height of the fitted Gaussian peak in tties , the slope 
parameter of the ARGUS background shape, and the nor- 
malizations of the signal and backgrounds by one stan- 
dard deviation around their central values, resulting in 
an uncertainty of ±0.002 ps""'^ in Am^ (e). 

the At distribution of the 



As discussed in Sec. VI B 1 



background is described by the combination of a prompt 
component and a lifetime component. To estimate the 
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FIG. 25: Time-dependent asymmetry A{At) between unmixed and mixed events for hadronic B candidates with rriES > 
5.27GeV/c^, a) as a function of At; and b) folded as a function of \At\. The asymmetry in a) is due to the fitted bias in the 
At resolution function. 



systematic uncertainty due to this choice, we add an ad- 
ditional component, with its own separate hfetime, that 
is allowed to mix; the observed value of Ato^ changes 
by 0.001 ps~^ (f). Similarly, adding an additional Gaus- 
sian distribution to the At background resolution func- 
tion changes Anid by no more than 0.001 ps~^ (g). 

Finally, the composition of the background changes 
slightly as a function of tties, since the fraction of back- 
ground due to continuum production slowly decreases 
towards the B mass. As a result, the At structure of 
the background could change as well. To study this de- 
pendence, we split the niES sideband region into seven 
mutually exclusive, lOMeV/c^-wide intervals, and repeat 
the Arrid fit with each of these slices in turn. The vari- 
ation of Anid is then extrapolated as a function of the 
position of the sideband slice relative to the B mass. We 
correct the value of Amd by —0.002 ps~^ obtained from 
this extrapolation, and assign the statistical uncertainty 
of 0.002 ps^^ of this procedure as a systematic error on 
Anid (h). 

A small fraction (about 1.5%) of the events attributed 
to the B^ signal by the fit to the tties distribution con- 
sists of B^ events, mainly due to the swapping o f a soft 



with a charged pion as described in Section [II C 1 



The uncertainty on this peaking fraction is propagated 
to Amd, yielding a systematic error of 0.002 ps~^ (i). 



3. External parameters 
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FIG. 26: a) Correlation between the event-by-event error on 
At {(7 At) and the mistag rate in the Kaon category from Monte 
Carlo simulation; b) Dependence of mistag rate on aAt after 
scaling the mistag rate by \/'}2Pt- 



An error in the boost of the T{AS) system (0.1%) or in 
the knowle dge o f the z scale of the detector, as described 
in Section V A , could bias the Amd measurement be- 
cause these parameters are used to reconstruct the decay 
length difference Az and to convert it to the decay time 
difference At. The uncertainties on these quantities are 
propagated to Amd and lead to systematic uncertainties 
of 0.001 ps-i (1) and less than 0.002 ps'^ (j), respec- 



tively. In addition to these, we assign the difference of 
0.001 ps^^ (k) in the value of Amd obtained by using 
the Az to At conversion described in Eq. |2^ instead of 
Eq. |l| as a systematic error. Finally, in the likelihood 
fit, we fix the lifetime to the PDG value The 
present uncertainty on this value of ±0.032 ps leads to a 
systematic error of =f0.006 ps~^ (m). 
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TABLE XIII: Results from the likelihood fit to the At distributions of the hadronic B decays. The value for Am^ includes 
small corrections as described in the text. The first major column contains the fit results, while the second major column 
contains the correlation coefficients with respect to ArUd for each fit parameter. 



Parameter 


Fit Resuh 


Correlation 




Run 1 


Run 2 


Run 1 


Run 2 


Amd [ ps""^] 


0.516 ±0.016 




Signal Resolution Function 


Si (core) 


1.37 ±0.09 


1.18 ±0.11 


0.25 


0.16 


61 (At) lepton (core) 


0.06 ±0.13 


-0.04 ±0.16 


0.08 


0.00 


61 (At) kaon (core) 


-0.22 ±0.08 


-0.25 ±0.09 


0.03 


0.00 


61 (At) NTl (core) 


-0.07 ±0.15 


-0.45 ±0.21 


-0.00 


0.00 


fei(At) NT2 (core) 


-0.46 ±0.12 


-0.20 ±0.16 


0.01 


0.03 


62 (At) (tail) 


-5.0 ±4.2 


-7.5 ±2.4 


0.04 


0.06 


/2(tail) 


0.014 ±0.020 


0.015 ±0.010 


0.06 


0.07 


/a (outlier) 


0.008 ± 0.004 


0.000 ± 0.014 


-0.09 


0.01 


Signal dilutions 


(/?), lepton 


0.842 ± 0.028 


0.24 


\d), kaon 


0.669 ± 0.023 


0.30 


(D), NTl 


0.563 ± 0.044 


0.11 


(D>, NT2 


0.313 ±0.041 


0.11 


AD, lepton 


-0.006 ± 0.045 


0.02 


AD, kaon 


0.024 ± 0.033 


0.01 


AD, NTl 


-0.086 ±0.068 


0.00 


AD, NT2 


0.100 ±0.060 


0.00 


Background properties 


r, mixing bkgd [ps] 


0.853 ± 0.036 


-0.01 


/(r = 0), mixing bkgd, lepton 


0.05 ±0.10 


0.01 


/(t = 0), mixing bkgd, kaon 


0.42 ± 0.05 


0.01 


/(r = 0), mixing bkgd, NTl 


0.33 ±0.08 


0.01 


/(r = 0), mixing bkgd, NT2 


0.32 ±0.08 


0.01 


Background resolution function 


Si (core) 


1.211 ± 0.043 


1.131 ±0.046 


-0.00 


0.00 


61 (At) (core) 


-0.135 ±0.031 


-0.015 ±0.038 


-0.00 


-0.00 


/a (outlier) 


0.022 ± 0.004 


0.036 ± 0.007 


-0.01 


0.02 


Background dilutions 


(D), lepton, T — 


0.0 ±2.9 


-0.02 


(D>, kaon, r = 


0.52 ±0.08 


-0.03 


(D), NTl, r = 


0.67 ±0.27 


-0.01 


(D), NT2, r = 


-0.05 ±0.13 


-0.00 


(D), lepton, r > 


0.34 ±0.13 


0.02 


(D), kaon, r > 


0.26 ±0.06 


0.04 


(D), NTl, r > 


-0.13 ±0.11 


0.01 


(D), NT2, r > 


0.12 ±0.031 


0.01 



4-. Monte Carlo validation of measurement technique 

Candidate selection criteria, or the analysis and fit- 
ting procedure, could potentially cause systematic biases 
in the measurement of Am^. These potential biases are 
estimated by repeating the analysis with a large sam- 
ple of Monte Carlo events, which are generated with the 
full GEANT3 ^ detector simulation. In the Monte 
Carlo sample, the fitted result for Am^ is shifted by 
+0.007± 0.003 ps~^ from the input value. A correspond- 
ing correction with this central value is applied to the 
fitted result with data, and the uncertainty is assigned 
as a systematic error (n). 

The main cause of this bias is a small correlation be- 



tween the mistag rate and the At resolution that is not 
modeled in the likelihood function. This correlation is 
seen most readily in data for Kaon tags and is shown for 
simulation in Fig. p6^. We find that both the mistag rate 
for kaon tags and the event-by-event error cta* depend in- 
versely on v^pf , where pt is the transverse momentum 
with respect to the z axis of tracks from the -Stag de- 
cay. Correcting for this dependence of the mistag rate 
removes most of the correlation between the mistag rate 
and (TAt, as can be seen in Fig. [2^b. The mistag rate 
dependence originates from the kinematics of the physics 
sources for wrong-charge kaons. The three major sources 
of mistags are wrong-sign D'^ mesons from B decays to 
double charm, wrong-sign kaons from decays, and 
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TABLE XIV: Result of fitting for Amd in the entire Sflav 
sample and in various subsamples. The difference in the fitted 
value of Arrid versus the result from the fit to the full Bflav 
sample are reported. 



Sample 


Arud — Amdiall) 


Decay mode 


D*-Tx+ 


-0.029 ±0.030 


D*-p+ 


+0.017 ±0.039 


D*-a+ 


+0.066 ± 0.063 


D-7r+ 


+0.022 ± 0.030 


D-p+ 


-0.031 + 0.038 


D-a+ 


-0.033 + 0.041 


D*-X 


+0.000 + 0.025 


D-X 


-0.005 ± 0.023 


Tagging category 


Lepton 


+0.005 + 0.026 


Kaon 


+0.002 ± 0.023 


NTl 


-0.032 ± 0.044 


MT2 


+0.12 ± 0.10 


-Srec st ate 


Brcc = B" 


+0.015 + 0.023 


-Brer = B" 


-0.003 + 0.023 


Bu,,. state 


-Stag = -S" 


+0.019 + 0.023 


-Stag ~ -B" 


-0.007+ 0.022 


Data sample 


Run 1 


-0.012 + 0.022 


llmi 2 


+0.019 ^ 0.025 



TABLE XV: Systematic uncertainties and contributions to 
statistical errors for Am<j obtained with the likelihood fit to 
the hadronic B° sample. 



Source 


a{Amd) 


Signal properties 


(a) SVT alignment 

(b) At outlier description 


0.004 
0.002 


(c) Bcamspot position/size 


0.001 


(d) 0"^/ rotpiirt'iiic'iit 


().()()3 


Background properties 


(c) Background fraction 


0.002 


(f) Background At structure 


0.001 


(g) Background At resolution 

(h) Sideband extrapolation 


0.001 
0.002 


(i) Peaking B'^ background 


0.002 


External parameters 


(j) z scale 


<0.002 


(k) z boost (parameters) 
(1) z boost (method) 
(m) B° lifetime 


0.001 
0.001 
0.006 


Monte Carlo studies 


(n) Signal MC statistics 


0.003 


(o) Tag-side? D composition & lifetime 


0.001 


(p) Rigiit/wroug tag resolnticni diil'erenct^s 


0.001 


Total systematic error 


0.010 


Statistical error 


0.016 


Contribution due to resolution function 


0.005 


Contribution due to mistag rate 


0.005 


Total error 


0.019 



kaons produced directly in B decays. All these sources 
produce a spectrum of charged tracks that have smaller 
Pt than B decays that produce a correct tag. The 
At resolution dependence originates from the l/pf de- 
pendence of az for the individual contributing tracks. 

Since the effect is small and well described by the 
Monte Carlo simulation, we have chosen to treat the im- 
pact of this correlation as a correction, rather than build- 
ing the effect into the likelihood function. We include 
additional systematic errors related to the tag-side prop- 
erties that could affect the accuracy of the description of 
this correlation in the simtilation. In particular, the -D", 
and meson branching fractions, the D meson 
lifetimes, and the wrong-sign kaon production rates in B 
meson decays are all varied. These studies lead to an 
assigned systematic error of ±0.001 ps~^ (o). 

In addition, we consider the possibility that correctly 
and incorrectly tagged events could have different resolu- 
tion ftmctions. Based on Monte Carlo studies of the vari- 
ation in the fitted value for Arrid with and without allow- 
ing for independent resolution functions for correctly and 
incorrectly tagged events, an uncertainty of ±0.001 ps~"^ 
is assigned to this source (p). 



C. Validation studies and cross checks 

1. Monte Carlo studies 

A high-precision test of the fitting procedure was per- 
formed with fast parameterized Monte Carlo simulations, 
where 2000 experiments were generated with sample size 
and composition corresponding to that obtained from the 
actual data. The mistag rates and At distributions were 
generated according to the model used in the likelihood 
fit. The full fit was then performed on each of these ex- 
periments. The resulting distribution of pulls (defined as 
the difference between the fitted and generated value of 
a parameter divided by the statistical error as obtained 
from the likelihood fit) has a mean —0.038 ± 0.022 and 
standard deviation 1.012±0.023, consistent with no mea- 
surement bias in either the value of Am^ or its estimated 
error. 



2. Simple counting experiment 

If the mistag rate is known, the time-integrated frac- 
tion Xd = \^d/ i^-^x^d) mixed events can be determined 
from the iJflav sample by counting mixed and unmixed 
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events. The value for Xd obtained by this means, af- 
ter correcting for the mistag rates obtained from the full 
time-dependent fit and assuming the PDG value for the 
lifetime, leads to a value of Am^ — Xd/rga that dif- 
fers from the likelihood-fit result by —0.003 ±0.013 ps^^, 
where the quoted error is the difference in quadrature of 
the statistical errors of both measurements. 

Due to the choice for normalization of the likelihood 
Cmix, the time- integrated ratio of the number of mixed to 
unmixed events contributes to our measurement of Am^. 
Alternatively, it is possible to normalize the likelihoods 
of mixed and unmixed events individually, in which case 
Amd is determined solely from the shape of At distri- 
butions. The value of Arud determined by the At dis- 
tributions alone differs from the full measurement by 
0.003 ± 0.015 ps^^, where the quoted error is given by 
the difference in quadrature of the statistical errors of 
the two measurements. 



3. Cross check with Tgo 

If we allow the value of r^o to float in the Am^ fit 
the value of Am^ increases by 0.008 ± 0.007 ps^^, and 
the lifetime is found to be 1.51 ± 0.03 ps, consistent with 
our recent measurement |l^. We have also performed a 
series of fits with fixed values for Arud and Tgo in order 
to determine the dependence of Arud on r^o, which is 
found to be 



Arud 



0.516-0.279 



Tbo 



1.548 ps 



ps 



(42) 



VIII. CP VIOLATION IN NEUTRAL B DECAYS 
A. Likelihood fit results for sin2/9 

The value of sin2/3, the dilution factors Pi, the At 
resolution parameters Oi, and the background fractions 
and time distribution parameters are extracted with an 
unbinned maximum-likelihood fit to the flavor-eigenstate 
Bfiav and Bcp samples as described in Section We 
also demand a valid tag and At determination for the 
event, based on the algorithms described in Sections ^ 
and 0. The looser requirements |Ai| < 20 ps and (XAt < 
2.4 ps are applied to the proper time difference measure- 
ment. The fit results are summarized in Table XVI to- 
gether with the correlation of the parameters with sin2/3. 
The mistag fractions and vertex parameters are predomi- 
nantly determined by the -Bflav sample. The CP asymme- 
try and parameters describing the background for the CP 
events are determined by the CP sample. The value of 
sin2/5 obtained from the combined rjcp = — 1, rjcp = +1, 
and J/-0 K*^ CP samples is 

sin2/3 = 0.59 ±0.14 ±0.05, 



where the first error is statistical and the second system- 
atic. 

The rriEs distribution for events in rjcp = — 1 modes, 
separated into tagging categories, is shown in Fig. |2^. 
The signal probability ff^ig{mEs) for the rjcp = — 1 sam- 
ple is determined from these fits as described in Sec- 
tion 
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FIG. 27: Distribution of rriES for ricp ~ —1 candidates in sep- 
arate tagging categories (Lepton, Kaon, NTl and NT2), overlaid 
with the result of a fit with a Gaussian distribution for the 
signal and an ARGUS function for the background. 
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FIG. 28: a) Distribution of At for tagged events in the full CP sample. The upper (lower) panel is the sum of -B° (-B'')-tagged 
events in the r\cp = — 1 and J/xj} K*^ samples, combined with the (_B'')-tagged events in the J/ip K1 sample. Corresponding 
B°- (lower panel) and B^-tagged (upper panel) distributions for the b) ricp = —1, c) J/ip K°, and d) J/ip K*'^ samples are also 
shown. In all cases, the data points are overlaid with the result from the global unbinned likelihood fit, projected on the basis of 
the individual signal and background probabilities, and event-by-event At resolutions, for candidates in the respective samples. 
Therefore, the curves correspond to sin2/3 = 0.59, rather than the fitted value obtained with the individual subsample. The 
probability-weighted At spectra of the background candidates obtained from the fit are indicated by the shaded areas. 



Table XVII summarizes the event yields and sin2/3 



values determined for the full CP sample and various 
subsamples. Results are provided by CP channel, tag- 
ging category, versus tag, J/tp decay mode and 
data-taking period. The consistency between the six CP 
modes is satisfactory, the probability of finding a worse 
agreement being 8%. The large observed asymmetry in 
B'^ ^ XciKg causes the likelihood for this channel to 
become negative in certain regions of At. The likelihood 
of each of the selected candidates is of course positive. 
Fast parameterized Monte Carlo studies show that sin2/3 
is unbiased if the likelihood is not required to be positive 
for all values of At and that the probability to measure 



such a large asymmetry is about 1%. The observed asym- 
metry in the number of B'^ (160) and B'^ (113) tags in 
the J/^K'^ sample has no impact on the sin2/3 measure- 
ment. The results obtained with the full Bcp samples for 
Run 1 and Run 2 are consistent at the 1.8 sigma level. 
The y ields a nd fitted values for sin2/9 are also listed in 
Table KVII for the high purity r]cp = — 1 sample alone, 
along with a similar breakdown into subsamples; again, 
no significant variation is seen. 

The distribution of events as a function of At for B^ 
and tags is shown in Fig. p8| a for the full CP sam- 
ple. For this purpose, only those events with tties > 
5.27GeV/c^ in the rjcp = —1 and J/ipK*° samples or 
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Ai? < lOMeV in the rjcp = +1 sample are included. 
Overlaid on the data are the projections of the signal 
and background At distributions obtained from the fit, 
where the latter is normalized to the projected back- 
ground level. Figure p8[b-d shows the corresponding At 
distributions for the rjcp = — 1, rjcp = +1 samples and 
5° ^ J/%PK*'^ {K*° K°TT°, K° ^ n+TT-). The super- 
imposed likelihood curves show the quality of the fit for 
each subsample. The value of sin2/3 obtained by fixing 
all other parameters to results obtained with the full CP 
sample and then fitting for sin2/3 in bins of At is shown 
in Fig. p9|a. The values obtained for sin2/3 are all consis- 
tent, demonstrating that the oscillation as a function of 
At has the expected behavior. The observed asymmetry 
A{At) is shown in Fig. p9| b and c for the rjcp = ^1 and 
ricp = +1 samples respectively, along with the projec- 
tions from the fit results. 

The average dilutions and dilution differences for B'^ 
and tags obtained from the fit to the S'^ flavor eigen- 
state and full CP sample, and the corresponding tagging 
efficiencies, are summarized in Table XVII] . We find a 



CO. 



1.5 



total tagging efhcicncy of (68.4 ±0.7)% (statistical error 
only). The lepton categories have the lowest mistag frac- 
tions, but also have low efficiency. The Kaon category, 
despite having a larger mistag fraction (17.6%), has a 
higher effective tagging efficiency; one-third of events are 
assigned to this category. Altogether, lepton and kaon 
categories have an effective tagging efficiency Q « 22.4%. 
The neural network categories increase the effective tag- 
ging efficiency by ~ 4% to an overall Q — (26.1 ± 1.2)% 
(statistical error only). These mistag fractions are very 
similar to the mistag fraction that are obtained from the 
Bflav sample alone (see Table XIII ). The small differences 
are due to the correlation between the mistag fractions 
and the At resolution function parameters. 

Based on a large number of fast parameterized Monte 
Carlo experiments with the same number of events as our 
full Bcp and Bflav data samples, we estimate a probabil- 
ity of 27% for finding a value of the maximum likelihood 
lower than that observed. These same studies, based on 
samples with the size and composition of the data, show 
that the expected statistical error is 0.132 with a spread 
of 0.005, in very good agreement with the observed error 
of 0.137. 



B. Systematic error estimation 

Just as for the Arrid measurement, systematic errors 
can usefully be grouped into signal description, including 
detector reconstruction effects, background description, 
fixed external parameters, and statistical limitations of 
Monte Carlo validatio n tests for the fitting procedure 
(discussed in Section VIII CI). A summary of th ese 
sources of systematic error is shown in Table XIX for 
the various CP samples. In the following, the individual 
contributions are referenced by the lettered lines in this 
table. 



0.5 



-5 



At (ps) 




At (ps) 



FIG. 29: a) Fitted value of sin2/3 obtained in bins of At by 
fixing all other parameters to the values obtained with the 
full CP sample; Raw asymmetry in the number of B" and 
tags in the signal region, {Ngo — Ngo)/iNgo + N^o), with 
asymmetric binomial errors, as a function of At for b) rjcp ~ 
— 1 and c) J/i(; K1 samples. The data points are overlaid with 
the separate fit results for the two samples. 



1. Signal properties and description 

The parameters of the At resolution function, the di- 
lutions and dilution differences are determined from the 
data sample itself with the likelihood fit. Thus, they do 
not contribute to the systematic error, but rather are 
incorporated into the statistical uncertainty at a level 
determined by the size of the data sample itself. Their 
overall contribution to the total error on sin2/? is 0.02, 
as determined from the difference in quadrature between 
the statistical error on sin2/3 from the full likelihood fit 
and from a fit with only sin2/3 allowed to vary. 

While the bulk of the uncertainties from these sources 
is thus incorporated into the statistical error, we assign 
additional systematic uncertainties due to the fixed form 
of the parameterization for the At resolution function. 
This form may not be flexible enough to account for all 
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TABLE XVI: Parameters for the combined likelihood fit to the Bcp and Bflav samples. The first major column contains the 
fit results, while the second major column contains the correlation coefficients with respect to sin2/3 for each fit parameter. 



Parameter 


Fit Result 


Correlation 




Run 1 


Run 2 


Run 1 


Run 2 


sin2/3 


0.59 ±0.14 




Signal Resolution Function 


Si (core) 


1.2 ±0.1 


1.1±0.1 


0.018 


0.020 


bi{At) lepton (core) 


0.07 ±0.12 


0.04 ±0.16 


0.008 


0.045 


bi{At) kaon (core) 


-0.26 ±0.08 


-0.18 ±0.09 


0.002 


0.021 


bi{At) NTl (core) 


-0.21 ±0.15 


-0.33 ±0.21 


0.004 


0.001 


bi{At) NT2 (core) 


-0.31 ±0.11 


-0.17 ±0.15 


-0.001 


-0.002 


62 (At) (tail) 


-1.7± 1.5 


-3.3 ±2.8 


0.001 


0.006 


/2(tail) 


0.08 ± 0.06 


0.04 ± 0.04 


0.009 


0.005 


/3 (outlier) 


0.005 ± 0.003 


0.000 ± 0.001 


-0.001 


0.000 


Signal dilutions 


(D), lepton 


0.82 ±0.03 


-0.042 


(D), kaon 


0.65 ±0.02 


-0.083 


(D), NTl 


0.56 ±0.04 


-0.015 


(D), NT2 


0.30 ±0.04 


-0.032 


AD, lepton 


-0.02 ±0.04 


0.010 


AD, kaon 


0.04 ±0.03 


0.005 


AD, NTl 


-0.11 ±0.06 


0.014 


AD, NT2 


0.12 ±0.05 


-0.008 


Background properties 


r, mixing bkgd [ps] 


1.3 ±0.1 


-0.001 


fir = 0), CP bkgd 


0.60 ±0.12 


-0.011 


/(r — 0), mixing bkgd, lepton 


0.31 ±0.10 


-0.001 


/(r — 0), mixing bkgd, kaon 


0.65 ±0.04 


-0.001 


/(r = 0), mixing bkgd, NTl 


0.62 ± 0.06 


-0.001 


/(r = 0), mixing bkgd, NT2 


0.64 ±0.04 


-0.001 


Background resolution function 


Si (core) 


1.5 ±0.1 


1.3 ±0.1 


0.004 


-0.003 


61 (At) core [ps] 


-0.16 ±0.03 


0.02 ± 0.04 


0.000 


-0.001 


/2 (outlier) 


0.016 ± 0.004 


0.017 ±0.005 


-0.001 


0.000 


Background dilutions 


(D), lepton, r = 


0.33 ±0.27 


0.003 


(D), kaon, r = 


0.45 ± 0.03 


0.008 


(D), NTl, r = 


0.25 ±0.10 


0.002 


(D), NT2, r = 


0.11 ±0.06 


0.003 


(D), lepton, r > 


0.33 ±0.14 


0.000 


(D), kaon, r > 


0.24 ±0.06 


0.000 


(D), NTl, r > 


0.05 ±0.14 


-0.001 


(D), NT2, r > 


0.09 ± 0.09 


0.000 



possible effects. In addition, tests of the assumption that 
the resolution function and dilution parameters are the 
same for the Bflav and Bcp samples are limited in preci- 
sion by the size of the available Monte Carlo samples. 



The resolution function, described in Section |V H , is 
one of several possible functional forms. In order to test 
possible biases induced by this particular choice, an al- 
ternative model has been considered where a Gaussian 
distribution is convolved with an exponential, with the 
effective lifetime in the exponential depending on the tag- 
ging category. No difference between the fit results with 
the two models is observed in Monte Carlo simulation. 
We assign as a systematic uncertainty the difference in 
the fit results observed in the data (Table pCIX , line a). 



The largest systematic uncertainties from the At behav- 
ior arises from possible effects that our model of the res- 
olution function cannot accommodate or completely pa- 
rameterize. These include residual uncertainties in the 
SVT alignment (b) and possible differences in the At de- 
termination for correctly and incorrectly tagged events 
(c). An additional uncertainty is assigned due to the 
treatment of the A< outliers (d). Fits with Monte Carlo 
samples of -Bflav and Bcp signal events show no signif- 
icant difference between resolution function parameters 
for the two samples. We assign a systematic uncertainty 
of ±0.003 due to the residual shift in sin2/3 between the 
two sets of fitted At resolution parameters (e) . 

An underlying assumption of the global fit is that di- 
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TABLE XVII: Result of fitting for CP asymmetries in the entire CP sample and in various subsamples. The yields are obtained 
with likelihood fits and are therefore background subtracted. 



Sample 


Nta.g 


Purity (%) 


sin2/3 


Im\/\\\ 


1^1 


CP sample 




oU 


n piQ-i-n 1 A 
u.oyztu. ±41 






J/tpK g (iVg — > TT TT ) 


316 


98 


0.45 ± U.18 


0.45 ± 0.18 


n n 1 _L All 

0.91 ± 0.11 


7 1,1, rro (Tro , ^o_o\ 

■J/W l^S K'^S ^ TT TT j 


fiA 




U. ( U zt u.ou 


u. / 1 =1= u.ou 


n ctK 4- n 07 


ll}{Zb)Ks (As ^ -K -K ) 


D ( 


HQ 

ao 


U.4 / dz U.4z 


U.48 ± U.4o 


i.zz zh U.OO 


XclAg ^Ag TV TV ) 


00 


Q7 


9 CQ4-0.55 


9 fi7 + n "iQ 
z.u / 31 u.uy 


n 71 +0 9"^ 




Z / 


or 


n 70 4- n "^/i 
u. ( u m u.o'i 










74- 


U.OZi HZ _L.UU 






Lepton 


ioU 


QO 

oz 


U.04 dz U.zD 






Kaon 


438 


79 


0.58 zb 0.18 






NTl 


79 


74 


0.89 zb 0.30 






NT2 


156 


80 


0.40 zb 0.65 








420 


79 


0.54 zb 0.19 








383 


78 


64 zb 20 






jyT/j — e^^e" 


385 


78 


0.49 ± 0.20 








418 


84 


70 ± 18 






Run 1 


533 


80 


0.49 ± 0.20 






Run 2 


270 


84 


0.82 ± 0.22 






T^cp ~ — 1 sample 




yo 


U.OD It U. iO 


U.OO It U.iO 


u.yo It u.uy 


Lepton 


/ 4 


iUU 


u.o4 it u.^y 


n r7 _i_ n 9ri 
U.O / ± u.zy 


A 77 _l_ A 1 /I 
U.I 1 ± U.i4 


Kaon 


271 


98 


0.59 ± 0.20 


0.59 ± 0.20 


0.98 ± 0.12 


NTl 


46 


97 


0.67 ±0.45 


0.57 ±0.46 


0.73 ± 0.29 


NT2 


89 


95 


0.10 ±0.74 


0.28 ± 1.29 


2.95 ±3.83 




234 


98 


0.50 ±0.22 








246 


97 


0.61 ±0.22 






J/'ip — > e+e~ 


219 


94 


0.54 ±0.22 


0.52 ±0.22 


1.00 ±0.15 




261 


98 


0.60 ±0.21 


0.63 ±0.21 


0.87 ±0.11 


Run 1 


310 


95 


0.37 ±0.20 


0.37 ±0.20 


1.16 ±0.15 


Run 2 


170 


98 


0.86 ±0.24 


0.96 ±0.26 


o.66;;!:l^ 


Control samples 




7579 


84 


0.00 ± 0.04 








6800 


86 


-0.02 ±0.04 








705 


95 


0.12 ±0.12 






B+ ^ J/ipK^''^+, ip{2S)K+ 


2031 


94 


0.07 ±0.07 







TABLE XVIII: Average mistag fractions Wi and mistag differences Awi — Wi{B'^) — Wi{B^) extracted for each tagging category 
i from the maximum-likelihood fit to the time distribution for the fully-reconstructed sample (-Bflav+Bcp). The figure of 
merit for tagging is the effective tagging efficiency Qi — e,;(l — 2wi)^, where £i is the fraction of events with a reconstructed 
tag vertex that are assigned to the i"* category. Si is computed for the ijcp = ±1 samples as well as the combined Bcp and 
Bflav samples. Uncertainties are statistical only. The statistical error on sin2/3 is proportional to 1/^/Q, where Q ~ "^2 Qi- 



Category 


rtcp = —1 


Vcp = ±1 


Bflav ± Bcp 


e[%\ 


e[%\ 


s[%\ 


w [%\ 


Aw [%\ 


Q [%\ 


Lepton 
Kaon 
NTl 
NT2 


11.0± 1.2 
38.9 ± 1.9 
6.9 ±0.9 
13.0 ±0.4 


10.4 ±3.0 
28.3 ±4.5 
4.8 ±2.3 
13.9 ±3.3 


10.9 ±0.3 
35.8 ± 1.0 
7.8 ±0.3 
13.8 ±0.3 


9.0 ± 1.4 
17.6 ± 1.0 

22.0 ±2.1 

35.1 ± 1.9 


0.9 ±2.2 
-1.9 ± 1.5 

5.6 ±3.2 
-5.9 ±2.7 


7.4 ±0.5 
15.0 ±0.9 

2.5 ±0.4 
1.2 ±0.3 


All 


69.8 ±2.7 


57.4 ±6.7 


68.4 ±0.7 






26.1 ± 1.2 



lutions and dilution differences are the same for the i?flav 
and Bcp samples. We assign the full difference as seen 
in Monte Carlo simulation as systematic error, ±0.027 
(f). In addition, the data sample was used to study 



any possible dependence of the dilutions on At. No sig- 
nificant effect was observed. However, a dependence of 
the dilutions on (Ta* has been seen, b oth in d ata and the 
Monte Carlo simulation (see Section VII B 4). Finally, it 
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TABLE XIX: Summary of contributions to the systematic error on sin2/3, /mA/|A| and |A|. Note that the last two measurements 
use only the -qcp = — 1 sample. 



Source 


CP Sample 


ncp ^ -l\J/ipKl\J/'tPK'"\ Full |/mA/|A|| |A| 


Signal parameters 


(a) At signal resolution model 


±0.009 1 ±0.01 1 ±0.07 1 ±0.009 1 ±0.003 


±0.003 


(b) SVT aligmnent 


±0.027 


±0.012 


(c) At for right /wrong tagged events 


±0.012 


±0.011 


±0.003 


(d) At signal resolution outliers 


±0.002 1 ±0.018 1 ±0.03 1 ±0.002 


±0.003 


±0.002 


(e) At signal resolution 


±0.003 


±0.009 


(f) Signal dilutions for CP vs. Bflav 


±0.027 


±0.011 


(g) Tagging Efficiencies 


±0.003 1 ±0.004 


±0.012 


Background properties: rjcp = — 1 


(h) Background fraction 


±0.006 






±0.005 


±0.006 


±0.004 


(i) CP bkgd peaking component 


±0.004 






±0.003 


±0.005 


±0.001 


(j) CP bkgd CP content (ARGUS) 


±0.015 






±0.015 


±0.015 


±0.001 


(k) CP bkgd CP content (Peak) 


±0.004 






±0.004 


±0.004 


±0.001 


(1) CP bkgd effective lifetime 


















(m) CP bkgd resolution 


±0.002 






±0.002 


±0.002 


±0.001 


Background properties: J/ipK'^ 


(n) Background fraction 




±0.075 




±0.01 






(o) AE distribution 




±0.04 




±0.007 






(p) Effective CP of backgrounds 




±0.020 




±0.001 






(q) Background composition 




±0.014 




±0.002 






(r) Background At and dilution 




±0.023 




±0.003 






Background properties: J/^jj K*" 


(s) Sample composition 






±0.08 


±0.001 






(t) Ri_ 






±0.08 


±0.001 






Background properties: Bflav 


(u) Background fraction 


±0.001 


±0.008 


±0.003 


±0.002 


±0.002 


±0.001 


(v) -Bflav bkgd mixing contrib. 


±0.001 


±0.002 


±0.001 


±0.002 


±0.001 





(w) Bflav bkgd peaking component 





±0.001 


±0.001 











External parameters 


(x) z scale and boost 


±0.003 


±0.001 


(y) Beam spot 


±0.002 


±0.006 


(z) B" lifetime 


±0.008 


±0.011 


±0.022 


±0.009 


±0.009 


±0.012 


(aa) Anid 


±0.015 


±0.012 


±0.082 


±0.013 


±0.015 


±0.001 


Monte Carlo studies 


(bb) Monte Carlo statistics 


±0.012 


±0.007 


Total systematic error 


±0.05 


±0.10 


±0.16 


±0.05 


±0.05 


±0.02 


Statistical error 


±0.15 


±0.34 


±1.01 


±0.14 


±0.15 


±0.09 



is possible that tagging efficiencies could be different for 
and 5° mesons. A separate study o f the re lative tag- 
ging efficiencies is described in Section VIII D , since the 
relative efficiencies form an important part of the direct 
CP violation search. The systematic error on sin2/3 due 
to this effect is estimated to be ±0.003 (g). 



2. Background properties 

The fraction of background events in the rjcp = — 1 
sample is estimated from fits to the rngg distribution. 
Varying this fraction within the stated errors and chang- 
ing the signal probability as a function of rngg results in 
a systematic error of ±0.005 on sin2/3 (h). The uncer- 



tainty on the fraction of peaking background contributes 
a systematic error of ±0.003 (i). Varying t he effe ctive 
sin2/3 assumed for the A RGUS [ A in Section |VIB l| ) and 
peaking ((5pcak in Section VI B 1 ) backgrounds in the CP 
sample from —1 to +1 contributes a systematic error of 
±0.015 (j) and ±0.004 (k), respectively. In addition, the 
contributions due to the uncertainty of the At resolution 
model (±0.002), and the effective lifetime (negligible) of 
the CP background, have been evaluated (l-m). 

For the ~* J/ipK^ channel, the signal and non- J/ip 
background fractions are varied within their statistical 
uncertainties (ilc) as obtained with the fit to the AE 
distribution of the sample. This contributes a system- 
atic error of ±0.075 to the B° J/^K^ sin2/3 result 
and ±0.01 to the final result (n). We also vary back- 
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ground parameters for the B° J/tpK^ sample, includ- 
ing the J/tp X branching fractions according to Table [v| , 
the assumed ijcp, the mistag rates and efficiencies, the 
At resolution function, and AE shape (o-r). The to- 
tal — > J/ij^K^ background systematic error, summing 
these contributions in quadrature, is ±0.09 for the sin2/3 
fit to the — > J/ipK^ sample alone and ±0.013 for the 
full sample. 

For the B° J/^pK*" {K'^n") sample, the value of 
R± as well as the sample composition are varied (s-t) 
according to Table |^. 

The effect of the uncertainty on background compo- 
nent in the Bflav sample on sm2f3 has also been eval- 
uated. The only significant sources of uncertainty are 
the fraction of background that mixes (v) and the signal 
probability distribution as a function of toes (u,w). 



3. External parameters 

The residual uncertainty on the physical z scale (x) 
and the boost parameters of the T(4S') center of mass (y) 
contribute systematic uncertainties. We fix the i?" life- 
time to the current world average values tbo = 1.548 ps 
and Arrid = 0.472 ps~^ ||l^. The errors on sin2/3 due to 
uncertainties in tb" and Auid are ±0.009 and ±0.013, 
respectively (z-aa). 



4-. Monte Carlo validation of measurement technique 

The analysis method has been studied with a high- 
statistics Monte Carlo sample. A fit result that is consis- 
tent with the generated value for sin2/3 was found. We 
assign a ±0.012 systematic error due to the statistical 
limit ation of the Monte Carlo sample size (bb). Sec- 
tion VIII C 1 describes this study in more detail. 



C. Validation studies and cross checks 

We have used data and Monte Carlo samples to per- 
form validation studies of the analysis technique. These 
tests include studies with parameterized Monte Carlo 
samples, full GEANT3 |^ simulation samples, as well 
as data samples where no CP asymmetry is expected. 



1. Monte Carlo studies 

The highest precision test of the fitting procedure was 
performed with fast parameterized Monte Carlo simula- 
tion, where 1000 experiments were generated with sample 
sizes corresponding to the observed -Bflav and Bcp events 
in data, including mistag rates, At resolutions, and back- 
ground fractions and time dependence. The full fit is per- 
formed on each of these experiments. The resulting pull 
distribution (defined as the difference between the fitted 



and generated value of a parameter divided by the statis- 
tical error as obtained from the likelihood fit) has a mean 
-0.029±0.032 and standard deviation 1.007±0.022, con- 
sistent with no measurement bias in either the value of 
sin2/3 or its estimated error. 

In addition, large samples of signal and background 
Monte Carlo events generated with a GEANT3 ||2^ de- 
tector simulation are used to validate the measurement. 
For these tests, we obtained the resolution function pa- 
rameters as well as the dilutions from a Monte Carlo 
sample of -Bflav events. Using these parameters, we fit for 
sin2/3 in Monte Carlo samples of CP signal events that 
correspond in number to the reconstructed data sample. 
These Monte Carlo events are generated with various 
values of sin2/3 (0.1 to 0.9) and different CP-eigenstate 
modes, corresponding to those used in the measurement 
with data. The mean and spread of the pull distribution 
for these Monte Carlo samples can be used to check for 
any measurement bias and to confirm the validity of the 
reported error. We find that the mean pull is consistent 
with zero and the spread is consistent with the reported 
error. A systematic error of ±0.012 is assigned to sin2/3 
due to the limited Monte Carlo statistics for this test. 

The effect of background has been evaluated by adding 
an appropriate fraction of background events to our sig- 
nal Monte Carlo sample and performing the likelihood 
fit. The background samples are obtained either from 
simulated B J/ipX events or AE' sidebands in data 
(|Ai?| < 120 MeV but outside the signal region). We find 
no significant bias for sin2/3 with the addition of either 
source of background. 



2. Cross checks with Tgo and Arrid 



Table XX shows results for sin2/3 if Ato^ and Tbo are 
allowed to float in the combined fit to the CP and Bflav 
samples. The fitted value of Aro^ is somewhat larger 
than that reported in Section VII . However, with no 
kaon veto applied to the tagging vertex, the correction 
for the bias introduced by known correlations between 
mistag rates and the At resolution is also larger. Tak- 
ing this into account, the two results are consistent within 
the independent statistical errors. Likewise the lifetime is 
found to be consistent with our recent measurement [ pT[ . 
We have also performed fits with Ato^ and r^o fixed to 
a series of values around the world average in order to 
determine the dependence of sin2/3 on these two param- 
eters, thereby finding that 



sin2/3 = [0.59 - 0.35 
- 0.45 



Arrid 



0.472 ps-i 

tb" _ 
1.548 ps 



(43) 
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TABLE XX: 
the sin2/3 fit 
sample alone. 



Results when Am^ and (or) Tgo are floated in 
to the full CP sample and the r^cp = ^1 sub- 



Fit 


sin2/3 


Amd (ps '-) 


Tgo (ps) 


All CP modes 


Nominal fit 


0.59 


± 


0.14 


0.472 


1.548 


Float Amd 


0.55 


± 


0.13 


0.533 ±0.015 


1.548 


Float Tgo 


0.60 


± 


0.14 


0.472 


1.53 ±0.03 


Float Amd and Tgo 


0.56 


± 


0.13 


0.542 ±0.016 


1.50 ±0.03 




Vcp 




— 1 modes 




Nominal fit 


0.56 


± 


0.15 


0.472 


1.548 


Float Arrid 


0.51 


± 


0.15 


0.531 ±0.015 


1.548 


Float TgO 


0.57 


± 


0.15 


0.472 


1.53 ±0.03 


Float Arud and Tgo 


0.52 


± 


0.15 


0.540 ±0.016 


1.50 ±0.03 



3. Asymmetries in data control samples 

Control samples in data where the reconstructed 
and meson decays to a flavor-eigenstate mode with a 
D^*'> or charmonium meson in the final state can be used 
to validate the sin2/3 measurement, since the asymmetry 
is expected to be zero in this case. For these samples, the 
At resolution function parameters and the dilutions are 
fixed to the values obtained with the i?fiav sample. The 
CP asymmetry and the fraction of prompt background 
(identical for each tagging category, as is the case for 
the fit to the CP data sample) are allowed to float. The 
measured asymmetries are all consistent with zero, as 
shown in Table XVIl. The observed At distributions 
for the and i3"-taggcd events in the Bflav sample 
is shown in Fi g. ^0^ , where good agreement is clearly 
visible. Figure pQb demonstrates that there is no visible 
asymmetry as a function of At. 

Control samples are also used to check the assump- 
tion that the At resolution function, which is primarily 
determined by the Bflav sample, can be applied to the 
charmonium decay modes in the CP sample. Figure ^ 
graphically compares the fitted At resolution function 
for the D'-*^X control sample with that of the 

_B+ ccX control sample. A la error envelope encom- 
passes the fit to the ccX sample, which has five 
times fewer events. The level of agreement is accept- 
able. The same comparison between the B° D'-*^X 
and B^ —f ccX samples was inconclusive due to the low 
statistics of the i?° ccX sample. 



4- Time- integrated measurement of mistag rates 

As described in Section I, a time-integrated technique 
can also be used to measure the mistag fractions in data, 
thereby providing a simple check of the likelihood fit 
method. The statistical precision of the time-integrated 
measurement is enhanced by restricting the sample to 
events in a single optimized At interval. Taking into ac- 
count detector vertex resolution, the optimal interval is 




At (ps) 



FIG. 30: a) Distribution in At for B°- and B°-tagged and 
b) observed asymmetry for events in the flavor-eigenstate B'^ 
sample. The projections of the likelihood fit for the S"- and 
-B^-tagged samples are shown in a) as the overlapping solid 
lines. 




At - At._..„ (ps) 



FIG. 31: Comparison of the fitted At resolution obtained 
with the data control samples — > D'-'-'X and B^ ccX, 
showing the fitted distribution for a) 5t = At — Attruo and 
b) the normalized difference 5t/cr/s.t. The one sigma error 
envelope from the fit to the B^ — > ccX sample (shaded 
region), overlaps the central value for the five-times larger 
B^ D'*'X sample (dashed line). 



found to be \At\ < 2.5 ps. Events with \At\ > 2.5 ps 
have, on average, equal numbers of mixed and unmixed 
events due to flavor oscillations, and therefore contribute 
nothing to the determination of the mistag rate. We refer 
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to this time-integrated technique using a single optimized 
Ai interval as the "single-bin" method and apply it to 
both the -Bflav sample described in Section [II C 1 and 
the semileptonc sample described in Section IIIE| . 

To correct for the presence of backgrounds, a term is 
added to Eq. |l| to account for the contribution of each 
background source to the fraction of mixed events in the 
sample: 



Xobs = fsigiXd + (1 



^Xd)w)+J2ffJXi3, (44) 





where /gig and are the fraction of signal and back- 
ground source f3, respectively, xp is the fraction of mixed 
events in each background source, and Xobs is the ob- 
served fraction of mixed events. In this expression, Xd 
must also be modified to represent the integrated mixing 
probability for |At| < 2.5 ps. Using the world-average 
values for ArUd and Tgo |ll| , and taking into account the 
At resolution function TZ(At), we find 



Xd = ^[1 



= -|At|/r 



cos(AmdAi) «> 7^(A^) d{At)/ 



-|At|/r , 



> 7^(A^) d(Ai)] = 0.079, 



(45) 



where the integral is performed over the range \At\ < 
2.5 ps and TZ{At) is modeled by a double-Gaussian dis- 
tribution with five parameters (one fraction, two biases 
and two widths) determined directly from data using the 
hadronic sample. Solving Eq. ^ for w, and using the 
calculated value for x'di the mistag rates are obtained: 



Xobs - fsigX'd ~ I]/3 fpXP 
/sig(l - ^X'd) 



(46) 



TABLE XXI: Yields, efficiencies, mistag rates w, and tagging 
separation Q — etag(l — 2w)^ as measured by the single-bin 
method in the hadronic Bflav event sample. A comparison of 
the mistag rates measured in the same sample with th e single - 
bin method, Wsh, and the likelihood fit, uiiikc (Table [KVIII ), 
are reported as the differences Ankc = ''^'sb — wuko between 
the two extraction techniques, normalized to the uncorrelated 
statistical and systematic errors. 



Category 


Yield 


Efficiency 


Mistag rate w 


Q 


Alike 






[%] 


[%] 


[%] 




Lept on 


1128 


11.0 ±0.3 


9.5 ±1.5 ±0.6 


7.2 


-0.8 


Kaon 


3687 


35.8 ±0.5 


17.8 ± 1.0 ±0.7 


14.8 


-0.4 


NTl 


819 


7.9 ±0.3 


22.0 ±2.2 ±0.9 


2.5 


±0.0 


NT2 


1428 


13.9 ±0.3 


34.3 ± 1.9 ± 1.1 


1.4 


±0.8 



All tagged events in the Sflav sample with \At\ < 
2.5 ps are used for a single-bin study. The combinato- 
rial background fraction in the signal sample is deter- 
mined from a fit to the rriEs distribution as described 
in Section VII. The signal region is defined as events 



with Toes > 5.27GeV/c^. The peaking background 
in this signal is estimated to be (1.3 ± 0.8)%. The frac- 
tion of mixed events in the combinatorial background is 
determined by tagging category with the sideband con- 
trol sample, 5.20 < toes < 5.27GeV/c^, and the mistag 
fraction associated with the B'^ peaking background has 
been measured directly in data. The number of tag ged 
events in each category is summarized in Table XXI. 



TABLE XXII: Yields, efficiencies, mistag rate w, and tagging 
separation Q — (-tagi^ — 2wY as measured by the single-bin 
method in the semileptonic event sample. A compari- 
son of the mistag rates measured with the single-bin method 
are reported as the differences Agampic = Wflav — be tween 
the mistag rates in the -Bflav sample, ?j;flav (Table XXI), and 
semileptonic samples, Wsi, normalized to the quadratic 
sum of statistical and uncorrelated systematic errors. 



Category 


Yield 


Efficiency 


Mistag rate w 


Q 


^sample 






[%] 


[%] 


Wo] 


\o] 


Lept on 


3046 


11.9 ±0.4 


8.7 ± 0.9 ± 1.4 


8.1 


±0.4 


Kaon 


10270 


36.2 ± 1.9 


19.5 ± 0.7 ± 1.2 


13.5 


-1.1 


NTl 


2127 


8.1 ±0.4 


22.3 ± 1.4 ± 1.2 


2.5 


-0.1 


NT2 


3967 


13.5 ±0.9 


36.0 ± 1.2 ± 1.3 


1.1 


-0.7 



A separate single-bin analysis is also performed with 
the sample of B^ D*~t^V events described in Sec- 
We use tagged events with \At\ 



tion 



HIE 



< 2.5 ps 

and evaluate the backgrounds for events in this time in- 
terval. The backgrounds and mixed-event fractions are 
evaluated separately for each tagging category. Back- 
grounds are larger for the semileptonic modes than for 
the hadronic modes and originate from a variety of 
sources. In the case of the combinatorial background, 
the estimate is obtained from the m{D^-K~) — m{D^) 
sideband. For the continuum background, off-resonance 
data is used after correction for the combinatorial com- 
ponent. The mixed-event fraction for BB background 
is estimated with generic BB Monte Carlo simulation. 
The mistag fraction of the last background component, 
the decay B~^ D*XiT7, has been determined with data. 
The estimates of the contributions of the various back- 
grounds are described in Section HI E . The number of 
tagged events in each category are summarized in Ta- 



ble XXIl 



We use Eq. ^ to obtain the mistag rates in each tag- 
ging categ ory sh own in Table XXI for the Bfiav sample 
and Table XXII for the semileptonic sample. The 



sources of syst ematic erro r on th ese results are summa- 
rized in Tables XXIII and XXIV respectively. 



Three sources of systematic uncertainties are common 
to both the hadronic and semileptonic samples. The first 
is the uncertainty due to the errors on the world-average 
values for the B° lifetime and Ato^ values. The second is 
due to the At resolution function, whose fit parameters in 
data are varied within errors. The third common uncer- 
tainty is related to the possibility that wrong tags have 
worse At resolution than correct tags. This effect has 
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TABLE XXIII: Sources of systematic error for the mistag 
measurement on the Bflav sample in the single-bin method. 



Type 


Variation 


Lepton 


Kaon 


NTl 


NT2 


r(B"), Amd 


±lcr 


0.005 


0.004 


0.003 


0.002 


Resolution 


see text 


0.002 


0.002 


0.001 


0.001 


Wrong-tag resolution 


see text 


0.003 


0.006 


0.007 


0.009 


Combinatorial bkgd 


±la 


0.002 


0.002 


0.005 


0.004 


peaking bkgd 


±1(7 


0.001 


0.001 


0.000 


0.000 


Total 




0.006 


0.007 


0.009 


0.011 



TABLE XXIV: Sources of systematic error for the mistag 
measurement from the semileptonic B'^ sample in the single- 
bin method. 



Type 


Variation 


Lepton 


Kaon 


NTl 


NT2 


r(B"), Amd 


±1(T 


0.006 


0.004 


0.004 


0.002 


Resolution 


see text 


0.001 


0.001 


0.001 


0.001 


Wrong-tag resolution 


see text 


0.003 


0.006 


0.007 


0.009 


Combinatorial bkgd 


±1(7 


0.001 


0.006 


0.004 


0.004 


Contiimum bkgd 


±1(7 


0.001 


0.003 


0.006 


0.007 


BB bkgd 


±1(7 


0.011 


0.005 


0.005 


0.004 


B+ -> D*Xeu bkgd 


±1(7 


0.003 


0.004 


0.002 


0.001 


Total 




0.014 


0.012 


0.012 


0.013 



been studied with Monte Carlo simulation, where we ob- 
serve a slightly larger RMS width for events with wrong- 
sign tags. From this study, scale factors comparing the 
right and wrong-tag resolution functions have been ex- 
tracted and then applied to the resolution function for 
wrong tags. 

The systematic uncertainties unique to each sample 
are due to the background components. These are esti- 
mated by varying both the background fractions f/^ and 
the fraction of mixed events associated with each back- 
ground source, x/3) by one standard deviation in their 
uncertainty. 

For the semileptonic sample, the systematic error due 
to backgrounds is the dominant source. The charac- 
terization of these various backgrounds is described in 
Section |III E|. For the combinatorial background frac- 



tion, a relative systematic uncertainty of 20% is added in 
quadrature to the statistical error to cover the range of 
results obtained with various m(D°7r~) — m{D°) fitting 
functions. 

The systematic error due to the continuum background 
is determined by varying both the background level and 
the mixed fractions. The BB background fraction un- 
certainty is obtained by combining the statistical uncer- 
tainty and the systematic error given in Section III E . 
The systematic errors introduced by uncertainties on 
the background from the decay B^ D*X£T' are ob- 



performed to verify that the mistag fractions are not af- 
fected by the presence of the extra pions in the decay 
_B+ D*X£P. An additional uncertainty due to the 
statistical precision of the Monte Carlo study has been 
added to the charged B mistag fractions measured with 
data. 

Table XXI shows the difference Auke = ifsb — ifiike 
between the mistag rates measured with the single-bin 
method in the -Bfiav sampl e, Wsh, and the likelihood fit 
result, wiike (Table XVIII ). The difference is reported 
in terms of the uncorrelated statistical and systematic 
errors for the two methods, when applied to the same 
data sample. The component of uncorrelated statistical 
error is estimated with a fast parameterized Monte Carlo 
simulation. It varies with category due to different event 
yields. The differences Agampic = wnnv — Wsi in the mistag 
rates measured with the single-bin method in the -Bfiav 
sample, Wfiav, and in the se mileptonic B'^ sample, Wsi, 

The quadratic sum of the 



are reported in Table XXII 



tained by varying the fraction described in Section III E 
as well as the mistag fraction of B~^ mesons measured 
on data. Studies with Monte Carlo simulation have been 



statistical and uncorrelated systematic errors is used to 
estimate the consistency of the measurements. 



5. Vertexing algorithm checks 

In order to verify that the results are stable under vari- 
ation of the vertexing algorithm that is used for the mea- 
surement of At, several less powerful alternatives to the 
default method have been considered: 

• Charmonium mass constraint for vertex fit: 

The mass constraint on the charmonium daughter, 
used in the selection of the events, is also applied 
in the determination of the vertex. 

• No mass constraint: The mass constraint on 
the candidate is not applied during the vertex 
reconstruction. 

• No Bremsstrahlung recovery: Only events 
without an associated Bremsstrahlung photon for 
the J/ij] daughter electrons are considered in the 
likelihood fit. 

• Charmonium daughters only: The vertex of 
the fully reconstructed B meson is reconstructed 
only with the tracks from the charmonium daugh- 
ter. 

• No converted photon veto: Pairs of tracks from 
gamma conversions are retained in the vertex fit. 

• ""At requirement: Only events with a^t < 1-4 ps 
are retained, as is required in the mixing analysis. 

• Boost approximation: The boost approximation 
(Eq. |2^) is used to convert the Az measurement 
into At. 

• Kaon veto: The more restrictive requirement 
from the mixing analysis that no kaons participate 
in the tagging vertex is applied. 
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• No Y(4S ) co nstraint: The algorithm described 
in section VA is simplified by dropping the T(4S') 
momentum constraint. 

• Dilution dependence on er^t- Dilutions for the 
kaon category are parameterized as a function of 
the error on At. 

A summary of the results obtained with these differ- 
ent configurations for the At determination is provided 
in Fig |3^. In all cases the variation of the measured 
asymmetry is consistent with the error assigned to the 
parameterization of the resolution function. 



6. J/iIj K1 background cross checks 

As a cross check, a likelihood fit was performed to 
J/ipK^ candidates in a AE sideband region (20 < AE < 
80 MeV) treated entirely as signal events. This sample is 
actually a mixture of B decay modes with an expected 
average rjcp of +0.04. The true value for sin2/3 in the 
Monte Carlo simulation is 0.7 and consequently the ex- 
pected result from the likelihood fit to the At distribution 
of the control sample is 0.03. The actual fits to side- 
band regions in data and Monte Carlo simulation find 
sin2/3 = 0.16 ± 0.18 and -0.03 ± 0.10 respectively, both 
of which are consistent with expectations. 

As another cross check, a sample of J/ifi if ^ events was 
selected in the data, where only the direction infor- 
mation was used, thereby emulating the K'^ selection. 
The purity and background composition of this control 
sample is very similar to that of the J/ip sample. How- 
ever, in this case, the subset of true J/ipK^ events can 
be identified with the normal J/'0 ifg selection criteria. A 
fit to the AE distribution of the full control sample finds 
(49 ± 3)% signal, which is in good agreement with the 
fraction, 47%, observed for the cleanly identified J/ipK^ 
subset. Likewise, a likelihood fit to the At distribution 
of the full control sample agrees well with the value of 
sin2/3 obtained with the true J/tpK'^ subsample. 



7. Graphical display of the asymmetry 

An elegant display of the CP asymmetry in the data 
can be obtained with the use of the so-called Kin vari- 
able, hereafter denoted as /C. It is also possible to verify 
directly the fitted value for sin2/3 from the ratio of ap- 
propriate weighted averages for /C. In particular, JC has a 
PDF with an asymmetry known to be linearly dependent 
with a slope given by sin2/3 regardless of the details of 
the analysis. 

Writing the PDF for events with a B° {B") 

tag in terms of the general functions Fi{At) and F2{At) 



jn' mass constr. 
no Ks mass constr. 
no brem. recovery 
charmonium only 
no conversion veto 

boost approx 
kaon veto 
no Y(4S) 

co(a^,) 

systematic error 



0.003±0.010 

0.000±0.005 

0.039±0.018 

0.009±0.025 

-0.012±0.009 

-0.014±0.016 

0.013±0.025 

-0.017±0.009 

0.016±0.025 

0.002±0.003 

±0.030 



-0.1 -0.05 0.05 0.1 
5sin2p 

FIG. 32: Results obtained with several alternative variations 
on the vertexing algorithm that impact the At measurement. 
The shaded band is the systematic error assigned to the pa- 
rameterization of the resolution function. The full range cor- 
responds to one statistical standard deviation. 



allows us to introduce 



/C(At) = ±i^2(At)/Fi(At), 



(48) 



where -1- applies to events with a B'^ tag and — with a 
B'^ tag. Ignoring resolution effects, dilutions and back- 
ground, the Standard Model expectation for the At dis- 
tribution of tagged B^ decays into CP modes (Eq. ||) 
gives /C(At) — —rjcp sin Anid At. When these efi^ects are 
included, we can still write 



T±{At) = Fi(At)(l-fA:(At) sin2/3). 



(49) 



although JC will be a more complicated function of At 
and could depend on kinematic variables as well. 
The distribution of events as a function of /C is 



dN 

lie 



dAt 



T+5 JC 



El 
Fi 



El 
Fi 



(l + /Csin2/3) / dAtx 



F^ 



(1 +/Csin2/3)*(/C), 



(50) 



J'±(At) = Fi(At) ± sin2/3F2(At) 



(47) 



where 4'(/C) is an even function of K.. It follows that the 
ratio of the odd to the even part of the distribution for 
/C is a linear function of K, with coefficient sin2/3. Thus, 
the distribution of K. can be used to test for the effect of 
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an CP violation simply by examining the dependence of 
the asymmetry 



A{JC) 



dN,c>o/dlC - dN,c<Q/dlC 
dN,c>o/dlC + dN^<Q/dJC' 



(51) 



The observed asymmetry for the CP sample is shown in 
Fig. ^ as a function of /C, along with an overlay of the 
expected linear dependence. The data agree with this 
hypothesis at the 55% C.L. 



< 




FIG. 33: Observed asymmetry A{IC) as a function of |/Cj, with 
an overlay of the expected linear dependence superimposed. 



From the expression in Eq. we find 



sin2/3 



± 



(sin2/3): 



(52) 



In averaging K. the even component of the ^'(/C) cancels 
out, while the odd component cancels in averaging /C^. 
This offers a method of measuring sin2/3 that is math- 
ematically equivalent to the result with the global like- 
lihood fit. However, it can only be applied when sin2/? 
is the one remaining free parameter. The moments of 
/C for the full CP sample give results that are numer- 
ically identical to the likelihood fit, thereby confirming 
the minimization procedure used for the fit. 

The fact that the mean value of K. is proportional to 
sin2/3 also allows a visual representation of the CP asym- 
metry. Fig. 1^ shows the distribution of JC in data, with 
events in the individual tagging categories indicated as 
well. The larger the value of K, for a given event, the 
larger the weight that this event carries in the measure- 
ment of sin2/3. Again, the CP asymmetry in the data is 
clearly evident in the distribution of /C. 



D. Fits results without assuming |A| = 1 

A more general description of the time evolution of 
neutral B decays to CP eigenstates contains a term pro- 
portional to cosArridAt (Eq. |^). The coefficient of the 
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FIG. 34: Distribution of the observable IC for the individual 
tagging categories. 



cosine term is expected to be negligible in the Standard 
Model, where |A| = 1. In order to search for a non- 
standard Model effect, we fit the rycp = — 1 sample for 
|A| and /mA/|A|. The latter is equal to sin2/3 if |A| = 1. 
The rjcp = — 1 sample has the advantage of having very 
little background, while the other CP modes have back- 
grounds that are both significantly larger and dominated 
by other B decay modes with possible direct CP contri- 
butions. 

The fitted values for |A| and /mA/|A| with the CP = 
— 1 sample and various subsamples are listed in Ta- 
ble KVII. The two CP parameters are almost uncor- 



related, with the coefficient between /toA/|A| and |A| of 
— 1.7%. The sam e system atic error studies as described in 
detail in Section VIII B| were repeated for the fit to the 
r]cp = — 1 sample for /mA/|A| and |A|. The estimated 



uncertainties from these sources are listed in Table XIX 



We have also performed detail ed cros s checks, simi- 

In particular. 



lar to those described in Section VIII C 



large samples of parameterized simulation, as well as full 
Monte Carlo samples, have been used to verify the fit- 
ting procedure. The Sflav sample has also been used to 
demonstrate that no bias is introduced in the measure- 
ment. The relative normalization of the tagged events in 
the two fiavors is in fact sensitive to the coefficient of the 
cosine term in Eq. ^, and therefore |A|. The systematic 
error introduced by the uncertainty on the parameters 
(etag)i and fii listed in Table ^ are uncorrelated between 
tagging categories. Therefore, they are added in quadra- 
ture t o obt ain the systematic error contribution listed in 
Table pax| (g). 

The final result of the fit with the rjcp = — 1 sample is: 



|A| = 0.93 ±0.09 ±0.02 and 
/toA/|A| = 0.56 ±0.15 ±0.05. 



(53) 



Thus, we find no evidence for direct CP violation in the 
ricp = — 1 sample and the value of /mA/|A| is consistent 
with the result from the nominal CP fit with |A| = 1. 
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IX. CONCLUSIONS AND PROSPECTS 

In 29.7 fb^^ of e+e^ annihilation data collected near 
the T(4S') resonance, we have obtained a new measure- 
ment of the time-dependent oscillation frequency 
with a sample of 6350 tagged flavor-eigenstate me- 
son decays that are fully-reconstructed in hadronic final 
states: 

Amd = 0.516 ± 0.016 (stat) ± 0.010 (syst) ps^^ 

This result is at a level of precision comparable to the 
most recent world average for Amd and lies about l.la 
above the combined value of 0.472 ± 0.017 ps~^ [0. It 
is also quite compatible with our own recent measure- 
ment with a dilcpton sample. The Arud study re- 
ported here confirms our understanding of B reconstruc- 
tion, flavor tagging, and At resolution in our data sam- 
ple. Our measurement contributes significantly to the 
precision of the determined value for Auid, one of the 
fundamental parameters constraining our knowledge of 
the CKM matrix, and remains dominated by statistical 
errors that will improve with more data. 

We have presented a measurement of the CP- violating 
asymmetry parameter sin2/3 in the neutral B meson sys- 
tem: 



sin2/3 = 0.59 ± 0.14 (stat) ± 0.05 (syst). 



(54) 



which establishes CP violation in the system at the 
4.1(7 level. This significance is computed from the sum in 
quadrature of the statistical and additive systematic er- 
rors. The probability of obtaining the observed value or 
higher in the absence of CP violation is less than 3x 10~^. 
The corresponding probability for the r]cp = — 1 sample 
alone is 2 x 10"'*. Our measurement is consistent at the 
1.9(7 level with the recently reported result from Belle of 
sin2/3 = 0.99±0.14 (stat) ±0.06 (syst) and with pre- 
vious measurements ||30|| . The observed value for sin2/3 is 
currently limited by the size of the CP sample, allowing 
for substantial improvement as more data is recorded in 
the next few years. 

We have also used the rjcp = — 1 sample to search for 
possible direct CP violation through interference of decay 
amplitudes. The direct CP parameter A is found to be: 

|A| = 0.93 ± 0.09 (stat) ± 0.02 (syst). 

This result is consistent with the Standard Model ex- 
pectation, where |A| = 1 and no significant direct CP 
violation should exist in charmonium decays. 

As already noted in Section I, measurements of CP 
asymmetries in B decays to charmonium can be used 
to constrain, with little theoretical ambiguity, the pa- 
rameters of the CKM matrix. In the Standard Model 
with three families, the CKM matrix V ^ incorporates 
three real parameters and one phase S generating CP 
violation if J 7^ or tt. The Wolfenstein parameteriza- 
tion of V takes advantage of the observed hierarchy 



in the matrix elements in terms of the expansion param- 
eter XcKM = \Vus\- The remaining parameters in this 
representation are denoted A, p, and 77, where CP viola- 
tion requires t] ^ 0. 

The parameter Ackm is determined from semileptonic 
kaon decays and nuclear f3 decays. Semileptonic B meson 
decays to charm are used to determine the parameter A. 
Constraints on p and 77 are obtained from CP violation in 
mixing in the kaon sector \eK\, the ratio |Ki6/Vch|, and 
the oscillation frequency Arud for B^-B^ mixing. The 
oscillation frequency Am., has not been measured, since 
B^-B^g mixing has not been observed yet. However, the 
observed amplitude spectrum A{Ams) improves the con- 
straints on p and rj. Together, these measurements pro- 
vide indirect constraints on sin2/3. 

Our overall knowledge of the CKM parameters is lim- 
ited by the relatively large uncertainties in some of the 
theoretical quantities, mainly due to non-perturbative 
QCD effects. In particular, the constraints on p and rj 
suffer from theoretical and systematic uncertainties in 
the determination of and from theoretical un- 

certainties in QCD parameters entering the prediction of 
\^k\, Anid^ and Am^. Recent analyses constraining the 
CKM matrix have been performed with different statis- 
tical approaches [g, |3|, |8[ |39[ |4^ . They mainly 
differ in the treatment of theoretical uncertainties and 
also in the choice of the input values and their errors. 

Due to the four-fold ambiguity in the value of (3 ob- 
tained from the sin2/3 measurement, there are four al- 
lowed regions in the p-r] plane. One of these regions is 
found to be in agreement with the allowed p-rj region ob- 
tained from CKM fits within the Standard Model. Fig- 
ure taken from Ref. |34 , shows our direct measure- 
ment and the indirect constraints in the p — 77 plane in 
terms of the renormalizcd parameters p = p(l — Ap^jj,^ /2) 
and ^ = 77(1 — Ap^j(j/2). The contributions of the indi- 
vidual measurements le^-l, iKib/Vch], Arud, and Amg HJ] 
are indicated, as well as the allowed region if all the 
constraints are considered simultaneously. Overlaid as 
the diagonally-hatched area are the regions correspond- 
ing to one and two times the one-standard-deviation ex- 
perimental uncertainty on our sin2/3 measurement. 

It should be emphasized that, beyond being a direct 
constraint on j3, the measurement of sin2/3 differs quali- 
tatively in its interpretation from the indirect constraints 
on fi obtained from \eK\, iKib/V'cbl, A'md^ and eventually 
Arrig. For sin2/3, the size of the allowed domain is de- 
termined by well-defined experimental uncertainties that 
are predominantly statistical in origin, while in contrast 
the region allowed by the indirect measurements is mostly 
defined by theoretical uncertainties, which makes a sta- 
tistical interpretation difficult. 

The current experimental uncertainty on sin2/3 has 
now reached a level of precision that offers significant 
constraint on the Standard Model. Over the next few 
years there will continue to be substantial improvements 
in precision of the sm2[3 determination, including mea- 
surements for other final states in which CP-violating 
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FIG. 35: Present indirect constraints on the position of the 
apex of the Unitarity Triangle in the (p, ^) plane, not in- 
cluding our measurement of sin2/3. The fitting procedure is 
described in Ref Q. Our result sin2/3 = 0.59 ± 0.14 (stat) ± 
0.05 (syst) is represented by diagonally hatched regions, corre- 
sponding to one and two statistical standard deviations. The 
individual indirect constraints lie between the pairs of solid 
hues that are connected by the double-ended arrows with la- 
bels. 



asymmetries are proportional to sin2/3. Beyond this, 



studies of time-dependent asymmetries in modes involv- 
ing b ^ u transitions have already begun ]4l|] and may 
provide additional constraints, although here the inter- 
pretation in terms of sin2Q! from the Unitarity Triangle 
is likely to be made difficult due to significant penguin 
contributions. Nevertheless, these measurements will be 
able to directly test the validity of the CKM picture as 
the origin for the observed CP violation in neutral B 
decays. 
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